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1.0 Introduction. 


Prediction methods based on probability density functions (pdfs) have played an important 
role in turbulent combustion for some time. The interaction of mixing and chemical reactions 
is an important factor in determining the performance of practical combustion devices. In 
the special application of the SSME pre-burner, this interaction can affect the internal fluid 
mixing and as a result significantly alter the exit temperature profile. To fully anticipate 
the effect of the interaction of mixing and combustion reactions in a design process requires 
extensive research in turbulent combustion in a range of Mach numbers relevant for the 
applications. The aim of this effort is to develop accurate prediction methods. 

Prediction methods for turbulent reacting flows developed in analogy to noureactiug 
flows were based on statistical moments of first and second order. The mean value of density 
(and other thermodynamic variables) were determined using an assumed form for the pdf of 
the scalar variables describing the local thermodynamic state. This approach is acceptable for 
reacting flows if only the expectations of the stable species and temperature and density are of 
interest and where the reactions are so fast that equilibrium is achieved. It is well known that 
assumed forms of the pdf are not flexible enough to represent truly the variation of the pdf 
occurring in a turbulent nonhomogeneous flow with finite rate chemistry ( Kollmann and C hen. 
1992). Hence, methods employing the pdf directly have been developed for turbulent reacting 
flows (Pope, 1985, Kollmann. 1990). They have several significant advantages, notably the 
ability to deal with the highly nonlinear source terms arising in combustion rigorously in 
closed form. The aim of the present project were the analysis of pdf methods and the 
development of closures for mixing and turbulent transport of single point pdfs in compressible 
and incompressible flows. 

The main parts of the project were the analysis of the foundations of pdf methods 


1 



including the recent development of mapping closures (Chen et ah 19S9). The closure models 
for turbulent mixing were analyzed in detail and pdf methods were extended to compressible 
turbulent flows. The application of a particular pdf method to supersonic turbulent jet flames 
burning hydrogen with air was used as a test case to evaluate the closure model. 

2.0 Objectives. 

The objective of the proposed research project was the analysis of single point closures based 
on pdfs and characteristic functions and the development of a prediction method for the joint 
velocity-scalar pdf in turbulent reacting flows. Turbulent flows of boundary layer type and 
stagnation point flows with and without chemical reactions were be calculated as principal 
applications. Pdf methods for compressible reacting flows were developed and tested in 
comparison with available experimental data. 

3.0 Research work. 

The research work carried in this project was concentrated on the closure of pdf equations 
for incompressible and compressible turbulent flows with and without chemical reactions. 


3.1 Foundations of pdf methods. 

The single point pdf equations, which are the central part of the prediction methods for 
turbulent reacting flows, can be deduced from the exact and closed transport equation for the 
characteristic functional containing all the statistical information on the complete flow field. 
The result of this derivation is the equation for the characteristic function, corresponding via 
Fourier transformation to the single point pdf. It follows that two equivalent formulations 
of single and multi-point pdf equations emerge as a consequence of the functional equation 
at special argument functions composed of Dirac-pseudofunctions. All nonclosed terms can 
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be given in terms of the characteristic function or pelfs and closure models can be set up in 
either formulation. This fact can be exploited to obtain equivalent expressions for a closure 
model. For instance, it turns out that the exact mixing term and, therefore, the pair-exchange 
model for it in the single point pdf equation, has the property to increase the width of the 
characteristic function (analogous to a random process with positive diffusivity for the pdf ) as 
turbulence decays, since the limit of zero fluctuations is given by unity as Fourier transform of 
the Dirac spike for the pdf. The detailed discussion of these results can be found in appendix 

I. 


3.2 Interaction of Turbulence and Chemical Kinetics. 

The objective of this part of the project is to provide a fundamental understanding of the 
physics inherent in various processes causing turbulence to interact with chemical kinetics. 
In view of the importance of various practical combustion processes that occur in turbulent 
flows, the emphasis is put on the influence of turbulence on chemical kinetics. 

The interaction of turbulence and chemical reactions occurs in turbulent reacting flows 
over a wide range of flow conditions. Various degrees of interaction between turbulence and 
chemical reactions can lead to different phenomena. Weak interactions between turbulence 
and chemical reactions may simply modify the flame slightly causing wrinkles of flame smface 
(Williams, 19 S 9 ). Strong interactions could cause a significant modification in both the 
chemical reactions and the turbulence. If chemical reactions cause small density changes in 
the flow, then the turbulence is weakly affected by the chemical process, but the turbulence 
may still have strong influence on the chemical reactions. However, the purpose of combustion 
is generating heat; therefore, one expects large density changes (i.e.. an order of magnitude) 
which can alter the fluid dynamics significantly. 

The research work was aimed at the investigation of mean reaction rates in nonpiem 
ixed systems. Rigorous bounds were established for the mean reaction rates in binaij and 
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multi-component mixtures for given fluctuation levels of composition and temperature. The 
combustion of methane with air was used as an application of the prediction model incor- 
porating finite rate chemistry. Mixedness parameters were evaluated and the bounds on the 
reaction rates were verified. The main conclusion was that, the quasi-laminar calculation of 
mean reaction rates is unacceptable if unconditional mean values are used. The details can 
be found in appendix II. 

3.3 Mapping methods for pdf equations. 

A new approach for the closure of pdf equations was suggested by Kraichnan (Chen et 
ah, 1989) during the grant period and an investigation of the applicability of mapping closures 
to turbulent combustion problems was undertaken. Kraichnan's idea to apply mappings as 
tool in constructing closures for pdf equations (Chen et ah 19S9, Kraichnan 1990, Feng 1991. 
Pope 1991, Valino et ah 1991) proved very successful for the ca r e of a single scalar variable in 
homogeneous turbulence. It was not clear how powerful this approach is for the case of more 
than one variable (Pope, 1991). Pope's (1991) method relies on the cumulative distribution 
function and the representation of the pdf of n > 1 random variables as the product of n 
conditional pdfs. The resulting closure is, therefore, dependent on the ordering of the n 
variables appearing in the conditions. The investigation of mapping methods in the context 
of the present project lead to the important result, that the use of the cumulative distribution 
function can be avoided altogether and that the mapping equations can be established directly 
using the pdf equation. No particular ordering of the variables is required. The detailed 
results are presented in appendix III. 

3.4 Prediction of supersonic turbulent flames. 

The prediction of turbulent supersonic nonpremixed flames was the central part of the 
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present project. The effects of compressibility, the interaction with random shocks created in 
the turbulent zone and shocks created outside the turbulent zone are crucial to the successful 
prediction of compressible turbulent flows. Pdf methods can be adapted to deal with these 
phenomena and a detailed investigation into pdf formulations for compressible reacting flows 
was carried out. It was found that even for equilibrium chemistry at least three scalar variables 
are necessary to fix the local thermodynamic state. In fact, it is advantageous to consider 
a fourth scalar to obtain a pdf equation with the familiar structure. Compressibility effects 
are dealt with using the generalized Langevin approach. A pdf closure including mixture 
fraction, the logarithm of the (dimensionless) density, the internal energy per unit mass and 
the relative rate of volume expansion was established. The supersonic flames of Evans et al. 
(1978) were used as test cases and several successful runs were carried out. The results are 
presnted in detail in appendix IV. 
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Abstract. Probability density function (pdf) methods provide a complete statistical description of 
turbulent flow fields at a single point or a finite number of points. Turbulent convection and 
finite-rate chemistry can be treated in closed and exact form with pdfs in contrast to methods 
based on statistical moments. The equations for pdfs at a finite number of points are indetermi- 
nate due to molecular transport and pressure-gradient terms which require pdfs of higher order. 
The theoretical foundation of pdfs methods are developed in this paper starting from the exact 
and linear equations on the functional level. The closure problem for single-point pdf equations is 
treated in detail and several closure models are analyzed. Turbulent combustion at low Mach 
numbers constitutes an important area of application and selected results for a turbulent methane 
flame are presented as an example. The extension of pdf methods to supersonic turbulent flows 
with and without chemical reactions are outlined. Progress in the numerical solution of pdf 
equations is reviewed briefly. In the concluding remarks, both the advantages and disadvantages 
of pdf methods are evaluated. 


1. Introduction 

Significant progress has been achieved over the last 10 years in the theory and application of 
evolution equations for probability density functions (pdfs) to turbulent flows at low Mach numbers. 
Pope (1985) reviewed the development up to 1985 and provided a detailed introduction to this 
subject. The present paper is concerned with the theoretical foundation and recent development of pdf 
methods. Pdf methods derive their justification from the basic fact that turbulent convection and 
chemical reactions can be dealt with in exact and closed form. This is in stark contrast to the 
approach based on statistical moments, which requires closure models for nonlinear processes such as 
convection or chemical reactions. Pdf methods succeed here because they transform certain nonlinear 
processes into linear terms with variable coefficients by converting the associated dependent variables 
in the basic laws into independent variables of the pdf. Hence, two of the most important closure 
problems encountered in moment equations are overcome by pdf methods. Furthermore, pdfs provide 
a complete statistical description of the fluctuations at a single point or a finite number of points in 
the flow field. However, the equation governing the evolution of the pdf at n points is indeterminate, 
because the terms accounting for molecular transport and the fluctuating pressure gradient require the 
pdf at n + 1 points. The closure problem for these two terms must be overcome to arrive at a 
determinate equation. The numerical solution of pdf equations was for some time considered next to 
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impossible, rendering .he interest in pdfs 

neered by Pope (19851 proved very succ^sful for solution of pdf equations is the 

cular the parabolic type). The reason for he difficulty tn ‘ “"‘"mem me thods. which may consist 

large number of independent variables o t e P in co independent variables. The reason for the 

of a large number of equations governing functions of J umerical ef r ort grows only linearily 

success of stochastic simulation techniques is the fact that the nummca ^ nows ranging 

with the number of independent variables s q and Wu , 987; Haworth and 

from incompressible turbulent shear ayers (see p , ’ n uctuat j ons a nd finite-rate chemistry 

5S ssr-* - - — - 

’’‘“of this paper is ,0 provide a 

methods. In Section 2 the basic laws govern g the exact and linear equations for the 

Eulerian and the Lagrangean frame or ater re ^ are discusset j. They form the theoretical 

characteristic functionals in Eulerian an agrang Eulerian and the Lagrangean frames 

basis for pdf and moment methods. Pdf equations in character istic functions, which follow 

are then derived as Fourier transforms of the equations pdf equatlon is the 

from the exact equations on the functiona ev . possible formulations for the nonclosed 

primary focus of the subsequent sections. and the closure models for 

terms in the pdf equation are discussed firs . eir p P . are rev iewed. Single-point pdf 

the molecular-transport and the fluctuating pressure g scales Hence, methods of incorpo- 

equations do not provide information on turbulen e g h prediction method. Section 4 is 

rating scale information are introduced in or er o areas . The most important 

devoted to the application of pdf methods an * turbulent nonpremixed methane flame 

applications are turbulent combustion ows. methods Then the extension to supersonic flows 

is presented in some detail to verify along with directions for future 

and the interaction of turbulence with shock (1985) review article. 

research. The numberical-solution met o was pr hrieflv in Section 5 Finally, in Section 6 

Hence only the most recent developments are discussed briefly in Section o. y, 

conclusions are drawn for the theory and application of pdf methods. 

2. Theoretical Background 

The analysis of turbuLn. flows Is restricted to Newtonian Jit £ £ SZ7JSZ i°e atmbied 
thermodynamic relations for Ideal gases are assumed to ^ ^hat pdf methods are par, of a 

formulation is discussed briefly and the pdf-transport equation is derived from . . 

2.1. Basic Laws 

The thermodynamic stale of a flowing mixture of ^ZsTrThis^t 

two indecent (mien jl I momentum, and energy and the thermodyna- 
of variables are governed by the balanc . evera i f ra mes* we consider their 

mic state relattons for ideal gases. The balances can be s« up " ” „ the flow held, and 

form in the Euler, an frame, where the flows otoervedat ** poim of the fluid, 

the Lagrangean frame, where the flow ,s observed x an d time ,. whereas 

^ttl x,.,„ and serves as transformer, on between the frames 

x = X(a ,t), a = X 1 (x,t), 

where X- denotes the positton a, time zero of the mater, at point, that is a. position x a. time . 
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(uppercase letters indicate dependent variables in the Lagrangean frame and lowercaselettersareused 
for the Eulerian frame). If the mapping X and its inverse X“ are tw.ee continuously ^erentiaWe 
then the partial derivatives in the Eulerian frame can be transformed into the Lagrangean frame and 
vice versa. The gradient, for instance, transforms according to (Euler relations (Truesdell, )) 

(1) 
( 2 ) 

( 3 ) 


and 


8 

ax a _ 2J Wi, “ 

3__J 

8a, 


dX„ dXy 8 
da , 8a m da s 


dXp 1 dx; 1 e 

8x„ 8x a dxt 


where J denotes the Jacobian determinant 


_ A dXs 8X„ dX„ 

J - 6 W*K> fa* dajj pa. 


(Note that repeated subscripts imply summation and that e a p y is the permutation tensor.) Repeated 
application of (1) or (2) leads to transformation formulae for second and higher derivatives. e wi 
need in particular the relation for the Laplacian, which is given by 

dX 0 d (\ dX- 8X V 8 


1 

A * “ 2J £xpy£s ' ,<a 8a n da a 8a d \j dap 8a y da 


(4) 


The time derivative in the Lagrangean frame plays a fundamental role, because velocity and acceler- 
ation are by definition given as 

(dX, 

' 8t 

It transforms to the Eulerian frame according to 


t) = , . 


(dV \ 


8 D 

+ ^ (X ’ () Z s m 


(5) 


'8 

Kpl}* _ \8t 

which is called the substantial or Stokes derivative. The transformation rules (1) (5) enable us now to 
set up the basic laws in both frames. 


Mass Balance. Mass is conserved and this statement translates into 


8t 8x x 


( 6 ) 


for the Eulerian frame, where p(x, t) denotes the density. Transformation to the Lagrangean frame 
leads to an integral of (6) given by 

R(a,0) 


R( a, 0 


= J, 


(7) 


where R(a, t) = p(x, t) for x = X(a, t). Equation (7) is therefore the mass balance in the Lagrangean 
frame. 

Species Balance. A mixture of n ideal gases is considered and its composition is described in terms of 
mass fractions Y t ( a, t ) = y f (x, t), x = X(a, f) i = 1, ... . n. Chemical reactions may occur and the i Y, are 
therefore not conserved, but may be consumed or produced according to a reaction mechanism 
consisting of many steps. At this point we only need to know that the rate Q, of production due to 
chemistry is a local function of the thermodynamic variables (no derivatives or integrals with respe 
to time or space/label appear in the &). The balance for the mass fraction y t in the Eulerian frame 

then given by 

Dy, d 

P 


Dt 8x„ 


png ] + 


i = 1 , .... n. 


(8) 


w. Kollmann 


252 

where r, denotes the Fickian mass dilTusivtty. Transformation 10 the Lagrangean frame is earned out 
using (1), (5), and (7) and results in 

8Y,_ 1 + j = 1, .... n, (9) 

IT ' da„ da „ da\ R 0 da fi da y daj 

where R 0 = 0)* 

Momentum Balance. Newton's second law leads to the balance equation for momentum. In the 
Eulerian frame it appears as 

Dv. dp , dx„ A (10) 


Dt 


= + 


Zfi 


dx n dx B 


+ Pfcn 


where is the stress tensor and /, is the external force per unit mas. Transformation to the 
Lagrangean frame is straightforward and results in 

(U) 


8K 

dt 


1 


ItC*'*** da, 8a m da, + 2 R 0 da , da a da, 


dX I dX 1 dT 3£ + F 

£ “' A ™ " !, da a da. 


dXp dX y dP 1 _ C*y c?A c CI *l 


, t — t (x is the stress tensor in the Lagrangean 

where P{ a, t) = p(x, f) is the pressure and T afi ( a, 0 t /,(x, j 

frame. Newtonian fluids satisfy the linear constitutive relation 


dv, dv 8 


dv. 


T *' = /i U^ + S 3 4 "« y 


( 12 ) 


between stress and rate of strain where p is the dynamic viscosity. 

Energy Balance an d Stale Relation, The firs, law of “ b “ "'Tp tn “«'»! 

ssrss - - 

given by 


DH D P | rj, 

<’d7 = d7 + ®‘H: 


(13) 


The specific enthalpy for a mrxlure of ideal gases is composed of the enthalpies of the components 


n U . 


(14) 


where h is the molal enthalpy and M, is the molecular mass of the ith component. The molal 
enthalpy h t in turn consists of the formation enthalpy (if and the sensible enthalpy 


hi = hf + r 

J r 0 


dT c p (T'), 


(15) 


where c p (T) denotes the molal specific heat at constant pressure. The dissrpation function * is defined 

by Sc. 


<D = T 


*p 


dx p 

and the energy flux q. consists of conductive, diffusive, and radiative fluxes 

. dT v ^Lr^i + / 7 B 
q '~ k dx . P J, M, T, dx. 

The system of equations is closed if the ideal gas equation 

p = pRT 


(16) 


(17) 


(IB) 


is included and the chemical sources fl, are spec, lied. 

be deduced from the first law of thermodynamics for a differential system or by 


The pdf Approach to Turbulent Flow 


253 


Eulerian form (13). The result is given by 

dt 

where 


1 l dXedXydQt 

Rdt ' ' 


_ ! dP 


da,„ da* 


I(D = 

R 


1 

2 T 0 


^a^y^drtvt ■ 


dX. dX y dV p 
’ da„ da,., da* 


(19) 


( 20 ) 


wmmmiwm 

igssss mm 

equations (Eulerian frame) 


D<P, 

1 Dt 


d 

dx„ 


p r £ 


# i 

dx„ 


+ pq b 


= l 


or (Lagrangean frame) 
dT: 1 


dX r dX a d (R 2 T ( 


dt 


2R 0 da„ da M da* \ R 0 da p da y da 


dX* dX„ d'V 


d 


+ Q h 


( 21 ) 


( 22 ) 


, ,\_nu n « = Xla fl are not identical with the sources in (8) and (9). Their 

with in a slightly different fashion. 

2.2. Characteristic Functionals 

A complete statistical description of a turbulent (low can be achieved if f 

(see Hopf. 1952; Hopf and Tilt, 1953; Lewis and Kratchnan. 1962; Foias, 1974, Vtshik « »!., 1919, 
Constantin et at . , 1985) 


m[d , v, </>,] = <^exp ji if) + (v, «) + £ M. %)] 

for variables in the Eulerian frame or 


M[d, x, <ft] = ( ex P y 


^rTfR, d) + (X, x)+ t (T f , <p,) | 


(23) 


(24) 


the flow domain fR x [0, T], The expressions (p, d), etc., indicate the scalar produc 


in the Eulerian frame and 


(t>„, o,) -- 
(X x , x.) = 


dxvjx, t)v a (x, t ) 


m 


daXJa, t)xj a, t ) 


(25) 

(26) 
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in the Lagrangean frame, where 9?(£) is the flow domain at time t and [0, T] is the time interval of 
interest. The angular brackets represent the mathematical expectation, which is defined as a functional 
integral over the space of all realizations of the turbulent flow field (see Vishik et al. (1979) for the 
definitions of the appropriate function spaces). This functional integration requires the existence of the 
probability measure (see Daletskii (1962) and Skorohod (1974) for functional integration). For incom- 
pressible flows and homogeneous boundary conditions, the existence of the probability measure has 
been established (see Hopf and Titt, 1953; Foias, 1974; Vishik et al., 1979), but not its uniqueness. 
For compressible and reacting flows this is still an open question. We assume in the following that the 
probability measure for compressible and reacting flows exists. With this assumption we can proceed 
to set up the transport equations for characteristic functionals following the method put forward by 
Vishik et al. (1979). Noting that the characteristic functionals are independent of time t and location x 
and label a, we have to form functional derivatives (see Averbukh and Smolyanov, 1962) in order to 
establish the dynamical change at a given time and location or label. The transport equations for the 
derivatives of the Eulerian functional m can be obtained without difficulty (see Lewis and Kraichnan, 
1962; Dopazo and O’Brien, 1974; Kollmann, 1987). The mass balance (6) leads to 

d dm _ . d d 2 m ^ 7 ) 

dt Sd(\, t) ‘ dx a 6d(x, t)dv a (x, t ) ’ 

where d/dd{x, f) denotes the functional derivative (Averbukh and Smolyanov, 1962) defined by 


J* 


dx 


T dt TTT\ h < x > '> = -j m L d + eh] I.- 0 - 
0 dd{x, t) de. 


(28) 


Momentum balance (10) leads to 


8 l m 


dt Sd(x , t)Sv a (x , r) l dx p 8d(x , *)&>«(** t)So fi (x, t) ^ dx J 8x p ^ ' u 8d(x, t)' 


8YI 8 T a 

_l_ 




+ tfat 


dm 


(29) 


where the external force per unit mass f a was assumed to be nonrandom and n denotes the pressure 
functional 

' T 

(30) 


■r 


n[d, <p u <p{] = ( p(x, f) exp<{/ | dx[(p, d) ■■•] 
and T aP denotes the stress functional 

T afi [d,o, s (z t0 (x,t) exp \i | dx[(p,d) + 


'i: 


The stress functional for Newtonian fluids can be given explicitly in terms of m\ 


71* = 


- iu 


8 


8 m 


8m 


-$v 


8 m 


dx p <3o a (x, 0 8x a 8v p (x , t) 3 dx y 8v y (x, t)J 


(31) 


(32) 


The viscosity \i was assumed to be constant. The thermochemical variables governed by (21) lead 
to (Kollmann, 1987) 


8 2 m 


8 


8hn 


dt 8d(x , t)8(pj(x , 0 


8x p 8d(x, t)8vp(x , t)8(Pj(x, t) 

8 8 


+ i- 


cx* 


J fi 


8 m 




8m 


cx„ S<Pj(x , t) / 

7=1,...,/, 


(33) 


^/^(/>i(x, t) ’ iS<pi(x, t)J Sd(x, t) ’ 

where the transport coefficient pTj was assumed constant. The sources cpj appearing in (21) 

become operators acting on m. If <jj is not a polynomial in the <pj, then it may turn out to be a 
pseudodifferential operator (see, for instance, Taylor, 1974) on the functional level in (33). The 
functional equations for the Lagrangean functional M are obtained in similar fashion as shown by 
Monin (1962) for the case of incompressible turbulence. Mass balance (7) leads to the equation 

S d S d d SM 


SM i ,. d 

,m 


da' a (3x 4 (a', t) da'/, dx n (a", t) da* <5x w (a*, t) Sd( a, f) ’ 


( 34 ) 
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where the limit is carried out for a', a", a* approaching a. The mass balance (7) is not an evolution 
equation and thus (34) is not of evolution type. The balance of momentum (11) leads to 

8 2 8 2 M _ in e , im A ^ 

dr 1 8d( a, 0)<>x 7 (a, f) ~~ ~~ 2 fc «*v e *»<° im s a ' n Sx f { a', f) da" <5x.,(a", t) 8a s 


+ 2 e p-/Z E dqto 


. e 

lim — - 


8T„ 


C0 


da' <5x y (a', t) da" d.x ; (a", t) da 6 


+ if ; 


SM 

<M(a, 0) 


(35) 


and the thermochemical balances (22) imply that 

d 8 2 M _ i n e ,. m _d_ d d d 

dt dd(a, 0|d^(a,7) “ ‘ m da; dx ; (a\ t) da; dx^a", t) 

8 8.8 8 8 8 8 dM \ 

* d^ \ J hm do| dx ; (a*,7) do** dx,(a**, t) 8a. 8<pj( a, 0/ 

+ i q/ - j = 1, (36) 

^ J \id<Pi(a, r)’ ’ i8<pi(a, t)J 8d{ a, 0) 

The transport coefficients Dj = R 2 V } IR 0 are again assumed constant in order to avoid unnecessary 
complications. Variable transport coefficients can be dealt with, but the resulting equations become 
rather unwieldy. 

The transport equations for the Eulerian and the Lagrangean functionals exhibit several properties 
of fundamental importance. First we note that the equations are linear in contrast to the physical 
balances, which are highly nonlinear. Furthermore, the system of functional equations is closed in 
both Eulerian and Lagrangean frames if we consider the first functional derivatives of m or M as 
unknowns. This follows from the fact that both the stress functional according to (32) and the 
pressure functional II can be expressed in terms of m or M. No explicit form for the pressure 
functional can be given, but the thermochemical relations imply that the pressure can be expressed in 
terms of p and i/q, ..., \jf, and thus there exists a relation between FI and m or M. Finally, we outline 
a procedure to establish the characteristic functional m or M from the solutions of the functional 
equations, which provide the first derivatives of m or M. We note that functional differentiation and 
appropriate combination of the solutions of the system of functional equations (27) (33) or (34)-(36) 
lead to a Poisson-Levy equation (see Feller, 1986) for m or M. The Dirichlet problem for this 
Poisson Levy equation can be solved analytically and the result shows that m and M can be 
represented as functional integrals with respect to a Wiener measure (see Theorem 3.5 in Feller 
(1986)). 

2.3. Finite-Dimensional Characteristic Functions and pdfs: Eulerian Frame 

The functional equations contain all the statistical information on the turbulent reacting or nonreacting 
flow. In particular, the transport equations for finite-dimensional characteristic functions follow from 
them. The governing equations for finite-dimensional pdfs are thus determined also, because pdfs are 
the Fourier transforms of characteristic functions. The derivation of the equations for finite-dimensional 
characteristic functions are outlined for the single-point case. We note first that the generalized 
argument functions 

d* = A 8(x - x°)8(t - t °), o? = M(x - x°)<5(r - t°), q>* = <Pj8(x - x°)8(t - t°), (37) 

where d, v x , <pj are parameters independent of x and t, produce the single-point characteristic function 
m l when applied to m, 

m[d*,u* = m i(d,u,<pi,...,^;x°,t°), (38) 


and m x is, according to (23), defined by 

m,(<J, v, <pi,...,<p ( ;x°, t°) = / exp <! i 


<?p(x°, t°) + (3 • v(x°, f°) + X <pjipj(x°, t°) 


j=i 


(39) 


Variational derivatives reduce to partial derivatives with respect to the parameters <?, u, < p p as for 
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instance in 


^-g- id*, u*, vt, . . . , (pn = d ' * • • - ’ x °’ f0) ' 


(40) 


Sd{x°, t°) 

Derivatives with respect to x and t appear differently for m and m. From the definition of the 
functional derivative (28) it follows that 

l 


dtdd(x,t) \dt 


T dx (p, d) + (v, ») + Z Wj, <Pj) 

o L J_1 


whereas it follows from (38) that 


dm l ^/dp_ o o 

£ t°)cxp{»‘[-]}). 


g(0 = id {jp (X°, t° ) exp { i [■"]}) + id, ^ ( x °> f ° ) ex P { 1 f " ^ 


+ i 


>=i 


at 0 


(41) 


(42) 


where 


[...]* <Jp(x°, t°) + 6 • v(x°, t°) + t <°) 


(43) 


holds. Using (40)— (42) it is easy to show that the following relation holds for time derivatives: 

& - + + i w - 

whpre r*l = U* o* <pf <Pi*l Note that differentiation has to be carried out first and then the 
arguments are a^pikd. ^.thermo, e, note ,ha, standard and functional different,*, on do no, 
commute. Similarly it follows for the substantial derivative that 


dm 

mot W 


d 3m, d d 1 m l _ * 

dt 5 ^ + d^JdSm t - mx°, ty i6o a (x°, t°) 

+ iv 


S S 2 m d Pm U*] 

dt° idd(x°, t°)iSv fl (x ( \t <T ) dx° iSd(x°, t°)idvjx 0 , t°)iSvyx°, t ) 


d 2 m 


d 


d 3 m 


+ i E iSd(x°, tySQjix", t u ) + dx° x iSd(x°, t°)iSu,(x°, t°)id<Pj(x°, t°) 


[*] 

(44) 


(superscripts of x°, t°, and hats are omitted from now on) 

d \ 3m, 


(» 

dt idd dX'iddido, M J '’ 


}d<p i 


' id (ft/ idd 


/ , 3d, a , 

= — id ( p 2 - — Wi ) + it > 0 

Fa, 


dp 8x 
dx a dx p 


\ \ i I 8 ( 2i fa 


ri, , (45) 


where 


mj = exp < i 


/ 

dp(x, 0 + rv(x, r) + £ 
j=i 


The transport equation for the pdf 

f 1 (d,v,<p 1 ,...,<Pi-,x,t)ddd\>d<p 1 ...dq> l 

= Prob {d < p(x,t)<d + dd, u, < u a (x, t)<o a + dv x , cpj < i /'/x, t) < q>j + d(p p a = 1 , . . 


3, ;=1 /} 

(46) 
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can be obtained by Fourier transformation of (45) or it can be derived directly from the basic laws 
(see Lundgren, 1967; Pope, 1985; Kollmann, 1987). It can be given in the follow, ng form. 


d^ + do,p-+ I 

"dXa j = l 0 (Pj 


dt 


d / 2 dv * A , A- 

~ dd\ dxj) do a 


dp 

8x x 


dig? 

dx R 


- p/.y 


' 8 / d ( r ^h\f 


where 


t 

f = S(p{x, t) - d)S(y(x, t) - o) d{^j(x, t ) - (Pj 

j~ 1 


(47) 


(48) 


r r.f m The eauations for higher-dimensional (multipoint) characteristic 

is the Fourier transform of m,. The equations lor mg* 

functions and pdfs can be obtained in the same manner see un ^ severa l important properties. 

The equations tor the character^ function », and the pdf f, P arac ,J s , P (unc . 

They were deduced from the closed and linear equ Viscous and diffusive terms 

tional, but they themselves are no, dosed at any jSns or pdfs into the 

as well as the pressure-gradjent mt^^ the P h her . dimensional equations. Hence, the 

S 23 ZZL. and pi. are a, ways indeterminate. 

2 . 4 . Finite-Dimensional Characteristic Functions and pdfs: Lagrangean Frame 

The transport' equations fo, finite-dim.nstonal chastens.* hmctions and pdfs in £ 
rer"d“^ single-point case. The 

generalized argument functions analogous to the Euler, an case would be 

e - in. - «»)4(r). xj = MO - *“)«' - '"»• if = ^ (a - ” 0|il ' “ ' 

where the parameters l x, « are independent of the ^ temsTf' iteMW 

taken a, the reference time a. the generalized 

density via the tn^tatrf ntmis bala 1 characteristic function M,(AM„ ...,*>• Jt » 

arguments (49) turns out to De me a g v . . information on the statistics of the 

important to notice at this point that M, does not contain the information ^ ^ statis _ 

“■ •» mod ' rj ,hs 

argument functions x* as follows. 


x** = x 3 t5(a - a °)6{t - t °) + - a °)^ f “ ^ 


_d_ 

'a? 


(50) 


and regard • d. as a parameter like ^ 

characteristic functional M implies the choice of the derivative dx./rt. q at a sing , e point as 

derivatives become unknown variables. Taking 
the characteristic functional of the modified arguments now leads to 

(51) 


<ta(a°, t°) + x a X a (a°, t°) + f<>) + £ ^ (a °’ f<>) _ 


M[d*,x**, <pf, = \^ ex P 

which is indeed the desired characteristic function M , (l 

and thermochemicaf scalars in fcUp-J- fcrV » “ !ZL with quantities a, 
S°o"n ::"d^d^"Sic e v^ « ” <X 0 ' 
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whereas the Lagrangean function M x depends on the augmented probabilistic set (<?, x, o, <Pj, (pi) 
and parametrically on (a, f). The time rate of change of the characteristic function M, can be 
expressed in terms of the functional derivatives of M[d, x, <p,, </>,] 


d dM x 


d 2 M, 


d 2 


,X ’ lU °dt 02 idd{a°,0)iSx x {a°, t 


s 2 m r _. . ^ A a 

o-^W + ' I ^7 


a 2 M 


at 0 iSd( a 0 , 0)i<5^(a°, t 


(52) 


iddidv 3 Vi louya , ,i / ;=i ^ 

where [*] denotes the arguments given by (51). The dynamic equations (35) and (36) for the functional 
derivatives of M can now be used to eliminate the time derivatives on the right-hand side of (52). The 
result can be given in the form (superscripts and hats are omitted from now on) 

a \ aM, 

K id<p l ' ' idtpj idd 

a d 


8 8M l . a 2 M, ' . 

dt idd tx ‘iddido„ ' ^ (PjQj \ - 


= 


‘ 2 Hm 


+ ?£/>rrf-Sn*, Hm 


a 


an 


aa' Sx p (a\ t) da" <5x.,( a", t) da d 


[*] 


a 


dZ 


a/l | 


^ (Pj£gt0y^&l1<o 
1=1 


lim 


<5x.,(a\ f) da" a <5x ? (a", t) (3a a 

a a a s a 

)fl' <5x ; (a', f) dal ax<s(a", /) Sa 4 

3M 


W 


aM, 


/aj 




. a a a <> o 

x ^1 *' m 5 fl * t) 5. x< j(a**, f) aa. <5(p,(a, f) 


[*], 


(53) 


where the limits are carried out for the labels a', a", a*, a** approaching a. The right-hand side of (53) 
contains the nonclosed terms generated by pressure gradient and viscous transport. The transport 
equation for M, can be considerably shortened (but not simplified) if a mixed Eulerian/Lagrangean 
form of the basic laws is employed, which is given by 


R n 


8K 

dt 


= -J 


dP 
dX \ 


+ J^f+R 0 F X 

( A a 


and 


= RoQj + J 


8X n 


3% 

Rr >e£ 


It follows that (53) can be recast as 


d dM t d 2 M { 

- ix. 


8t idd 


' iddido„ 


1 

i Z <PjQj 

j=i 


id<P i 


dP 


dZ 


if} 


J dX, +J dX, 


M 


d \ cM, 

id(pj idd 

_ SM X . / 8 

' W ° Fa 1dd + ' (f>J \ dX~ 


dVj 


RT >Y£) M ) 


,\ 


(54) 


where 


M = exp i < K(a, 0) J + X a (a, t)x a + F a (a, + X 

7=1 


and therefore 

M x = <M> 

holds. The transport equation for the Lagrangean pdf F x can be deduced from (54) using a Fourier 
transform. The result can be given in the form 

dt dx a d<Pj 


d_ 

8o„ 


dx„ dX„ 


„SF L _±d_._ 

F F *dv a jS d(Pj{ J dX . 


C?T:\ 


( 55 ) 
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where F is the Fourier transfer™ of »j and JeTSSf 

:re^^ix:rn^v^ = - in - szte 

rr r ~L^i« «,.,, > - . ,« 

density at the reference time zero. If we define 

d* = ^(a - a 0 )^(t - ,0 ) 

instead, we obtain the pdf Ff(d, x, «, <P» ^ f ) governed by 

d ^i + dVa d A + £ 
dt 8x* M d( Pi 


d 

Id 


R 


,8V. 

d> (, 


F )- 


d 

dv. 


dP 

'lx* 


+ 


eit, 

dX, 


dF a 


8F, 


‘ diK 


l ^ / 

j?i S(pj\dX a 


<!)' 


(56) 


Ctf \ VA X / ^Ct \\ ~ R, , 

o thp Fnlerian case (47) The relation between the Eulerian 
which has exactly the same structure as in the Eulerian t ) {) of (56) i s 

wiiil<u ii« j aeraneean solution r^ia, x, i), <p x , \ ' 

ScSd;^by^'098S, and Kollntann and “ Thi“s 

^nations and the Mowing sections are now 

devoted to the construction of closed pdf equations. 

3. Pdf Methods 

the prevtous section if was 

mechanical and thermodynamic state . and linear equation for the characteristic 

points in space and time) follows from the de tte rn l»«r«.r. iodeKr- 

functional. The equations for finite-dimensional p solvable set of equations. In 

as wen as ore, hods of 

closure, are discussed. 


3 1 Properties of the Single-Point pdf Equation 

The pdf „;x,t> is governed b, the transport e q ua.,on ,47,. whtch has the ton. of a 

conservation law a 1+4 $ _ ^ 


i(4M + /-(*./,)+ I 

dt ex,, k=\(yk 


()t CX a k = l ' J k 

»i m \ renresents the set of probabilistic and independent variables while F k 

where y k = {d, t), Vl , .... <P f ) rcpvcscnU s t P uatlon must satisfy two fundamental 

represents the corresponding fluxes. The solution oi mis c M 

conditions, Ul 


/i>0 


(58) 


and 


jdy/ l (y;x.r)= 1. <59) 

in order to qualify as a pdf. Dividing (471 b, i and integrating over the probabilistic variables y leads 


to 


Hence, (59) 


8_ 

dt 


dy A = o. 


is preserved provided that initial and boundary conditions conform with (58) and (59). 
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Multiplication of (47) with expressions of the form 


1 + 4 


n (* - <*»■ 


k = 1 


, r i , _ i 1 /-U 4) leads to the transport 

(a k > 0 is an integer) and integration over the range of the y k , k - 1, . - • A + h 

The latter group ( 


j-i 


dcpj 


is closed “ l ° n J : of combusting (lows) and can 

rTo f?:s^ 

— s- — 

co r ^ ^ “ 


d / , dv x A o 

aVWrmv <*.. 

represents the etTec, of volume expansion on the pdf. This term is aero if dens,., ,s constant or 

material points, because 

Dp 


„ dV„ 


: d, V = W, = <jPj, j - U • • * * 0 / 


(60) 


Dt 


6T a 

dx a 


This ease includes stratified Hows (where p ma, „f volume expansion ,s 

incompress, hie (lows (where p is cons, an. in space a*, ''“f J 1 ^” ”d supersonic speeds. The 
nonzero in turbulent reacting (lows and m turbulen lhe dens,., axis, 

conditional moment in (60) can then be viewe a ed and the pdf is moved toward 

now and for points x far 

away from boundaries. Pressure is then determ.ned by a Poisson equatio 

dv x <h) p 

Ap= - p ^eZ 


whose solution 


p(x 


p f dy dv a 8vp 




(61) 


™ j ... yl 8y f 8y a 

allows us to express <(dp/dx.)/> a, follows (see Hanjalic and Launder, 1972; Lumley, 1978; Pope. 


1985): 


dp A _ _A 

dxj 1 4n 


dy ^-\x - yr 1 F ^ 7 <t>(y)My)/>- 
J dx x dypdy, 


( 62 ) 


\ a / 

The expectation in the integrand involves the two-point pdf / 2 (..., *, y, 0. s,nce 

t)r,(y, .)/> « | dMWVM. v.V,. ■■■■*. *• * 

where the values of v(y, t) are denoted b, »■’>. From the fact that pressure is related via a differential 



261 


The pdf Approach to Turbulent Flow 

equation to velocity. follows .ha, the pdf flux due to the pressure gradient cannot be expressed in 

« cases: ftrst, the nonHuctuating case where 

appears in Cosed fonn, and, second, the "up attng case whe, the £ 

SSS ;"o a^Wiener procesri'then ^Sn^ tl^utn balance .us, be wrttten as a 
system of stochastic differential equations conta.ning 

/* dt = (Up dWp{t), 

. AW (U denotes the velocity increment and dW f (t) denotes the increment of a vector Wiener 

^ss^eTeLri W wilh aero Lap and uni, variance. The pdf equatton then con, a, ns 

where is the Sr (47> repiesenc 

^"^“i ,y and scaUr spa. They can be — in ^h a^way that then 


dv. 


dx„ / J 1 C% 


e_ 

dx, 




dx, 


r iVi 

p j dxJ dx f 

dp_ <lf\ 
dxp dd/ 


( ^ T \ ^ 

STte'-Vto' 



dv a cv y 


<«>„/>- 


JL- 

co.dd 


dp 

'dx. 


^010 A n 


’ a *dx a J / Mk=iC(pjC<p k 




v - („r gjjjg-A 

jTj d(pjdd \ dxp / 


- y — - — 

j=i d(pjdv x 


\ CXkOX 


(63) 


where 


®«V = /' 


gti. dv, 
dx p dxp 


and 


dV, dT k 
® ik ~ ^ dx,, dx t 


denote the rates of dissipation. The essenttal 7 1 f rU, *^ 1 ^^ I SS,r transport affects the 

in the underlined expressions. An inspect™, of those terms tsh ows tha ^ mrdccula P 

pdf /, differently in phystcal and veloctty-scaUr spac^The dtsmtatton of P^J^ di|r „ slv£ 

is smoothed exactly the same way as ve y fl at high Reynolds and Peclet 

transport. This contribution is, however, neghgtble fo J ^ |eadjng ord£r 

numbers. The underlined terms representing ranspo h time-inverse heat conduction 

in high Re/Pe-numbers flows and affect the pdf analogously to the time 

equation (i.e., heat conduction with negative diffusely) given by 


dh 

dt 


dv„dv, 


(D x pf y ), 


D a p = <4>.,l P = d,y = v, = <pj> 


where 
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is positive definite (note that the trace of is the dissipation function). The initial-value problem for 
this equation is not well-posed and it follows that closure models for the molecular transport 
cannot be based on this type of differential equation. The effect of the molecular transport terms m 
velocity -scalar space is in essence to reduce variances and convariances while leaving n °™^^ atlon 
and mean values unchanged. The pdf in freely decaying turbulent flows will approach Dirac- 
pseudofunction due to the molecular transport terms acting in velocity scalar space. 

3.2. Formulations 

The presence of nonclosed terms in the transport equations for finite-dimensional pdfs requires 

additional informa, ion in order to arrive at a finite and determinate system * Sna Th^am 
problem can be tackled in several different but not necessarily equivalent formulations. Three ap 
proaches are outlined, two of which have been the basis for successful modeling efforts. 

Formulation F Pdf. The nonclosed terms are analyzed as fluxes of the pdf in the multidimensional 
space spanned by velocity and scalar variables. The exact form of the nonclosed terms is given 
Tectly by the transport equation (47). The essential step in the analysis is the determination of the 
structure of the set formed by all realizable states in velocity and scalar spaces. The velocity space is 
usually 1 *! three-dimensional Euclidean space, but the scalar space can possess intricate boundaries. In 
particular chemically reacting flows lead to a description of the local state that requires many scalar 
va abts The range * of these variables is bounded and the boundary 391 of this range is determined 
^com l«?,3nships (see, for instance, the case of CH. combustion (Chen « al 989),. How =r 
the following fundamental restrictions are imposed on the possible forms of the boundary cX of the 
scalar range- The boundary cbR is an orientable, singly connected, and piecewise smooth hypersur a e 
(dimension / - 1) that encloses a convex body of nonzero /-dimensional volume. The enclosed volume 
does not have to be bounded. The condition of convexity is relevant for the m.xing models 
discussed later The existence of boundaries in velocity or scalar space imposes a condition on mode 
nondosed fluxes: the flux component normal to the boundary (where the norma 
exists) must be negative or zero if the normal is defined as positive outward (which is only possib 
orientable surfaces). Then 

X n k F k (y) < 0, y, on^Jt. (64) 

Jc = l 

If this condition is violated, unphysical states becomes accessible. 

Formulation II: Characteristic Function. Characteristic functions are the Fourier transforms of pdfs 
and therefore an equivalent formulation under mild conditions. They have een use rare y 
treatment of turbulent flows (see Kollmann, 1987). Characteristic funetions possess several - < ng 

properties which can be exploited for the construction of closure models. For the sake of simplicity, 
the case of a single probabilistic variable is considered. The distribution function My) and character- 
L,ic“unebon „,(C) can be decomposed uniquely into the sum of three dist.net eon.r, buttons (Jordan 
and Lebesgue theorems (Lukacs, 1970)) 

FAy) = *F d (y) + bF„(y) + cF s (y), a, b, c > 0, a + h + c = E 

and , /r . 

m,( 0 = am d ( 0 + bm ac (<>) + 

where F is a step function F is absolutely continuous and has a derivative everywhere, and F* is 
continuous but singular in the sense that its derivative is zero nearly everywhere. The decomposition 
of the characteristic function is completely analogous to the distribution function; m d is the charact 
istic function of the discrete distribution, hence it is almost periodic (see Lukacs, 1970) 


lim sup |m d | 

lei - 


i, 


m ac is the characteristic function of the absolutely continuous distribution, hence 


lim m ac (0 = 0, 

id-* 
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M C delhrot 1974' Levich et al 1984). This analysis is based on the fundamental but unproven 

i- Chapter 8 ot Monin and Yaplom „97 5 „. that the rate ot 

d,ssipa “ on n dv, &, 

~ \Re dx a dx # , 

becomes independent ot the Reynolds number for Re » 1 and approaches ? 
value as Re -*• x. The rate of dissipation is dominated by the enstrophy « « 

here as constant) n2 

<e> = <Re' 1 0 2 > + T ^-<Re _1 

w x dx.dxf 

because the dtfference with <r> is only a transport term. Hence the instantaneous rate of ^paho" 

„ nfr ot pH nn subsets of the flow field, where enstrophy becomes unbounded as Ke 

such a random variable. For instance, the discrete pdf 

/j(e; Re) s f 1 — 7-W) + ^ ^ £ “ Re < £ »’ 0 < < £ > < 001 


or the discrete characteristic function 

m,(C; Re) s 1 -T(i -exp{iC<E> Re}) 

^ ''^t'onaf to Re ^ variable "iT pToTccVthat assumes the value zero 

■■■■■I 

Formulation III : Langevin Equation. Single-point pdf equations can be simulatedundercerta.n 

m a, m,«d 

Eulerian/Lagrangean formulation: 

(65) 


( 66 ) 


dX a _ 

v„ 

dt 

d R _ 

K 1 W 

£ * 
QC 
1 

aT - 


i ep 

dt 

~RdX x 


vr„ i ", 1 L'T 7t ! F 

— = ~=^ r + RdX lt + *' 


( 67 ) 
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= ( D ^l 

dt R8X A J SX. 


+ Q P j= « — 


( 68 ) 


The right-hand sides of the baste laws (66) (68) are split into expectations and llttc.ua.ions as follows: 


^-T(R,<)'.» + b(r,||; 

pi / / SP‘ 8T; ty 

~ = A,(R, <P>, <T /)V >, <F.» + B x yR, 1 

s. i - -ii ^ 




7=1,.-,/. 


T he terms denoted with A A A, represent parts" which depend on 

and the dependent variables R, V 7 , %■ ’ *’ 1 terms of the dependent variables at a 

rr:“o^^ — * ^ 

T . Hence, the deterministic parts are given by 

, <?< K> 


A = -R 


dX„ 




R 3X 

1 8 


R 8X„ 


A =__iD:^) + e j (i'i 

T >-R8xA 1 8X„ 1 J 


contain ^(gSTE^ ^ - * 

dependent variables; hence 


B = - R 


8 v: 


B„ = 


8X a 

1 dP' 
P d X 


+ 


Bj ~ R 8X 


R 8X fi 

>*A 'dX.) 


18T - + k, 


Note that the random ,983) (65H«> 

dP'fcXx) coloured noise (see Soong, 1973, Linaenoerg e 

can be written in the form (Wk 

^ = 4 J .(Y,<Y» + B,(Y,(Y» 7r , 

8t 01 

, v = <R X V V t/t.) is the vector of the dependent variables and W k is a random different!- 
Tble vector representing ail the additional unknowns contai ined in ^ ^nnXd The 

and appears as ^ , A) d , + Bf dW t . j- I V + U (70 » 

This system is solvable if the initial condttions and the pr^re*^ r^ct£fhal°been determined. 

,n OTd “ ,o eslab “ 
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The single-point pdf F l (d, x, v, <*>, <Pu ^ 0 = F iO*’ • • - >' +7 ’ a ’ l} Sdl 

(71) 


Fi(y; a, f) = 


dy'Ff(y'; a, t|y; a, r + Ar)F,(y'; a, r). 


where fj deno.es .he pdf of V a, (». , + A,) condih oned upon V -£« (*• '>- The conditional pdf Ff 
can be expressed in terms of the conditional characteristic function M, 


Ff ( y; a, t + At) = 


27T 


dy' exp( — iyiy a )^i(y; l \y ; a? 1 + 


the derivatives of the characteristic function at the origin, we obtain (L = I + /) 

F 1 ,y;.. + A.,= £ o -| ( n , -- ; f|;<V W >• 

Moving .he first term of the series to the left-hand side, di.tding by Al, and letting At -0 leads finally 
to the pdf equation 


where 


and 


£ i n'-vrSii"-. 

dt n [ =0 n L =0 j=l "/* ( 'SJ 


p (fya, t + At) - Vj(a, t))” J |Y(a, f) = y 




(72) 


(73) 


(74) 


are called the ^urselveT to^Mdorn^ processLs^iMTO) 1 such tha^all derivate moments 

I"n ^ a"; e only two possibilities for the order of 

Soong (1973)): it is two or infinity. Only the first case corresponding to all A mu ....„ L Li- > 

being zero, is of interest. The pdf equation for this case follows from (72) 


^ ± { m OW. »> 0} = ^ {»*<* a ’ 


(75) 


Cl cyj -w 

type. The increments of !T k satisfy (see Pope, 1985, Keizer, 1987) 

<AVP k (a, t)> = 0 

and <AMi(a, t)AlT,(a, {)> = 2D k ,(*)t. 

Y). It follows (see, for instance. Chapter 7.3 in Soong (197.)) t a 

H k ( y; a, t) = ^4 k (y; a, t) (76) 


and 


H w (y; a. 0 = B km( y; a, OAnnW^y; a, t) 


(77) 
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hold The pdf equation (75) is now solvable, in principle, if initial and boundary conditions are set up 
properly. It is well known (see Soong, 1973; Keizer, 1987), that the pdf equation (75) represents a 

Markovian process if W k is a Wiener process. r ... 

The conditions on the W k are now further relaxed by only assuming that the increments of \i k 
bounded. We consider jump processes independent of the continuous random processes acting on t e 
v The pdf equation for this case can be deduced from (71) (see Pope, 1985). All that needs to be 
done is to set up the transition pdf F? for jump processes on the right-hand side of (70) If we denote 
by Af/r the probability for a jump during the time interval At and denote by T(y,t t y, t + ) e 

pdf for a change of Y from y' to y if a jump occurs, then it follows that the transition pdf is given 

by 

(78) 


Ff(y'; a, r + At) = ( 1 - yWy - y') + ^ T(y'; a, t|y; a, r + At). 


Then (71) leads to 


3Fi 

dt 


dy' F,( y'; a, t)T(y'; a, f|y; a, t) - F,(y; a, t) 


(79) 


Note that r can be regarded as the time scale of the jumps. Since continuous and discontinuous 
changes of Y are mutually exclusive, we can add the corresponding contributions to the change ot the 
pdf and obtain the pdf equation for the system of stochastic differential equations 


dYj = Aj + B jk dW k + dJj, j = 1, •••,(/ + 7), 
where dJ } is the increment due to the jump process, as follows: 

SF, 


t + k {U,r,) f,,yWlyl “ 


(80) 


(81) 


This transport equation is solvable if initial and boundary conditions, as well as the properties of the 
continuous and discontinuous random processes, are specified. 

33, Closure Models for the Pressure-Gradient Flux 


The pressure gradient induces the flux 

f dp 


R = 




/ = 


dp' 

cbc. 


P = d, v = u, % = <Pj )/i + 


(82) 


in velocity space. The flux due to the fluctuating pressure gradient cannot be expressed in terms of the 
single-point pdf and therefore requires a closure model. It was shown in Section 3.1 that the pressure 
gradient for incompressible flows depends on the velocity at all locations in physical space, and, 
furthermore, it can be shown (see Hanjalic and Launder, 1972) that the correlation involving dp /i K 
consists of three parts given by 


dp' 

dx„ 


P_ 

471 


dy 


dv' R dv 


dx a \\x-y\ Jdy Y dy 


*-£■- [dy- 
8 J 


1 


dx a \\x-y\J 3y y dy„ 




vi\. 


H(x) 


(83) 


where H(x) denotes the harmonic function required to satisfy the boundary conditions. The first 
contribution in (83) is called the “return to isotropy,” the second contribution is called the as 
response” and the last boundary term (the terminology was introduced in the context of second-order 
closure schemes (Hanjalic and Launder, 1972; Lumley, 1978). The most advanced closure for (82 was 
developed by Haworth and Pope (1986, 1987) for incompressible flows with the total conditional flux 
including (82) and the viscous dissipation. No closures for the boundary term in (83), w ic is 
required for wall bounded flows, have appeared so far. 

The Closure Model of Haworth and Pope (1986). This model incorporates both the pdf flux due to the 
fluctuating pressure gradient and viscous dissipation. It was developed for incompressible fluids using 
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the Langevin approach (formulation III). The system (65,-168) is now speetalized to 

_ y 

3t (84) 

9n_IJi + vAK, + /v 
Si r ex. 

The corresponding Langevin., ype stochasttc diflerenttal equations (model, ng the dynamics of X. and 
v p) were designed by Haworth and Pope as follows: 


dX t = F a dt , 


JV, - (vA< O + C -' <K « - < v >» dl + ,C »‘ ) '' 2 dW - 


(85) 


where d\V a is the increment of an isotropic Wiener process 

(dW a } = 0, (dW t dWp) = dt S x p, 

and c denotes the expectation of the rate of disstpation of the kinetic energy of turbulence. It follows 
from (84) and (85) that 


1 8P’ 

R8Xa 


+ vAF; + F' S G 0 (V t - <F,» dt + (C 0 £) l/2 dW a 


( 86 ) 


represents the closure model. Two fundamental assumptions lead Haworth and Pope (1986, to this 
model: 

(,) The effect of the hue, nations of the surrounding lluid depends linearil, on velocity and is 
(2) The'!L^uc™rm is consistent with Kolmogorov’s scaling law for the inertial subrange. 

These two assumptions cannot be just.Tted on a 

with the quadratic dependence of the llucluatmg pre no, a Markov, an process, but the 

ponents as is apparent from (83). Turbulence is, , stnc > ^ P g ’ processes (Monin and 

fluctuations in the inertial subrange are closely approximated by P of the closure 

Yaglom, 1975). Considering the second assumption we note M * due t0 the 

model (Wiener process) represents t e s ir ™ g exte rnal forces. This agitation is not restricted to 
fluctuating pressure gradient, viscous stress , account subrange scaling 

the inertial subrange of the spectrum, whereas the model (86 ^ kcs l " ° a bounded vanables 

only. Finally, it is important to notice that implemented in linear 

such as the thermochemical scalars becaus . f realizable and relaxes to a Gaussian 

*- - ih ~ 

moments, i.e., 


— a i + a 2 T 


kt> + H 


l x0yd' 


d<Vy> 

dx s 


,t . do (t — tM and H s is a linear function of the anisotropy tensor 

where z is the turbulent time scale (t = fc/e) ana is a mica 






— 3^«/l 


(87) 


( 88 ) 


containing nine constants. Applying all exact symmetry and reduction properttes, Haworth and Pope 

S=L5 Mr-CSi ” — - - 
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model was applied to self-similar free shear layers and Haworth and Pope (1987) found that good 
agreement with experiments could be achieved for homogeneous and nonhomogeneous flows, if the 
constraint (Speziale, 1983) of proper transformation in the limit of two-dimensional turbulence was 

relaxed. 

Closure Model for the “ Return to Isotropy " Part of the Pressure-Gradient Term. It is stralght ^™^^° 
construct a closure model for (82) that simulates the return to isotropy aspect of the (luctua g 
pressure gradient, as Pope (1985) has shown. This is accomplished by pairwise interaction of material 
points with velocities u' and u" in a sufficiently small neighborhood such that: 

(1) o'" = ^(i)' -f u") remains constant and 

(2) the difference Au = u' — u" is reoriented randomly but |Au| remains constant. 

The reorientation is carried out with uniform probability on a sphere with radius |Au| that is centered 
at «" in velocity space. This interaction model can be set up in the pdf formulation as follows 
consider a volume Au in velocity space centered at u and a time interval At and let the pdf, /,(»; x > 
be approximated by N(x) elements. Let JV(u) be the number of elements with < M*. t) < o, + An* 

for a = 1, 2, 3 and 

f ( d ; x t) as N -* oo and An -> 0 . 

N 

The change of the number of elements in Au at u during At due to the stochastic reorientation of 
velocities is then 

-1 AJV(u) = AN + (u) - AN («), 

A t 

where A, •.'*!») is the number of elements added to the -olume An at n and, analogously, AtV (n) is the 
number of elements removed from Au at o durtng At, The number of added and removed elemen s 
can be established if the probability of interaction per unit time and the probability that he inte 
action of two elements produces an element in Au at u or removes an element from this volume t, are 
set up. Let t denote the frequency of interaction (time scale of interaction) and let / 2 (u , u , x, f)(Au) 
be the probability of finding elements in Au at o (,) and elements in Au at u J ; then the assumption 
of statistical independence of finding elements in Au at u<‘> from finding elements in Au at u for 
| u <‘> _ v)^*| > | Aw | is invoked. It follows that the equation 

f 2 ( u“', u 0) ) = 

holds in analogy to the assumption of molecular chaos in the context of the Boltzmann equation (see 
Chapter 2.7 in Keizer (1987)). Finally the notion of transition pdf is introduced and 

T(u (i| , u 0 ' 1 -<■ u)Au 

denotes the probability that the interaction of an element «"'* with an element vP produces an element 
in Au at u. It follows that 


(89) 


AN><") = \ - fe I /. lv U) )fM k> )T(^\ u<‘> • 

2 z i j k 


► u (,) ) 


+ y £ /, (U°') fi (V (k> ) T( v (k \ U UI - u li) ) >( Au) 3 
j k J 

holds, and by restricting T to the symmetry T(u (J) , v {k) -» u) = T(\) {k \ u 0) -*• u), we obtain 

A/V + (u (i> ) = - E X /, (vU')Mx>' k) )T(» U) , « (4> - 

T j k 

and likewise (for the removed elements) 

A JV-(u«») = - I I fi ( w<J, )/i (u (i> ) T(v u \ u (i » - u ,fc >). 

T j*i 
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Passing to the limit An -» 0, N -> oo we get 




du' 


do" /, (o')/, (»")' T(»>" - ») - /i< M \- 


(90) 


This is the general form of the effect of pairwise interaction of elements on the single-point pdf It it 

very similar to the elfec. of a jump process in he interaction of 

however, a subtle conceptual difference between (90) and m Since ^ k dilrerent ial 

two elements, those elements cannot represent lwo ,! ..o hor Honcc 190) is outside the scope of 

conation since different realizations can now be green 

in^e^i^on^e tmnsition pd^T The two geometrical properties characterizing this model lead to 

TV, »'-»)- ^ - i( ”' + ”' )l - |l1 ' “ 

!t is easy to see that T it i the e pdf ' »'* <«**'* e 9“ iva ' 

[emTo^s'crmt^m.;™ S isotro/y- L„ibn,io„ to the pressure-gradien, correlation on 
the level of second-order closures (see Rotta, 1951, Pope, 1985). 

3.4. Closure Models for the Molecular Transport Terms 

The molecular transport terms (see (63)) were shown to ^afTect ^ ^ ^negligible except in 
velocity-scalar space. Stnce the part account, ng for the effect n for the effect of 

the close vicinity of fixed wall b ^ nda " es ’ w ^ scalar spa ce. The closure for the molecular 

molecular transport on the single-point pd J , Y 

transport terms is usually called a mixing model. 

r ■ MM Dooazo (1975) put forward the closure assumption that the pdf flux in velocity-scalar 
Linear Models. Dopazo (ivo, pui un « rP i a tion corresponds to quasi- 

theoretical %£££.“* caTi/shown^ * the shore time limit ofthc genera, noniincar closure 
model discussed below (Kosaly, 1986; Kosaly and Givi, 1987). 

Nonlinear Model). Mode,, based on the interaction in t veloci iS.lTS.IS 

volumes (elements) depend in a non inear as probability of finding those elements in a 

that the interaction of two or more elements requires the ° 8 ^ prevails 

given neighborhood of the physical space. . product of the pdf/, at the chosen velocity 

and the probability of finding elements is give y P eauation (90) for the case of 

and scalar values. The interaction model has already been es ^ b | ,s ^ ed * c ise. It follows from 
pairwise interaction. It remains to determine the transition pdf space (see (63) ), 

the properties of the molecular transport terms ac ing o values and reduce variances and 

that any closure model must preserve normalization a Dooazo 1979' 

covariances. Several authors have suggested expressions for the , runs, .ton pdf T (see Dopazo, , 
Janicka et al., 1979; Pope, 1982), which can be given in the common form 

(92) 


T( y\ y" -*• y) = 


da /4(a)^[y - (i - «)y' - i*(y' + y")]’ 


where . is a random variable governed b, the pdf -4,«) This variable « 

- - 

developed. 
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CWs (1963) Model. Thts model was originally derived for .be pairwise infraction of droplets. It is 

gi, '" by d<«) » .SC - D (93) 

and r(y'.y'-*y)-d[>-t(>' + J')]- (94) 

The definition of d(«) corresponds to complete mixing in 

properttes and f interact, they emerge with J( + , > f ” '^‘^Zlncing a continuous pdf if 
Z\ ^collection of' Dnac pseudofunctions, as would be the case lor initial,, 
totally segregated flows (see Pope 1982). 

tmn t m Vrt n) 10791 The basic idea of this model is to randomize 

* - 7 - sris n 

produces a continuous pdf for initially totally segregated flows: 

/1(a) =1. ' 

but all standardized moments 


Mm = 


<)' m > 

<>’ 2 >" 


m > 2, 


, f ->4 in the limit of decaying homogeneous turbulence (Pope, 1982). Hence, freely 

S:«™b„lence does no, approach a Gauss, an pdf w„h this mistng model. 

Pope's ,19*2) Model. Pope noted , ha, the Curl and t 

their mixing history. He suggested biasing t c samp including an additional probabilistic 

^ <* » “ d lh ' ' im ‘ s 

variable age in the pdl. Hence, j x \y, ) 

between mixing interactions and /,( y; x, t) can be recovered as 

/,( y; x, t) = Jds/,*(y,s;x, f). 

The pdf f s (s) of the age variable is obviously 

/ S (s; x, t) = j dy /,*( y, s; x, t). 

normalized 


1 


ds z(s)/ s (s) = 1 


and the pdf of finding two elements with y' and y" is now biased according to 

/**( y'; X, t)/**(y"; x, t) dy' dy ", 


where 


ds z(s)/*(y, s; X, f) 


/,**( y; x, () = 

is the pdf of the elements sampled for mixing. The mixing model appears now as 

d A = if f dy' j dy"/r*(y')/**(y'') T (y'. y" - y) - /**&>} • 


(96) 


( 97 ) 
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The choice of the sampling bias z(s) enabled Pope (1982) to produce standardized moments that 
remain bounded in the limit of decaying homogeneous turbulence. 

3.5. Time Scales 

The closure models for the pressure correlations and the molecular-transport terms require a turbu- 
lent time scale t. It is clear that the single-point pdf /, of density, velocity, and thermochem.cal scalars 
does not contain information on time or length scales. Two methods have been dev.sed to prov.de the 
necessary information on time (or length) scales. 

The Mean Rate of Dissipation. The first method provides the scale information by including a 
separate transport equation for the rate of dissipation. 

dvL ( dv' 


£ = ( V 


dx* \8x 


+ 


g))- 


which in turn determines a turbulent time scale 


t = C 


(98) 


(99) 


where k = K' ; X> denotes the kinetic energy of turbulence and C is a constant of order one. The 
continuous distribution of scales in the spectrum is therefore described by a single and global scale t_ 
The exact transport equation for the dissipation rate 8 is dominated by a sens,tlve ^^J e ^ 
production, destruction, and turbulent-transport terms, which are all nonclosed. The closure o 
equation has been carried out in the context of second-order models (see Hanjahc and Launder, 197. 
Lumley, 1978; Wilcox, 1989). It represents the least-justified part of closure models on the moment 

and pdf level. 

Pdf Closures Including Scale Information. Scale information can be included in the pdf in two different 
ways Meyers and O’Brien (1981) and Sirignano (1987) suggested including a scale-determining 
variable such as the gradient of scalars in the set of probabilistic variables of the single-point pdf. 
follows that the single-point pdf then provides the scalar time scale 




r<(Vf ) 2 > 


( 100 ) 


since both numerator and denominator are integrals of the pdf of <t> and V*. Carrying this idea over 
to velocity, it becomes clear that the inclusion of all strain-rate components amounts to adding six 
probabilistic variables, which is rather expensive in any numerical method of solution. Adding the 
dissipation rate as a probabilistic variable leads to closure problems similar to the moment equations. 
This approach has been suggested by Pope (1985), who gave a Langevm formulation that ensured 
relaxation toward a log-normal pdf for the nonnegative variable dissipation rate m the absence o 
inhomogeneity. The second possibility of dealing with the scale problem is to consider multipoint 
pdfs Ievlev (1973) and O’Brien (1980) discussed multipoint pdf transport equations and closure 
methods based on mathematical realizability conditions. Spatial multipoint pdf equations have received 
little attention so far due to the enormous computational cost for their solution. However, two-time 
pdf equations can be solved numerically and yield turbulent time scales. This is in fact the second way 
of including scale information in the pdf as pointed out by Pope (1985) and Kollmann and Wu 1987). 
The two-time pdf F 2 (v\ t'; u, x, t) satisfies a transport equation identical to (56) as shown by Kollmann 
and Wu (1987). The solution of this equation can be carried out using stochastic simulation tech- 
niques Time histories of velocity are recorded and their statistical evaluation leads to integral time 

scales t because 


Rat l(* 


leads to 


., t, t') = j'du' ’I dv o a v' / ,F 2 (v', t’; v, x, f) 

-r 


r(x, t) 


R Jx, t + C) 
Raa(*,t'0 


( 101 ) 
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Note that the dissipation rate e can be recovered from 

e(x, t) = C • 


(103) 


W-“' » ilh 

measured velocity moments up to order three. 

4. Applications and Extensions 
Pdf methods have found the widest application m turbulent 

OOTl-lSw (1988). Hence on. y the case of non- 
premixed turbulent combustion will be considered here. 

4.1. Turbulent Non-premixed Combustion at Low Mach Numbers 
Turbulent non-premixed combustion is created by the interaction of one or several fuel I streams with 

Missis 

therefore restricted by 

.. r^, f* follows from 1471 bv integration over the values d of the density, 

and the transport equation for / t follows from ( ) y 8 , ■ , . and se tting s £, we 

Denoting by p [ 9l , .... «) the local relation with the thermochemical scalars and setting , 

obtain the Eulerian pdf equation 


p(<Pu---,<Pi) 


df? 

dr 

8 

dv 


+ 


o M\+ i ^-(p(<Pu--'<pM<P la 
*8x,\ M 8( Pj 

*)t) -&*(&{*%)')• 


dp 

dx. 


dx/j 


where 


j 

/ = <5(v(x, t) - u) n <%<*’ f) “ ^ 
j-1 


(104) 


(105) 


chemical models. 

Equilibrium MoMs. Chemical equilibrium prevails « 

turbulent reacting flows with fast reactions, whic is e c Damkohler number). It follows 

- - «* - 
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energy for a given mixture fraction, pressure, and enthalpy (sensible plus formation enthalpies). 
Turbulent combustion at low Mach numbers allows two further approximations: 


(A) The pressure is constant in thermodynamic relations. 

(B) Enthalpy is a linear function of mixture fraction. 

The first approximation is based on the observation that, for low Mach numbers and small mean 
pressure gradients, pressure fluctuations are small 


O 




« 1 . 


Hence we can neglect the pressure variations due to fluid mechanics in the thermochemical aspects of 
the reacting system. This does not hold for supersonic and rapidly expanding and compressing flows, 
where the pressure fluctuations can be of the same order of magnitude as the mean pressure itself. The 
second approximation can be seen to be a consequence of the energy equation. If the Mach number is 
much less than unity, the dissipation function in the energy equation (13) can be neglected; if the 
pressure is approximately constant, the substantial derivative of the pressure in the energy equation 
can also be neglected. Hence 

Dh ^ dq a 
^ Dt ~ dx a 

is obtained. This form of the energy equation reduces to the same convection-diffusion equation for 
enthalpy as for the mixture fraction under the standard restrictions (Bilger, 1976) of the equilibrium 
flame sheet model. It can be argued (Bilger, 1976) that a linear relation between h and mixture 
fraction { exists sufficiently far away from boundaries, where h and £ do not necessarily satisfy such 
a linear relation. Again, this approximation does not hold for supersonic flows where <J> becomes 
significant and the enthalpy ceases to be a conserved quantity. It follows now that the pdf equation, 
for turbulent flows involving fast reactions and low Mach numbers, contains only a single conserved 
scalar (mixture fraction) to account for combustion. The value of the mixture fraction uniquely 
determines the density, temperature, and composition so that expectations of these variables can be 
obtained by a straightforward integration. Pdf methods have been applied successfully to this case 
and the closely related turbulent transport of passive scalars (Kollmann and Janicka, 1982), but the 
full power of pdf methods has not been exploited in the transport of conserved scalars. Furthermore, 
the condition of high Damkohler numbers is too restrictive for many fuels (in particular hydrocarbons 
(see Bilger, 1988)). 


Nonequilibrium Models. The inclusion of chemical nonequilibrium for reaction systems proceeding 
with finite rates requires careful analysis of the usually complex system of reactions in order to obtain 
a tractable system. The tools available for the reduction of complex systems of reactions are partial 
equilibrium for selected steps and constrained equilibrium, where the progress of reactions is assumed 
to take place in a series of quasi-equilibrium states subject to a set of constraints which are controlled 
by the rate-limiting reactions (Keck, 1978). These tools are not perfect and considerable insight into 
the detailed reaction mechanism and calculations of laminar flames are required to establish a 
reasonably accurate simplified mechanism (Kee and Peters, 1987; Peters and Williams, 1987; Rogg 
and Williams, 1988). Pdf methods were applied successfully to hydrogen-air flames (Chen and 
Kollmann, 1989a), CO-air flames (Pope and Correa, 1986), methane-air flames (Masri and Pope, 
1989; Chen et al , 1989), and propane-air flames (Jones and Kollmann, 1987; Chen and Kollmann, 
1989b). Other reactants were treated by Givi et al (1985) and Arroyo et al (1988). We consider the 
case of turbulent methane combustion in jet flames in some detail to show the properties of pdf 
methods following Chen et al (1989). The combustion of methane with air is described with the 
four-step mechanism of Peters and Kee (1987), which requires five scalar variables ^(x, f). They are 
defined as follows: 

T t = £ (mixture fraction), *¥ 2 = n CH4 , 

^3 = rt CO> ^4 = H, 

^5 = *H 
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(n denotes the number of moles per unit volume). The pdf equation (104) is integrated over velocity 
space 

(. dt dxp pi d(pj 


<P>- 


= «P><®, 


GX„ 


( * d / d ( dij/j 




(106) 


and the density-weighted pdf /, defined by 


/ 1 = 


PiVu <Pi) 
<P> 


•••. <Ph x > 0 


(107) 


is used for essentially formal reasons. The turbulent fluxes and the time scale are determined from a 
second-order closure model (Dibble et at., 1986). The closure for the flux is given by 


= %>fi = c '<p>- g v &}£- 


(108) 


The mixing model (scalar-dissipation model) is the nonlinear interaction model of Dopazo (1979) and 
Janicka and Kollmann (1979) described in 3.3, equation (91). The rates of the four-step mechanism 
determine the source terms , ■■■, <Pt), which control the motion of material points in scalar space 
due to combustion. The boundaries of the set of realizable states in scalar space are rather intricate 
(see Chen et at., 1989). They satisfy all the conditions laid out in Section 3.2 for the pdf formulation. 
A detailed comparison of first- and second-order moments with the experiments of Masri et at. (1988) 
was carried out by Chen et at. (1989). A sample of two-dimensional pdfs and characteristic functions is 
presented in Figures 1-18. The flame considered is a turbulent jet flame burning methane with air 




Figure 1. Pdf of mixture fraction and temperature in a turbu- 
lent methane-air non-premixed flame at x/D = 20 and r/D = 
1 . 11 . 



Figure 2. Pdf of CO mass fraction and temperature in a 
turbulent methane-air non-premixed flame at x/D = 20 and 
r/D = 1.11. 
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Figure 13. Pdf of mixture fraction and temperature in a Figure 14. Pdf of CO mass fraction and temperature in a 

turbulent methane-air non-premixed flame at x/D = 20 and turbulent methane-air non-premixed flame at x/D = 20 and 

r/D = 2.07. r/D = 2.07. 
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from a coflowing stream. The inner diameter of the jet pipe is D, - 0.0072 m and the bulk exit 
velocity of the fuel methane is u 0 = 48 m/s. The reaction is stabilized with a pilot flame (C 2 H 2 and 
H 2 ) in the annulus between fuel jet and coflowing air stream (D„ = 0.018 m and u ail = 15 m/s). The 
solution of the pdf equation was carried out using a stochastic simulation technique (see Chen et al., 
1989). The cross-section at x/D, = 20 is selected and results at three radial stations are presented in 
figures 118. The pdf for mixture fraction and temperature in Figure 1 (r/D, = 1.11), Figure 7 
(r/D-, = 1.49), and Figure 13 (r/D, = 2.07) show the structural change of the pdf as the radial location is 
shifted outward. A significant amount of nonequilibrium is observed at r/D, = 1.49 in Figure 7, where 
the pdf spreads over a range of 500 K on the rich side of the flame. Much less spread is seen on the 
lean side in Figure 13. The pdf of the CO mass fraction and temperature in Figures 2, 8. and 14 
exhibits the change from a narrow pdf at r/D, = 1.11 to a broad and doubly ridged form at 
r/D, = 1.49 and finally to a club-like shape at r/D, = 2.07. The pdf of CO and methane mass fractions 
in Figures 3, 9, and 15 reflects the amount of coexistence of these components. The pdf is spread over 
a wide range of methane mass fractions on the fuel rich side at r/D, = 1.11, but gets gradually 
squeezed toward zero as the radial location is moved into the lean side of the flame zone. Figures 
4 6, 10 12, and 16-18 contain the characteristic functions in separate graphs for real and imaginary 
parts corresponding to these pdfs. It is noticable that the characteristic functions are much smoother 
than the pdfs (no smoothing of pdfs or characteristic functions was performed). Furthermore, it should 
be remembered that discontinuities of the pdf appear as oscillations with decaying amplitude in the 
characteristic function. The characteristic function is real valued if the pdf is symmetric. Hence, the 
imaginary part of the characteristic function is a measure for the skewness of the pdf. Finally, note 
that the primed variables in the figures are the independent variables in Fourier space corresponding 
to the unprimed variables in scalar space. The characteristic function for mixture fraction and 
temperature in Figures 4, 10, and 14 shows a ridge, which decays with distance from the origin with 
increasing speed as the radial location is moved to the lean side of the raction zone. The characteristic 
function for CO and temperature in Figures 5, 11, and 17 is centered around the origin (real part) 
with quickly decaying waves, except at r/D-,= 1.49 indicating stronger discontinuity in the pdf. The 
characteristic function for CO and methane in Figures 6, 12, and 18 shows slowly decaying waves in 
real and imaginary parts as r/D, = 1.11 and r/D, = 2.07, but much less waviness at r/D, = 1.49 which 
indicates a smoother pdf at this location. 

This example shows clearly the degree of chemical nonequilibrium that can be present in a 
hydrocarbon flame. The fact, that all thermochemical variables are bounded, together with the 
unusual forms for the pdfs obtained in this example, demonstrates that assumptions like the quasi- 
Gaussian one are untenable for scalar pdfs. 


4.2. Turbulent Supersonic Flows 

Compressible turbulent flows pose new and challenging problems, in particular in the area of 
supersonic combustion. Pdf methods can be applied fruitfully to this type of turbulent flows, but the 
research effort is still in the beginning stage (Farshchi, 1989). The basic theory for the pdf approach is 
outlined here and possible ways of its application to supersonic combustion problems are discussed. 
The case of compressible nonreacting flows of an ideal gas is considered first. 


Turbulent Supersonic Flows Without Reactions. The balances for mass, momentum, and energy set up 
in Section 2.1 are transformed to dimensionless variables. Mass balance appears in unchanged 
dimensionless variables ((6); only the Eulerian frame is used in the following), but momentum and 
energy balances contain several dimensionless parameters. The momentum balance is now given by 


where 


Dv* = 

} Dt “ 
Re = 


- y- + 1 + W- 

dx a Re exp 

(109) 

Po u oL p _ foL 

, D = Y ' 

Po u 2 

(HO) 


with p 0 > L, /i 0 , T 0 , k 0 , f 0 denoting the reference values. The stress tensor appears in unchanged 
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dimensionless variables (12). The energy balance emerges as an equation for the internal energy e in 
the form 

--(•/- DO + yMlp ) + ytt - 1 " 

where the ratio of specific heats y was assumed constant and 

M 0 = — , a 2 0 = yRT 0 , Pe = Re • Pr, Pr = ^ 

a 0 K o 

are the Mach, Peclet, and Prandtl numbers for the reference state. The energy balance can be recast in 
terms of the specific (dimensionless) entropy 


Ds 

P lh ~ 


P cdxATr b U R eT Per 2 dx. 


( 112 ) 


The dissipation function <D in dimensionless form is given by (16) and the heat flux is 

cT 




dx' 


where k(T) is the dimensionless conductivity. The ideal gas equation appears in dimensionless form 


as 


yMoP = — pe — 1. U 13) 

Cv 

Single point pdfs can now be defined in several ways for compressible flows (since mechanical and 
thermodynamic state are completely specified) if the velocity and two thermodynamic state variables 

are known. In the present case the following set of variables is considered (for a different set, see 

Farshchi (1989)): Velocity v, density p, internal energy e, and rate of volume expansion D = dvjdx, x . 
The tranport equation for D follows immediately from the momentum balance 

DD \ d /I 8v x dv, + g _ JLfl ) (114) 

Dt Re dX a \p exp ) dxp dx a dx x \p dx^J 

The definition of/ (see (48)) is now modified as follows: 

/ = d(v — v)S(p — d)S(e — u)S(D — £), (115) 

and the pdf equation for compressible turbulent flow without reaction emerges as 



This form of the pdf equation contains information on a special time scale, namely the time scale of 
volume expansion. The inclusion of the inverse D of this time scale is the reason for the evolution- 
type structure of the pdf equation (116). The rate of volume expansion D could be eliminated, but 
then a pdf equation of the hyperbolic type would result. 

Turbulent Supersonic Flows with Chemical Reactions . Turbulent flows at supersonic speed, with 
combustion reactions, lead to several new problems. The mixing of fuel and oxidizer becomes a 
problem of central importance because of the reduced residence times at high speeds and the effect of 
compressibility to reduce the turbulence level (see Papamoschou and Roshko, 1988; Papamoschou, 
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,989, and .hereto* .he intensity of mixing- Turbulent Dow at "SSSI 

significantly by the in, erection with s oc • — turbulent shear layers. Pdf 

shocks (called shocklets (Johnson et al., 19 )) ..... . com bustion. We consider first 

methods can be adapted to include the effects o p . 4 , (Equilibrium Models) that three 

the case of infinitely fast react, ons. It was shown m Sect on 4d jbqu.l,^ pressure can 

variables determine the local state in t is case. mix u ^ conse ’ rved due to frictional heating in high 
vary significantly in supersonic flows an py number subsonic flames, 

- — - — miMre 
fracUon {; * { , x. „ - <3(. - - *m - *w> - c>«« - -i». " ,7 > 

.deal gas re,a.,o„s aUow replaceuren, of presto a^d en,Mpy ^ 

transport £ rST —immic properties follow fro. ,e pdf /, by 
integration. The mean temperature for instance is given by 


<T> = ^dn j" tiff ^ drj T(n, <r, r])J\ (ft, a, tj). 


where T ( , <r, ,) denotes .he local relation o, .e.pera.nre £2Z 

subsection ‘However^no dosuremoS have' ten developed so far for (116,. o, i.s coun.erpar. for 
reacting flows, and much work remains to be done. 

s „«cte f~. The interaction of 

Of practical importance (Bilhg and ^ugge^ 1 ^ So fi* ^ and mean shock properties 

1985) have been used to predict mean fields in e 8 aDDlication of pdf methods to 

and location were calculated with shock-capluring q ^ ncar dls continuities for finite 

y-rs.tr sr. — - ■ — 

<D(x, t) is considered. The equation 


£0) BF 

dt + 5x 


S = 0, <D(x, 0) = <D 0 (x), 


(118) 


where m*,t) is the ^ U )^ a j da " ll 'f 984 | f ' if'the'innia'l'co'ndnion 'ir chosen randomly from a set of 
^r^I^tici properties of the solution can be described in terms of the 

pdr/ ' ,v;x ' ,): /.6(®<x, .)-«>>, /, = </>• 

„ the „ux F and the source , are local functions of ®. and if the solution • remains a, leas, once 
continuously differentiable, then the pdf equation follows from (47), i.e., 

(119) 


+ + |~(S(<f>)/,) = 

dt d<& dx dtp 


g 2 F g<J> » 

dO 2 dx' 


However, the right-hand side is nonzero if the flux F J^^dm^st^Ted foMhirci^The theorj 
sutted for the analysis of discontinuities and a d. a fixed location x. The 
of stochastic differential equations allows analysis of the scalar <p 

temporal increment of <J> is then ^ = ^ + + dW (120) 

“w. 
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equation (118) leads to 


d<t> = 


dt + Sdt. 

dx 


The source term S dt can now be viewed as the sum 

Sdt = S 0 dt + dW (121) 

of the differentiable contribution S 0 dt and a continuous but nondifferentiable part dW. The flux term 
requires a closer look. If F(<D) depends on 1) in a nonlinear fashion such that discontinuities form in 
finite time from smooth initial conditions, then there exist random and discrete-time instances when 
discontinuities cross the fixed location .x. Hence, the derivative of the llux is the sum of a singular and 
a continuous part 

= I MMt - '.) + 

ex 



where [<D] f denotes the jump height at time and ( dF/ex) c denotes the continuous part. Hence, 

dP = X Mi dt S(t - r ; ) (122) 

<t 


is the increment 
identified as 


of a jump process and the increments in the stochastic differential equation are now 


along with (121) and (122). The pdf equation for the solution process of (120) can be deduced from 
(81) as follows: 


3/i 3 

“r h -r — 

dt cep 




1 c? 

2 dep 


2 {Bf ,] 




+ <1 I d<t>' j\m 7(<D' <D) - /, 


(123) 


if dW is specialized to a Wiener process. It becomes clear by inspection of this equation that the 
discontinuities crossing a given location x affect the evolution equation for the pdf in integral form 
appropriate for jump processes. This integral requires the probability of a jump from 4>' to O, denoted 
by T(d>' -♦ <t>), and the time scale t for the appearance of discontinuities at x. The pdf equation derived 
from (120) does not provide this information, because only a single location x is considered and 
x-derivatives constitute therefore new unknowns. 

This example showed that the appearance of shock waves with random location and strength 
produces an integral contribution to the pdf equation. It can be expected that the time scale and the 
transition pdf 7(0' -* <t>) are functionals of the flow variables. The explicit form of this functional 
relation is unknown at present. 


5. Methods of Numerical Solution 

The pdf /,(d, u, tp t , . . ., <p m ; x, t) is apparently a function of a large number of independent variables. 
Classical methods of numerical solution such as finite-difference or finite-element algorithms become 
prohibitively expensive because the numerical effort for the solution grows rapidly with the dimension 
of the domain of definition (number of independent variables). However, stochastic simulation tech- 
niques can be shown to grow in numerical effort only linearily with the dimension of the domain of 
definition. Hence they offer the possibility of a numerical solution of the pdf equation for a significant 
number of variables (up to about ten with current computational capabilities). The development of 
stochastic simulation techniques for the solution of the pdf equation is essentially due to Pope (1985). 
The basic idea is to represent the pdf by a sufficiently large number of notional particles, whose 
motion in physical and velocity-scalar space is governed by modeled transport equations. The closure 
assumptions introduced for these dynamical equations (see Section 4) correspond to the closure 
assumptions constructed for the pdf equation. The numerical algorithm is based on a fractional-step 
method, where in each fractional step a numerical operator corresponding to a distinct physical or 
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chemical process is applied. The motion o( ' hl ; . rtta^the Sond panicles in 

The pdf is then gi.en approximately as a histogram of he ^opernes ; ol £ ^ ^ 

sufficiently small neighborhoods m physical space A detailed ten „987) and 

in Pope (1985). Several variants of this baste , onS! of a Lte-difference algorithm for 
Chen and Kollmann (1989b) developed a hybri m . „ jmu | at ( on technique for the scalar 

a second-order closure , fcr the veloc.U technique to Hows in 

-ate method for ,to-„pe 

stochastic differential equations (Haworth and Pope, 1986). 

6. Conclusions 

Single-point pdf methods can be deduced 1 from i exao ^and closed 

pdf a, a single point or a. a I Bn* taken a. a special argumen, 

characteristic function which is in fact essentially equivalent ways to formulate 

(unction. This has ... important implications: them are two equally equi ^ ^ 

the turbulence problem at a single point or aim statistjca | treatI nent of turbulence based on a 

S^^SZ=t= rr^nden, va„ab,es. This ,s fact 

the main advantage of pdf methods over moment methods. Lagrangean frames and the 

Pdf or characteristic function methods can be set up “ “^nsidered in detail. The 

choice of the frame is a matter of convenience ® because molecular transport 

equation for any pdf at a finite num er P terms reauire information given at least 

(viscous and diffusive phenomena) and t e pressure-gra bJ j whic b was formulated in 

one additional point for Iheir description. Hen e, » „ P ”s ,'te,io"s a third formulation using 

terms of pdfs and characterisiic functions. Unde, * discussed and their 
stochastic ditferential equations can e given- satisfy, at best, the mathematical 

=r r.«s of’, he exact terms the pdf 

^^df^methods have ^^^^a^^Th^combus^ion ot ^than^with^iT^a^umbuf^ 

combustion are currently an active resear * f results for pc jfs and characteristic functions 

- the 

bounded scalar variables used in such a case. ™ m „ rP< v S ; h le turbulent flows and the interaction 

Further developments of pdf methods may of compressible lurbulen, 

sSetami ZEEt ^“scV^'and the .Tec, of discontinuities on the pdf was shown 
to appear as an integral contribution in the pdf equation. 
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Abstract 

The interaction of turbulence and chemical kinetics is examined here with the emphasis on 
the influence of turbulence on chemical reactions. Both nonpreinixed and premixed flames are 
considered. In particular, a nonpremixed methane turbulent jet flame is used to elucidate the 
complex nature of the interactions between turbulence and chemical kinetics. Furthermore, 
this example provides a useful evaluation of the mixing properties predicted by a probability 
density function (pdf) method. 
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1.0 Introduction 

tw' fal r f a, ti ° nS " m,rs in ,urbulem flow., over 

chemical regions clnTead o ifferen '°T **"" 1,."^'“ l>e ' We " n ,urM »™ 
and chemical reactions mav simnlv ni |V’ TT™' "«* '““actions between turbulence 
(Williams 10S9I St. •. p od, fr lke bailie slightly causing wrinkles of flame surface 
chenuv^eae, * T™' 10 ' 18 a significant modification in both the 

thefow 2 n r, h T ,urbuleu 7- If rf"-'"ical reactions cause small density changes in 
mavs^bat Urb u UlfnCe 15 Weakl >' b >' *he chemical process but the turbulence 

SSST ° n * m r' rraf ' i0 “ S ' H ™ «be Impose of combustion 
which can alter the ft ' j T one ex P ects lar 8 e density changes (i.e., an order of magnitude) 
the entr flmd dynamics significantly. It has been observed experimentally tha 

^~jr:z;z™ n z t - ba ™ »>»— *** 

et at. 1985; Dimotalds 19891 O !l 7* expected m constant density flows (Hermanson 

a point that chemi l • n the other hand, strong turbulence can strain the flames to 

flame to extinguish 8 r<?aCt '° ns can no lon *cr keep up with the mixing process causing the 

flame extinction cant^Trior ^7“* n ** y 09881 haVe —W ' ha < 
between turbulence and chemistry. ‘ ‘ ,,,d,Cat,n S a 5,I ™S interaction 

it is IseZ^o !rl a f" <l ,r iUlt | ify ‘ he , COm P ,ex '““actions between turbulence and chemistry. 
ovlSlamcteta, 5 f.l ^ “ d Hme ^ “ turbulent reacting flows. An 

lh t“T bet "' een turbulence and chemical reactions can be 
scale) versus the R . - i t" 1 ' ° , numb “ tke ratio of flow time scale and reaction time 
Based on the length 00 7 num >er over tbe whole range of length scales (Williams. 19S9). 
O^elx^emt crith he a' ° *7“ and '“bulence. two extreme regimes are identified. 
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of the Tt , t b kngth 18 'ticutihed as the distributed reaction regime. The nature 

to be eX™ L^f I 7 ,m T, betWeen theSe t "'° “ r “ber complex, and i,““ 

operation condition Tu 6 y ’ P ractlcal combustion systems involve a wide range of 

operation conditions including the intermediate regimes. 

inherent in^varioi ? ^ ^ P8Per *. S t0 provide a fundamental understanding of the physics 
^^o^:r CeSSeS Causi P® turbulence to interact with chemical kinetics. In View 
present naoer i d practlcal chenuca l processes that occur in turbulent flows, the 

“ a backgro H r d ma,nly u t0 tlle °f turbulence on chemical kinetics. To 

s^rntLs he b ,h T eS in tUrb,tlrat acting flows, Section 2.0 reviews anti 

on the P 1 • t 10 P iT CS aWS for tbcuucal reactions in mixtures of ideal gases based 
discussed. enan frame ' The corrt ' s P° ndill S transformation to the Lagrangian frame is briefly 

reactions First t0 ^7 "T' 1 ! dlscussl °n s on the influence of turbulence on chemical 

tratd and he de 7% * a -acting system (two species) is illus- 

pdf transport eciu^ 66 ° ^ 1S characterized by the segregation parameters. Then the 

Lution demLT^ . r *, Smg ' e , SC,l * r “ ^loml and an analytic solution is derived. This 
solution demonstrates the close relation between chemical reactions and the scalar dissipation. 



Next, nonpremixed flames are considered and the appropriate measures of mixedness are in- 
troduced. The essential issue for modeling nonpremixed flames is highlighted by examining 
the closure problems in the current modeling methods. In particular, the mixing models in 
the pdf methods are examined in depth and future developments are indicated. The analysis 
of premixed flames is limited to the flamelet regime. Two proposed theories are reviewed: 
the Bray-Libbv-Moss flamelet theory (Bray et al.. 1985) and the coherent flame theory of 
Marble and Broadwell (1977). The regime of flame sheet combustion clearly indicates the 
fundamental importance of surfaces embedded in turbulent flow field. Surfaces relevant to 
combustion flows are introduced and classified according to the underlying transport mecha- 
nisms. The relative progression velocity is shown to be dependent on a number of processes 
including diffusion, chemical rates and the local scalar gradients. 

In Section 4.0, the intrinsic topology of surfaces embedded in a three dimensional space 
is discussed. The effect of chemical kinetics on turbulence is briefly described. Section 

5.0 is devoted to the discussions of nonpremixed turbulent methane jet flames, which have 
been studied extensively by experiments and numerical simulations. The results from a 
stochastic simulation of the joint scalar pdf equation permit us to evaluate the mixedness 
parameters introduced in Section 3.0. It is shown that the mean chemical reaction rates 
can be larger or smaller than their corresponding quasi — laminar values by several orders 
of magnitude indicating the strong influence of turbulence on chemical reactions. The last 
section summarizes the main findings from this study. 

2.0 Basic Equations for Turbulent Reactive Flows 

The current analysis of the interaction between turbulence and chemical reactions is restricted 
to Newtonian fluids in gaseous phase, to which the thermodynamic relations of ideal gases 
is applicable. Given compositions, two independent (intensive) thermodynamic variables, 
and velocity, one can determine the thermodynamic state of a reactive mixture. The gov- 
erning transport equations for compositions, temperature, and velocity are dictated by the 
conservation laws of mass, energy, and momentum. In the Eulerian frame, these conservation 
equations can be written based on an observer fixed at an arbitrary location x in the flow 
field as function of time t. For some aspects of turbulent combustion, the physics can be 
better described based on the Lagrangian frame (following the fluid particles) than on the 
Eulerian frame. The relations that bridge the Eulerian and the Langrangian frames will be 
given at the end of this section. 

Mass Balance: 

Conservation of mass leads to the following transport equation for the density p(x. t) on the 
Eulerian frame 




(pv a ) = 0, 


where the repeated indexes mean summation over all possible states. 


Species Balance: 


( 1 ) 


3 



We consider a mixture of N ideal gases with its compositions described in terms of mass 
fractions Y,(.r , f). When chemical reactions occur, the mass fractions. }*,. are not conserved, 
but consumed or produced according to their net production rates Q , , which are determined 
by the reaction mechanism. Note that the net production rates, Q t , depend only on the 
local thermodynamic variables; that is, Q, do not contain time derivatives of thermodynamic 
variables or their integrals with respect to time or space. The transport equation for the mass 
fraction of the /-th species is then given by 


DY, 
’ Di 


dJa 


+ pQ / = 1..JV 


where J' a denotes the diffusive flux in the ev-coordinate. For multi-component reacting sys- 
tems, the diffusive fluxes can be expressed in terms of functions cont aining gradients of species 
concentrations and their binary diffusion coefficients. If the i-th species is sufficiently diluted, 
its diffusion fluxes can be approximated by the Fick’s diffusion law. 


J'a = -pTi 


dYj 

dx a 


(3) 


where T* is the diffusivity of the i-th species. For the purpose of this paper, we will use 
this approximation as it provides a very simple formula for calculating the diffusive fluxes in 
reacting flows. 


Momentum balance 

Newton's second law leads to the balance equation for momentum. 


P 


Dv, 


Dt 


dp 

dx a 


+ 


d.r j 


+ pfa , 


( 4 ) 


where r a $ is the stress tensor and f a is the external force per unit mass. For Newtonian 
fluids, the stress tensor obeys the following constitutive relation: 


Tad = P( a — + 


dv Q dv 3 2 dv y 

dx 3 + dx Q 3 ali dx 1 


), 


where p is the dynamic viscosity. An important feature of equation (5) is the linear re- 
lationship between the stresses and the rates of strain which are expressed as the velocity 
gradients. 


Energy balance and state relations 

Application of the first law of thermodynamics to a differential control volume leads to the 
energy balance equation in the Eulerian frame. This equation can be expressed in terms 
of several equivalent forms depending on the choice of the thermodynamic variables. If the 
specific enthalpy h{x_,t) is chosen, the energy conservation equation can be written as 


Dh Dp dq a 

P ~Dt = ~Df + ~ JT a 


( 6 ) 


4 



If the specific internal energy e 
equation for e 


= h 

De 


P 


Dt 


— pj p is used, one can derive the following transport 


dl'a , , dq a 
-p— + $ - 


d.r 


d.r a 


(7 


where $ and q a denote the dissipation function and the energy flux in the o -coordinate 
respectively. The dissipation function $ is defined by 


$ = 


dv, 


7~Ct J 


d.r, 


( 8 ) 


which represents the heat generated by mechanical dissipation due to viscous friction. The 
energy flux q a is expressed as a combination of conductive, diffusive, and radiative fluxes. 


q Q = -k 


dT 

dx Q 


_ V* Alt i n R 

P / 4 i r 


i=i 


M, dx c 


( 0 ) 


For multi-component systems, with the ideal gas assumption, the mixture enthalpy is 
simply the sum of individual specific enthalpy weighted by its concentration 


h 



( 10 ) 


where h, is the molal enthalpy and M, the molecular mass of the i-th component. The molal 
enthalpy hi consists of two parts: the formation enthalpy h® and the sensible enthalpy 

T 

h t = h° t + / dr'cpcn, (ii) 

J To 

where c p (T ) denotes the molal specific heat at constant pressure. The above system of 
equations is completed if the ideal gas equation 


P = r*T±± (U) 

,=1 A1 ‘ 

(where -ft denotes the universal gas constant and M , the molecular mass) is included and the 
chemical sources Q, are specified. 

So far we have presented the basic conservation laws in terms of mass fraction 1 enthalpy 
h , density p and velocity r. This set of variables may not always be the most convenient 
ones, and a linear or nonlinear combination of these variables could be used for the treatment 
of turbulent reacting flows. For low Mach- number flows, it has been shown (Pope, 1985), by 
using a Taylor series expansion of the state relations, that chemical sources are, to the lowest 
order, independent of pressure fluctuations and that, the substantial derivative of the pressure 
in the energy equation (6) can be neglected except under the condition of strong pressure 


o 



variations. Consequently, for low Mach number flows, the set of thermo-chemical variables 
t) obey a similar type of transport equations as 


Dil\ 

> 

Dt 


d 

d.i'o, 


pTi 


dvj_ 

d.r n 


+P ( U 


i = 1 
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For combustion processes at constant pressure, the total number of thermo-chemical variables 
is simply / = N + 1. 


Lagrangian frame 

Studies of Chemical reactions in turbulent flows can also be carried out in the Lagrangian 
frame (Borghi, 1988), in which the observer follows an arbitrary material point of the fluid and 
monitors the evolution of this material point. This approach has some advantages because 
it is the natural frame for phenomena that are dominated by time history. Hence, the 
transformation rules between the Eulerian and Lagrangian frames will be given, and the 
structure of the conservation equations for the thermo-chemical variables <5, in the Lagrangian 
frame will be discussed briefly (Monin, 1962). In the Lagrangian frame, the independent 
variables are time t and a, which is a variable used to identify the material point. The 
common choice for a is the position of a material point at the initial time to- The position 
of a material point in the Lagrangian frame is denoted by X_(a, t) which can be used as a 
transformation function between the Eulerian and Lagrangian frames in the following manner 

L — nnc/ a = X~ l (x,t), 


where X -1 denotes the initial position of the material point which moves to the position x at 
time t. Here, the upper case letters denote the dependent variables in the Lagrangian frame 
and the lower case letters correspond to the same variables in the Eulerian frame. If the 
mapping function X. and its inverse function ^V -1 are both twice continuously differentiable, 
both the time and spatial derivatives in the Eulerian frame can be transformed into the 
Lagrangian frame, and vice versa. For instance, the gradients of a variable can be transformed 
from one frame to the other frame according to (Truesdell, 1954) 

d 1 (14) 

dx a 2 Q 0r) q Quj d a6 


and 

d J dXj'dX; 1 d 

da, ~ dr, dr „ dr,’ 

where J denotes the Jacobian determinant 

1 dXf, dXrf <9A' W 

^ O 3 ^ 6 71 uJ O o 

6 oa Q aaj oa ^ 


(13) 


(16) 


Note that repeated subscripts imply summation and that is the permutation tensor. 

The transformation formula for the second and higher derivatives can be derived by using 
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Eqns. (14) and (15) repeatedly. In particular, the relation for the Laplacian is .summarized 
here as 

_ _L dx < dx « d ( 1 dx < dx * d \ 

* 2J €a,f3 Q (lr) Q a ^ Qf tfi yj Qa 3 Q a ^ da n ) ' 

The time derivatives in the Lagrangian frame are of fundamental interest because for a given 
material point, velocity and acceleration are. by definition, the time rates of change of position 
and velocity respectively. The latter terms appear in the Eulerian frame as the substantial 
or the Stokes derivative 




+ 1 1 


d 


•T, 


D_ 

Dt' 


(IS) 


The transformation rules (15)-(18) allow us to derive the transport equations for thermo- 
chemical variables in the Lagrangian frame based on those in the Eulerian frame. For instance, 
the transport equation ( 13) for $ ,(r, f ) can be transformed into the Lagrangian frame by using 
Eqns. (17) and (18). and the result is 

£ii = _L dx < dx » d f &rjdXcdXtdvA 

dt 2R 0 da^ da^ da g \ i? 0 da# da 7 da a ) " 


i = l../. (10) 

Note that the nonlinear convective terms disappear, but the diffusion terms become highly 
nonlinear as X_ is a dependent variable in the Lagrangian frame. In deriving Eqn. (19). we 
have expressed the mass conservation law in the integrated form as 

i Q) _ t ( 9Q \ 

R(a,t) 

where R{a,t) = p(x_.t) for r = X_(a,t). With this relation, the Jacobian J can be eliminated 
from the transport equations. 


3.0 The Influence of Turbulence on Chemical Reactions. 

In turbulent flows, chemical reactions proceed in an environment changing randomly in time 
and space. The fluctuating nature of mechanical and thermodynamic variables can have sig- 
nificant influence on the progress of chemical reactions. In this section, we will analyze this 
influence in detail by examining the mean reaction rate for a single irreversible chemical reac- 
tion involving only two reactive species. Under such a circumstance, the transport equation 
for the probability density function (pdf) of a single thermo-chemical variable can be solved 
for a non-decaying homogeneous turbulence. With the help of this solution, one can clarify 
the role of the correlation between the scalar and its gradients in describing the interactions 
between chemical reaction and turbulence. 

The progress of chemical reactions in nonpremixed flames depends strongly on the degree 
of mixing between fuel and oxidizer. Turbulence can greatly enhance the mixing process and 
thus increase chemical reaction rates by several orders of magnitude. Due to the practical 
importance of nonpremixed flames, we will discuss in depth the mixing phenomenon and the 
theories that describe this important process. Furthermore, special attention will be given to 
the mixing models which are currently used in the pdf methods. 

In premixed turbulent flames, the influence of turbulence on chemical reactions may 
be described as the combined effects of convection and distortion of the reaction zone. We 
will restrict our attention to the special case of thin reaction zones, which can be treated 
as flame sheets. As the concept of treating turbulent flames as an ensemble of laminar 
flamelets alleviates the need for detailed modeling of chemical reactions, we will provide 
some discussions on the kinematics and dynamics of surfaces moving through the turbulent 
flow field. Several types of surfaces relevant to turbulent combustion flows will be considered 
and the topological and geometrical properties of these surfaces will be analyzed. 


3.1 Mean Chemical Reaction Rates in Turbulent Flows. 

Let us consider the transport equation (2) for chemical species i. Density- weighted average 
of this equation leads to the mean transport equation for Y t as 




( 21 ! 


The most interesting quantity that will be addressed here is the mean chemical reaction rate 
w, as it depends strongly on the temperature and species fluctuations. We begin with an 
analysis of this dependence for a bimolecular, irreversible reaction between two species A and 

B 

A + B C. 


The instantaneous kinetic source term w, is written as 

U' A = u'B = ~k Y aYb ■ 


( 22 ) 


For simplicity, let us consider the case with a constant k first. Average x>f Eqn. (22) yields 
the mean kinetic source as 

w A = -HY a Yb + 1'A'Vb). 
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It is clear tha t the influence of turbulence on the chemical reaction is reflected in the cor- 
relation g . Application of the Schwartz inequality leads to the upper bound for this 
correlation as 


ynyn 

1 A 1 B 


~ 1 A 1 B * 


With the inequality 

Y" 2 < V'(l-y') , 1 = .YB, 


which is the consequence of Hausdorff ’s theorem ( Akhiezer, 1965) on the moments of bounded 
random variables, one obtains the following relation 


yPfl' 2 < Ya(1 - Ya)Y B (1 - Yb)- 


Hence, the upper bound for the mean chemical reaction can be derived as 


u>a\ < k(Y A Y B + \J Ya{ 1 - Ya)Yb( 1 - Yb))- 


(23) 


Therefore, turbulence can increase the mean sources compared to the corresponding 
quasi-laminar values However, the opposite effect is also possible; that is, turbulent 

fluctuations can also reduce the mean reaction rates. If the correlation Y ^ I q is negative, it 
follows that 

|« a | < \w(Yi)l 

which is the upper bound for the mean chemical source term. We can obtain a lower bound 
for |um| as follows. If the reactants A and B be totally segregated, the joint pdf of 1 A and 
Y b is given by 


f(VA< ys) = (l - VtMy.-tMl - ub) + - yA)^{yB)- (24) 

This pdf represents a special situation where species A and B do not coexist at the same spatial 
location, but they may have nonzero mean values. It can be shown that the fluctuations are 
maximal in this case. Physically, the flow consists of randomly distributed regions with 
only species A or species B but not with both, and the two species are separated by an 
infinitely thin interface. Straightforward integration of the instantaneous chemical reaction 
rate weighted by the pdf over the entire domain produces the following correlation 

yp’s = ~vaYb. 

It is clear that w 4 = 0 means zero reaction rate, because no mixing at the molecular level 
has taken place. We conclude that the mean reaction rate can have a wide range of values 
(but limited by the upper and lower bounds) which may be radically different from the 
corresponding quasi-laminar values. 
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Next, the ease with a temperature dependent chemical reaction rate will be considered. 
We assume that k is given by the following Arrhenius form 

k{T) = k 0 ex (25) 

One parameter to measure the influence of turbulence on chemical reactions is the ratio of 
the mean reaction rate and its quasi-laminar counterpart as defined by 

R = — -T-. (26) 

k(T) 


This expression can be further written in terms of of Te/T and the pdf of T as ( BorghL 1989) 

R = ^ dTf(T)g(T:f,T E ) = g, (27) 

Jo 

where the monotonically increasing function g is defined by 

g{T) = exp{-5(^ - 1 )} = exp(5)-exp(-^), (28) 

and it has the following properties: ^(O) = 0, g(T) = 1 and 5 ( 00 ) = exp{Tg/T) > 1. Due to 
the exponential dependence of chemical reaction rate on temperature, g(T) increases rapidly 
with temperature. The rate of increase is given by the derivative of g, 

^(T)= psIT), (29) 

which has the value of Te/T 2 at the mean temperature T. Under the condition of a large 
activation temperature with a low mean temperature, strong temperature fluctuations lead 
to large ratios R 1 , because the product of g(T ) and the pdf f(T ) increases drastically 
when T > T. It is also noted that the ratio R could be less than unity, but a reduction 
in the mean reaction rate can occur only when the pdf f(T) is highly biased toward the 
low temperature side. These properties of the ratio R illustrate the effect of temperature 
fluctuations on the mean chemical reaction rate. The combined effect of temperature and 
composition fluctuations on the mean kinetic source terms will be addressed in Section 5.0 
for a turbulent diffusion jet flame. 

For high Reynolds number flows, the pdf transport equation for the set of / thermo- 
chemical variables governed by Eqn. ( 13) in the Eulerian frame or Eqn. ( 19) in the Lagrangian 
frame can be written as (Kollmann, 1989) 


(p){ 


dh 

dt 




I Q 

+ E T 1’- • • i¥>f)/l)} = 


10 



( 30 ) 


d 

dx a 

i i 


((p)K\*j = r J )f i) 


d 2 


(p) n<-, n — *})$. i )’ 


) dsp k 

j = l k=\ rj rK 

where the density-weighted pdf f\ is defined by 

J 1 = T~\ Jii'Ph ■ ■ ■ -pr-L’t)- 

(P) 

and the scalar dissipation rates in the conditional expectations are defined by 

_ dVjdVj 

e ‘ J ~ d.r a d.r Q 

with Ti = Tj = T. For homogeneous turbulent flows, this equation reduces to 


j = i ^ 3 


i i 


-«EE 


a 2 


U 1 1 a ^ t 




3i: 


:32i 


(33) 


The special case with a single scalar variable ( / = 1 ) is of particular interest. Integration over 
the scalar interval [— 00 , 9 ] leads to 


^ + 9(*)/. + £{<*iil* = *>/.}=«. 04) 

where F\ is the distribution function associated with the pdf f\. For statistically stationary 
turbulence, the time derivative can be eliminated from Eqn. (34) and the result is 


Qi'p) f xi'P ) + y {( f ll|r')/l(y’)} — 0- 
dp> 

This equatuion can be solved analytically, and we obtain the following general solution 


33) 


(«ilb)/l(*?) = Cn ■ exp{- j 
where the constant C/v is determined by 


/ Q(v') ! (36) 






TV)’ 


N' 




00 


(eub) 


(37 


11 



so that the integration of fi over [— 00 , 00 ] is unity. This solution has some interesting prop- 
erties. First, we note that the influence of turbulence on the scalar pdf / 1 (^?) is expressed in 
terms of the conditional expectation of the scalar dissipation rate {en|>?). Applying scaling 
arguments to the scalar fluctuations, one concludes that e n and p are statistically indepen- 
dent at high Reynolds numbers based on the Kolmogorov's hypotheses. If this is true, it 
follows that the pdf fi(p) and its width are determined essentially by the source Q(p) and 
(cn). A source Q(p) that is linear in p can lead to the Gaussian pdf as time goes to infinity as 
expected in a homogeneous turbulence. However, this is true only for an unbounded variable 
p. For bounded variables, such as mass fractions or concentrations, both the source Q and 
eii must be correlated with p so that the pdf is confined to its allowable domain. This should 
be true even in the limit of infinite Reynolds number. 


3.2 Nonpremixed Reacting Systems: 

For nonpremixed turbulent flames, the effect of turbulence on chemical reactions is primar- 
ily through the enhancement of mixing. Mixing is defined as the inter-diffusion of different 
components on the molecular level. The mass fraction (mole fraction or concentration) of 
component i changes due to the unbalance of diffusive fluxes passing through the bounding 
surface of a control volume. The flux defined in Eqn. (3) and its divergence in Eqn. (2) 
axe determined by the spatial gradients of mass fractions at the current time t. Hence, the 
Eulerian frame is most appropriate for the description of molecular diffusion. Chemical reac- 
tions on the other hand depend only on the local thermodynamic state, and, therefore, they 
axe best described in the Lagrangian frame. It is clear that there is no definite preference 
to the Eulerian frame or the Lagrangian frame for the description of the combined effects 
of convection, diffusion, and reactions. So far, most theoretical treatments have been based 
on the Eulerian frame, but several recent approaches use the Lagrangian or a mixed Eule- 
rian/ Lagrangian formulation. Here we will consider the diffusion process in the Lagrangian 
frame first, and the mixed formulations will be given next. 

The diffusive flux in the Lagrangian frame can be expressed as 


J'a 


Art, , AkAAAl 

2 Uo ' ^da„ da^das' 


(38) 


which contains the Lagrangian deformation tensors, dX Q /das, and the Lagrangian gradi- 
ents, d/da a . The Lagrangian deformation tensors can be further expressed in terms of the 
Lagrangian deformation rates as follows: 


dX a 

dc 


— + 



, dV 

dr- — (a, t). 

da a 


Substitution of this expression into Eqn. (38) leads to 


J a — »jw( T 


r . dv 3 

/ dr ~a — 
Jo daj) 
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(39) 


{S yuJ + 


/' 


, aV -‘, u dY ‘ 

ilr- — (a. >■))■*— 

oa^j oaf, 


This equation shows that the influence of velocity fluctuations on the diffusive flux is through 
the deformation rate histories along the pathline of a material point a. It is evident, that 
in turbulent flows, the fluctuations of flux J l a are due to continuous changes of the mass 
fraction gradients and the deformation rate tensors. Note that the mass fraction gradient is 
time independent if the mass fraction Y’, itself is a material property (i.e.. no diffusion and 
no sources). Hence, the temporal fluctuation of dYJda a is solely due to the diffusion process 
and various production sources, such as chemical reactions. If the density changes are small, 
the fluctuations of the Lagrangian deformation rate tensors are dictated essentially by the 
conservation of momentum and mass. For large density fluctuations, significant modification 
of dXa/dajj is likely to happen due to the large fluctuations in the thermodynamic variables. 

From system dynamics point of view, the evolution of turbulence can be described as a 
point in the phase space spanned by the Lagrangian position 2 l(«) and the pressure P{(±). If 
the Reynolds number is sufficiently high, there exists a region in the phase space that attracts 
all the states of turbulence irrespective of their initial conditions. This phenomenon exhibits 
several interesting properties, which also appear in the strange attractors of low dimensional 
dynamical systems. In particular, the turbulent flows are shown to be very sensitive to their 
initial conditions. Consequently, the magnitude of the deformation rate tensor, dX a /da $ , 
can be very large because material points that are initially close to each other can drift far 
apart due to turbulent motion. 

The mixed formulation of the transport equation for the mass fraction Y t is given by 


R 



dX 




'dX, 


(40) 


This equation shows that following a material point, the time rate of change of 1 , is balanced 
by the diffusion process written in the Eulerian frame and by the chemical reactions expressed 
in the Lagrangian frame. The mixed formulation reduces to the kinetic rate equation for a 
closed reacting system in the absence of diffusion. The interpretation (but not the value) of 
the diffusive term is frame dependent. In the Eulerian frame only the mass fractions Y , are 
dependent variables, and hence they are stochastic; in the Lagrangian frame, both Y , and 
X are dependent variables, and hence both are stochastic. This suggests that modeling the 
processes in Eqn, (40) can be carried out in the Lagrangian frame with ratios of random 
variables as in Borghi s M.I.L. (1988) model. 

Complete mixing is achieved if there is no scalar fluctuation. Therefore, nonzero scalar 
fluctuations indicate imperfect mixing, and they can be used as an indication of the degree 
of mixing (the mixedness) (Bilger 1976). In the Eulerian frame, the transport equation for 
the scalar variance can be expressed as following 


. ,dY" 2 , . .. dY" 2 i 


dt 


d.v a 


dx Q 


((p)vW 2 ) 
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where 


Wi 


(p r 


<3}', dYj 

' dx a dx a 


(42) 


denotes the density-weighted expectation of the scalar dissipation. As the scalar dissipation 
is always positive, it is clear that the scalar dissipation term reduces scalar variance; hence, 
it promotes mixing. This can be better described in terms of the mixing parameters defined 
as 


Y"Y" 

a ‘J = —fTrT < i <J- 


YiYj 


(43) 


It is instructive to consider a binary reacting system which is homogeneous at the macro 
scales but not at the micro scales. For such a system, the following mathematic constraints 
are satisfied 

Y" + Y” = 0 , i\ + Y 2 = 1, 
d 

— — =0, and Y\ = constant . 

OX a 


Using the first two conditions, the mixing parameters for a binary system can be written as 


yn 2 

Of 12 = “ r — , 0 < Of 12 < 1 . 

Y'dl-Yi) 


(44) 


If the two species are totally segregated (i.e., no mixing at the molecular level), a 12 is simply 

equal to one as the maximum value of Y/' 2 is Y^ ( 1 — Y\). In the opposite case of perfect 
mixing, a 12 = 0. Using the transport equation (41), a time evolution equation for 012 can 


be derived as 


<9ai2 

dt 


-^Y'iU-Fi), 


which shows that a 12 is always decreasing in time as the mixing process progresses at the 
molecular level. For reacting systems of multi-components, the values of a tJ are not simply 
bounded by zero and unity, and they may even change signs in the flow field according the 
correlations Y^Yj 9 (see Section 5.0). 


Mixing Models 

The pdf methods (Pope, 1985; Borghi, 1988; Ivollmann, 1989) offer an attractive framework 
for theoretical and computational studies of turbulent reacting flows since they do not require 
modeling of the mean chemical reaction rates. We will restrict our discussions to the single 
point pdf method. In particular, the joint scalar pdf method will be considered because 
the central issue in this paper is the influence of turbulence on chemistry. It has been 
demonstrated that the effect of turbulent mixing can be described by 


dh 

dt 


l l 


= -££ 


a 2 


}Xi kXx d P]Pk 
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which shows the joint scalar pdf is transported in the scalar space due to mixing and, more 
importantly, the shape of f x becomes narrower in time. It is worthy noting that the transport 
equation (45) represents a time-inverse diffusion process in the scalar space with the diffusion 
speed determined by the conditional expectation of the scalar dissipation rates. Consequently, 
the initial value problem posed by Eqn. (45) is difficult to solve by the traditional finite- 
difference methods as the numerical errors tend to grow exponentially. Fourier transformation 
of Eqn. (45) leads to the equation for the characteristic function m x 


(tt) 

V / mu j = l k-l 


where m x is defined by 


m x (kj)= / ••• / dy i • • • d^ t f x { y?, t ) exp( i ^ kjfj ), 

J J j=i 

and the asterisk denotes convolution. Restricting our consideration to / = 1 and assuming 
(en|v?) independent of <p, we obtain an analytic solution for f x as 

/ dkm x {k, 0) exp {ik<p + k 2 j f/r(e n )(f )}• 

This solution indicates that if mi is initially non-Gaussian, high wave-number components 
can be amplified exponentially in time. This implies that the numerical errors in the high- 
wave-number range will grow exponentially, and eventually the solution becomes unstable. 
The major conclusion is that closure models for mixing should not be based on the time- 
inverse diffusion equation. Consequently, most mixing models developed so far are of integral 
form. 

Integral Mixing Models 

Any closure model for the mixing term in Eqn. (45) should possess as many as possible the 
important features of the exact term. Strict mathematic requirements demand the model to 
satisfy the following realizability conditions: 

1. The pdf remains non-negative. 

2. The pdf remains normalized. 

3. The domain of definition of the pdf remains unchanged. (The values of pdf do not migrate 
outside the allowable domain, which are often dictated by the conservation laws.) 

These realizability conditions ensure the solution of the modeled pdf equation satisfy all the 
basic mathematic constraints of a pdf. The physics of mixing must be introduced as additional 
requirements for the closure model. The exact equation (45) can be shown (Janicka tt al ., 
1979; Pope, 1985) to reduce the higher moments of a pdf such that in the limit t — ► oc 
the Dirac pseudo-function is produced. A closure model that is able to fulfill the above 
requirements can be constructed as follows. Three basic assumptions are made here: 
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1. Mixing proceeds by pairwise interaction of fluid volumes (elements, notional particles) 
that are small compared to the volume of the flow domain. 

2. Local chaos prevails so that the probability of finding n > 2 fluid volumes with specified 
properties in a given (small compared to the flow domain) neighborhood is the product 
of the n pdf values for these properties. 

3. The time scale for the mixing process is independent of the scalar values involved in the 
mixing process. 

It can be shown that a closure model that satisfies all these three assumptions has the following 
form (Kollmann, 1989) 

(w) =7(/ R rf£'/ K ^'7i(£')/,(£")T(£'. £ ''-£) 

-/.(£)). (- 16 ) 

where 3? denotes the domain of allowable scalar space. The essential properties of the mixing 
model are contained in the transition pdf T and the time scale r. A general form of the 
transition pdf T can be cast as (Pope, 1985) 

y* 1 1 

r(£',£"->£) = J daA{a)6[<£- (1 - a)^ - -«(£' + £ )], (47) 

where A(o) is a pdf defined within [0,1]. The random variable a controls the amount of 
mixing taking place during the pairwise interaction. The construction of the mixing model 
is now reduced to the specification of A(a) and the time scale r. If we set 

A(a) = 6(o-l), (48) 

Eqn. (47) reduces to Curl’s (1963) droplet interaction model, which is computationally 
efficient but has well known deficiencies (Kollmann, 1989). Dopazo (1979) and Janicka ( t al. 
(1979) suggested 

.4(a) = 1, (49) 

which randomizes the extent of mixing and overcomes the deficiency in Curl s model. Sub- 
stitution of Eqn. (49) into Eqn. (47) yields a new form for T as 

7W' t!>) = [ I £" “V l _1 for £€[£'.£"] (50) 

— — 10 otherwise. 

Both Eqn. (49) and Eqn. (50) indicate that an equal probability is assigned to any value in 
the interval ^p n ]. This model has been applied to a wide range of flows (Pope, 1985: Jones 
and Kollmann, 1987; Chen and Kollmann, 1989, a, 1990), but it causes the higher normalized 
moments, 


(m > 4) to diverge as time goes to infinity in decaying homogeneous turbulence (Pope. 1982). 
A possible remedy for this defficiencv has been suggested by Pope (1982) as follows. If the 
probability of finding elements with the values y> at a given location is biased with the age of 
the element (time elapsed between mixing interactions normalized with an appropriate time 
scale), then bounded limits for the normalized moments / i. m can be achieved. 

Time Scales 

The time scale r for the mixing event to take place depends on the turbulent flow field and 
the scalar field. If one assumes that the time scale ratios between the turbulent flow field and 
the scalar field, 


R = 


( f n) 


( 31 ) 


are constant and independence of kinetic sources, then the time scale r can be given by 



e 


(52) 


with C denoting a constant of order unity. This assumption w 7 as found to be reasonably 
good for nonreacting free-shear flows, but it becomes questionable for reacting flows (Borglii, 
1988 ). There are several possible approaches to improve the model of time scales. First, the 
transport equation for the scalar dissipation can be included in the closure model (Lumley. 
1978 and Dibble et al. 1986 ). Therefore, the effects of heat release on the scalar dissipation 
(or on the time scale) can be included in the model (Dibble ft al., 1986 ). However, this 
approach has not been successful. Alternatively, the time scale information can be carried by 
one of the variables in the pdf. This can be done by using the two-point pdf approach (Ievlev. 
1973 ; Pope, 1985 ; Kollmann and Wu. 1987 ) or by including the scalar dissipation rate in the 
pdf (Meyers and O'Brien, 1981 ; Pope, 1989 ). Both approaches are currently under intense 
development. 
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3.3 Premixed Reacting Systems 

The theory of laminar premixed flames is well developed (Williams. 1985). Laminar premixed 
flames with a high activation energy are thin, and their structure consists of a preheat zone 
with negligible combustion and a thin reaction zone (Pope 1987). The laminar flame speed 
ui and the thickness e i can be determined entirely by the diffusivity T and the chemical 
time scale r c (Borghi, 1988). Flames with low activation energies are more complex, but 
the notions of flame speed and flame thickness are still valid. With the help of u i and 1 1 . 
the structures of turbulent premixed flames can be analyzed in the Klimov- Williams diagram 
(Borghi, 1988; Williams, 1989), where the ratio k l ^ 2 /u L over It/eL °r the Damkoehler number 


D a 


U 

k^T c 


(53) 


over the Reynolds number 


Re = 


kU t 


V 


(54) 


is plotted ( k denotes the kinetic energy of turbulence and It the turbulent macro-scale). 
Several important regimes can be identified. For large values of U/ei with low to moderate 
turbulence levels, turbulent premixed flames can be described as an ensemble of wrinkled 
flames. The effect of turbulence is essentially through the increase of flame front area per 
unit volume without significant changes in the flame structure. As turbulence intensity 
increases, pockets of fresh gases can form inside the burned product if the rate of pocket 
formation (determined by turbulence parameters) is approximately the same as the rate of 
pocket consumption (determined by chemical and turbulence parameters). This regime is 
called the corrugated flame regime (Peters, 1986). As turbulence intensity increases further, 
thick flames with their thickness larger than the micro scale of turbulence can form, and 
turbulence can significantly modify the local flame structure. 

As discussed above, the effect of turbulence on premixed combustion ranges from distor- 
tion and wrinkles of thin flame fronts to complex interactions with thick combustion zones. 
Predictions of premixed turbulent flames have been explored by Pope ( 1987) with pdf meth- 
ods, which are particularly well suited for the latter regime. However, traditional moment 
closure methods have been useful in treating the flame sheet combustion regime. Two models 
will be discussed briefly. 


The BLM Model 

The BLM model (Bray, Libby and Moss, 1985) proposed the concept that the progress of 
the reactions can be described by a single scalar variable c with certain assumptions. If the 
reaction zone is thin, then the pdf of the reaction progress variable can be approximately 
given by 

/i(c;x,f) = o(x,f) 6 (c) + fi{x_,t)8{c - 1) + 7 (55) 
with positive a,l3 , 7 and the mathematic constraint 


ck -|- 3 4" 7 — 1- 
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The coefficient at is the probability of observing reactants at the point ( x, t)\ 3 is the proba- 
bility of finding products, and 7 represents the probability of a state within the reaction zone 
(flame sheet). The most important assumption in this model is 

7 < 1. (56) 

which is valid for the thin flame sheet regime. The BLM model makes the full use of this fact 
so that the statistical moments can be expanded in terms of 7. It follows that 


0 = 1 — 3 + 0(y) 


and 


at 


1 — c 

1 + TC. 


+ 0(7), 


where r is the heat release parameter. The model consists of the transport equations for the 
mean progress variable c and its variance. These equations need to be modeled, and their 
solutions determine the local pdf for c to first order in 7. Here, only the equation for c is 
considered 


(P) 






d 

dx Q 


(( P)v'a c ") + {p)u>. 


(57) 


which contains two terms that need modeling, the turbulent flux v”c" and the mean kinetic 
source w. Detailed modeling of the turbulent flux will not be repeated here and it can be 
found in Bray, Libby and Moss (1985). Since u?(0) = u?(l) = 0, one needs to estimate the 
contribution from the continuous part of the pdf f c (flame sheet contribution). A special 
form of f c can be constructed if the reaction zone consists of randomly convected, wrinkled, 
but unstrained laminar flamelets (Bray et al ., 1985). More recently, Bray et al. (1989) have 
proposed to model the mean reaction rate directly. Two approaches have been considered 
for this direct closure method. First, the mean reaction rate is represented as the product 
of the crossing frequency of the interface at a given point and the chemical reaction rate 
per crossing. Second, the mean source terms can be treated as the product of the average 
number of flamelets per unit length (in the neighborhood of the given point ) and the chemical 
reaction rate per unit length. In both proposals, the essential issue becomes the modeling of 
the topology and geometry of an interface embedded in a turbulent flow. This aspect will be 
discussed in the next section. 


The Coherent Flame Model 

The basic concept of the coherent flame model proposed by Marble and Broad well (1977) is 
the notion of laminar flamelets, which are transported and distorted by the turbulent flow 
field, but retain their identifiable structures. This concept provides an alternative way to 
model the mean reaction rate w,. The key parameter in calculating the w, is the mean flame 
surface area density £ y defined by (Darabiha et al., 1989) 


v 


/ = 


lim 

*V'— 


bS_ 

0 6V' 
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where 8S denotes the surface element and 6V is the volume element centered at x. This 
definition deserves further explanations. First, the reaction zone is represented by a surface 
on which the instantaneous reaction rate assumes a prescribed (nonzero) value. For finite 
Peclet and Damkoehler numbers, this surface is a differentiable manifold. Let s consider 
an arbitrary point (x_,t) in the flow domain and let Ki_(r) be a sphere of radius r (volume 
element A'V) which is centered at x,. If the sphere intersects a flame surface, then we can find 
a radius r small enough so that the intersection of the surface with the sphere looks like a 
plane. Based on the properties of a differential manifold, there exist two possibilities for the 
instantaneous value of T/. One possibility is x being on the surface; then 


rrr 

±ZL r 3 



as 


r 


0 . 


The other possibility is x being outside the surface; then there exists a radius r 0 such that 


for all r < r 0 . Hence, T, f can have only the values oo and zero. It is not clear whether 
or not S / has a mean value, because the value oo must associate with the zero probability 
in order to produce a nonzero mean. However, random variables can be constructed with 
the properties of E / leading to a nonzero mean but no higher moments (Kollmann. 1989). 
Therefore, the properties of E f are dependent on the characteristics of 8S/8\ as the volume 
8V is shrunk to the point x under consideration. 

The transport equation for the mean surface area per unit volume can be derived from 
the equation for the surface element with the assumption that the correlations on the right 
hand side exist, and the result is 


(p) 


dt f 

dt 



d 

dx a 


((/»)<’"-/) - { p)n a n i )S a 3^p 


(58) 


where n is the normal vector of the surface and s a $ denotes the strain rate. The source 
term poses an intricate closure problem as turbulence can change the flame surface area by 
straining, extinction, and mutual annihilation (Darabiha et al., 1989). The proposed closure 
model for the mean reaction rate is then given by 


ih, = c Di L/, (59) 

where vq, is the volume consumption rate of species i per unit flame area and it is obtained 
from calculations for a laminar flame. Several refinements of Eqn. (59) can be made by 
including the dependence of voi on the strain rate and temperature, but they will not be 
addressed here. 
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3.4 Theory of Embedded Surfaces 

It becomes clear from the previous sections that a certain class of surfaces embedded in tur- 
bulent flow fields plays a fundamental role in determining the progress of chemical reactions. 
Experimental data obtained from mixing layers and jets indicates that the topology of sur- 
faces embedded in turbulent flows can be rather complicated (Dahm and Dimotakis. 1985: 
Dimotakis, 19S9). A closer examination of the topology and geometry of surfaces is therefore 
warranted. Surfaces can be defined implicitly by 


$(£,<) - T 0 = 0, 


where ^(jr, t) is a variable relevant to combustion. For nonpremixed combustion, 'F can be 
the mixture fraction and $<, its stoichiometric value. For premixed combustion *F can be the 
reaction progress variable and being its value at the maximum reaction rate. The variable 
*F is governed by the transport equation (2) which includes effects of convection, diffusion 
and production. Equation (2) leads to the following classification of iso-surfaces: 

(A) If the flux J' Q and the source Q , are zero, then the surface is materially invariant. If the 
velocity v a (x,t) is sufficiently smooth (at least once differentiable), then the topology of 
the surface is preserved in time, because the solution of Eqn. (2) provides a diffeomor- 
phism of the surface. However, the geometrical properties (curvature, torsion etc.) of 
the surface can change drastically. 

(B) If the flux J' a is nonzero but the source Qi is zero, then the surface is not materially 
invariant but moves through the fluid with a speed determined by the local diffusive flux. 
The Fick’s law Eqn. (3) for the diffusive flux implies that J' a is normal to the surface. 
It is possible that the topology of surfaces can be changed due to diffusive reconnection 
(Ashurst and Meiron. 1987). 

(C) If both J' a and Q, are nonzero, then in addition to the relative motion through the 
fluids, the surface can change its area due to the source term Q t . Both topological and 
geometrical changes are possible. 

(D) The surface moves with an arbitrarily defined velocity relative to the fluid. This is the 
most general case, in which both self-intersection and loss of orientability of the surface 
axe possible. Furthermore, only the progression velocity on the surface needs to be 
defined. 

Relative Progression Velocity of Iso-surfaces 

The relative progression velocity at which the iso-surface moves relative to the fluid in case 
(C) can be expressed in terms of the diffusive flux and the source strength. To this end, we 
introduce the indicator function J(x, t) 




1 for ^(x, t) > "Fo 
0 otherwise. 


(60) 


and define the relative progression velocity V(. r, f) as 

Vn a = v Q - v s a . 


(61) 
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where v Q denotes the fluid velocity, is the velocity of the iso-surface, and n Q is the normal 
vector of the iso- surface defined by 


It can be shown (Byggstoyl and Kollmann, 1986) that 

o J- 

- = -v° a n Q 


and 


01 

dx n 




Combining Eqns. (63) and (64) one obtains 


01 * 91 n 

dt +l ' a 0x a 


and 


DI 

~Dt 


= V r |V'I>|<’>('l> - ¥„). 


From Eqn. (60), we can also express DI/ Dt as 


(62) 


(63) 


64) 


(65) 


( 66 ) 


DI - T _ D<H 
— = 6 $ - ^ — 
Dt Dt 
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Eliminating DI/ Dt from Eqn. (66) with Eqn. (67). we get an expression for V as 


V = — ^-(-f^ + Q,). *(x,<) = * 0 (68) 

|V#| dx a 

with |V^| ^ 0. Equation (68) shows that the relative progression velocity depends essentially 
on the gradient of 'J'fx, t) at the location of the iso-surfaee. The treatment for case (D) is 
different because the relative progression velocity is given and an expression for u* or its mean 
value is derived. This problem has been studied by Iverstein et al. (1988) for the special case 
with a constant V in homogeneous turbulence. 

Topology of Embedded Surfaces 

Properties of surfaces that do not change under continuous one-to-one mapping are topo- 
logically invariant (i.e., a homeomorphism). We note first that the iso-surfaces that are 
interested in turbulent combustion are ‘closed’, i.e., they do not have boundaries (holes) as 
a consequence of the smooth (differentiable) variation of the defining scalars. Now let s con- 
sider the homeomorphisms defined only on the surfaces. The relevant topological properties 
are (Seifert and Threlfall, 1934; Rushing, 1973): 
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(1) Connectedness: The number of closed, connected surface components that form the 
complete surface. 

(2) Orientability: An orientable surface possesses at each point a unique normal vector 
or a unique orientation (sense of rotation). Closed, non-orient able surfaces which are 
embedded in a three-dimensional Euclidean space must intersect themselves. 

(3) Genus: The number of closed cuts that can be made without disintegrating the surface 
into disconnected parts. For instance, a sphere has a genus value of zero and a torus has 
a genus value of one. 

Surfaces can be shown to be homeomorphic to spheres with attached handles (tori) and 
crosscaps (Moebius bands attached to circular holes in the sphere). Furthermore, surfaces 
can be triangularized and then analyzed by calculating the Euler characteristic function 

\ = V-E + F , (69) 

where V denotes the number of vertices, E is the number of edges and F is the number of 
triangles. \ is a topological invariant (hence the same for the original and the triangular- 
ized surface) and changes with the number of handles and crosscaps. It can be shown that 
two surfaces, either orientable or non-orientable, are homeomorphic if they have the same 
characteristic function (Seifert and Threlfall, 1934). 

Surfaces embedded in a topological space, such as in a three dimensional Euclidean 
space, can also be classified according to their structures. This leads to the consideration 
of homeomorphisms defined for the embedding space and for the surface. Surfaces with the 
same orientability and characteristic can be further classified according to the types of knots, 
links and braids present (Moran, 1983). Links and braids may be of particular importance in 
turbulent shear flows with and without combustion. Flow visualizations in plane mixing layers 
show that braided vortical structures are the linking mechanisms between the large-scale 
structures and the lateral vortices. Furthermore, in these braided vortical structures, there 
are tube-like domains, where viscous interaction occurs. Therefore, topological change can 
take place if these braided structures are brought together sufficiently close by the convection 
process. 

Gibson (1968) pointed out the importance of extremal sets of scalars in turbulent flow’s 
and showed that they appear as isolated points (local extrema and saddle points) and lines 
(saddle lines). Consequently, Kerstein (1982) used this result to construct a model for an 
iso-surface as the boundary of the Voronoi tesselation created by the extremal points. The 
Voronoi tesselation is essentially a simplicial complex and, therefore, it is amenable to the 
methods of algebraic topology. Here, a fascinating problem arises which warrants future 
investigation: the relation of interface topology to the distribution of extremal points in 
turbulent flows. Closely related to the geometric topology of interfaces is the measure and 
dimension of such surfaces. Fractal concepts have been proposed (Gouldin, 1988; Srini\ asan 
et al., 1989) to explain the variation of surface area with scales, but the question is far from 
settled (Miller and Dimotakis, 1989). 
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4.0 The Influence of Chemical Reactions on Turbulence 

The dynamics of turbulence can be affected by chemical reactions if the heat release due to 
chemical reactions is large enough to cause noticeable density changes. This implies that 
significant density fluctuations can be expected for strongly exothermic gas-phase reactions. 
The turbulent velocity field can be modified by the chemical reactions via two processes. First, 
density fluctuations can modify the acceleration of fluid particles; that is, the light-weight 
fluid particles are accelerated faster than the heavy-weight fluid particles even under the 
same pressure gradient and the viscous stresses. Second, the pressure field can be modified 
via the state relation, which links pressure, density and temperature together. These two 
processes are difficult to model and they will not be considered here (Pope, 1985; Borghi. 
1988; Williams, 1989; Bray et al ., 1989). 


5.0 Pdf Modeling of Nonpremixed Methane Jet Flames 

The interactions of turbulence with finite rate chemistry will be examined here, particu- 
larly, in turbulent nonpremixed methane jet flames. As our current computer capabilities 
do not permit a detailed chemical scheme to be incorporated in turbulent combustion mod- 
els, simplified reaction mechanisms are needed. Peters and Kee ( 1987) developed a simplified 
mechanism for methane-air combustion based on the assumptions of a steady state for certain 
intermediate components and the partial equilibrium for two reaction steps. The simplified 
mechanism contains the following four global steps 

CH a +2H + H 2 0 ^ CO + 4tf 2 , (/) 

CO + H 2 0 ^ C0 2 +H 2 , (II) 

2 H + M^H 2 + M, (HI) 

0 2 + 3 H 2 ^2 H + 2 H 2 0 ( IV ) 

With the assumption of equal diffusivity, this mechanism requires five scalar variables 'I' _, ( j\ t) 
to determine the local thermodynamic state. They are chosen as follows: 'I'i = £ (mixture 
fraction), ty 2 = ncm, ^3 = ncOi ^4 = n, ^5 = n//, where n denotes the number of moles 
per unit mass. To explore the pdf methods, we use a combined scheme which consists of 
a Reynolds stress closure (details are in Dibble et al., 1986) and the joint scalar pdf model 
for Eqn. (30). The Reynolds stress closure provides the mean transport properties and the 
time scale that are needed in the pdf model. In return, the pdf yields the mean density. The 
turbulent flux in Eqn. (30) is modeled by a gradient type formula 




We incorporated the nonlinear interaction model to model the effects of molecular mixing 
as discussed in Section 3.2 for Eqn. (46). The effects of chemical reactions on the pdf 
are modeled by by moving the pdf position in the scalar space according to the chemical 
reaction rates of (I) to (IV) (detailed expressions can be found in Peters and Kee, 1987). 
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The boundaries of the set of realizable states in the scalar space are rather intricate but can 
be defined mathematically (Chen et al ., 1989). A detailed comparison of the first and second 
order moments with the experiments of Masri et al. (1988) can be found in Chen et al. 
(1989). The numerical solutions enable us to evaluate the mixedness parameters defined in 
equation (43) and to compare the mean reaction rates with their corresponding quasi-laminar 
rates. 

The mixedness parameter a, } is plotted at x/D — 20 for several combinations of reactants 
as shown in Fig. 1 to Fig. 3. We note first that the are not positive definite for multi- 
component (more than two) mixtures in contrast to those in a binary mixture. It follows 
from the definition of that negative values indicate that both components exceed the 
local mean, whereas positive values indicate that one of the components exceeds the local 
mean and the other component is below the the mean. Positive mixedness values indicate, 
therefore, the lack of mixing and negative values correspond to the state of good mixing. 

Fig. 1 reveals the degree of mixing between CH4 and several other components. It is 
clear from this figure that CH4 is not well mixed with the intermediate components and the 
product H2O in the inner parts of the flame, but they become mixed in the outer parts of 
the flame. Therefore, reaction (/) proceeds in the forward direction if CH4 is well mixed with 
H and H2O and in the backward direction if CH4 is well mixed with CO and H2. However, 
the forward rate is larger than the backward rate over the cross section (see Fig. 4). The 
mixedness of CO with the components involved in step II is shown in Fig. 2. It is clear that 
mixing is good for all three components in the outer parts of the flame. 

The mixedness of O2 with the active components involved in step IV is positive through- 
out the jet indicating the lack of mixing. The comparison of the mean reaction rates with 
the quasi-laminar rates for the four steps I to IV is presented in Fig. 4 to Fig. 7 showing 
clearly the strong influence of turbulence. In particular, the second step in Fig. 5 and the 
third step in Fig. 6 illustrate that turbulence can greatly enhance the average rates. On the 
other hand, turbulence can also reduce the average rate compared to the quasi-laminar rate 
as evident in Fig. 7. 
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6,0 Conclusions 

The interaction of turbulence and chemical reactions is a complex phenomenon that can cause 
significant modification in both the turbulence and the chemical reactions. For chemical reac- 
tions with negligible heat release, this interaction has only one direction; that is, turbulence 
will modify the chemical rates, but the reactions have no influence on the flow field. However, 
the purpose of combustion is to generate heat in a short period of time; therefore, combus- 
tion processes usually cause strong density variations. Hence, the turbulence will experience 
strong influence by the chemical processes and vise versa. 

The influence of turbulence on the chemical reactions is the main topic in this paper. 
First, the basic governing equations of the dynamics in turbulent combustion are outlined 
in the Eulerian frame. The corresponding transformation to the Lagrangian frame has been 
presented with the aid of the Lagrangian position field. It was shown that in the Lagrangian 
frame, the transport equation for a thermo-chemical variable does not contain the nonlinearity 
due to convection, but the diffusive flux appears as a highly nonlinear process depending on 
the time histories of the Lagrangian strain rate. 

The influence of turbulence on the chemical processes is expressed by the statistical 
moments that appear in the mean reaction rate. Upper and lower bounds for the mean 
reaction rate were obtained for a binary mixture. It was shown that the mean reaction rate 
can have values radically different from the quasi-laminar values (reaction rate at the mean 
properties). The influence of temperature fluctuations on the mean chemical reaction rate is 
shown to be significant for reactions with large activation energies at a low mean temperature. 
Furthermore, it was shown that the pdf equation for a single reactive scalar can be solved for 
a homogeneous turbulence. The solution indicates that for bounded scalar variables, a scalar 
and its dissipation rate must be correlated in order to satisify the realizability conditions. 

The analysis of nonpremixed flames illustrated the important role of mixing in deter- 
mining the progress of chemical reactions. Mixedness parameters were introduced and the 
mixing models for the single point pdf methods were discussed. The main conclusions drawn 
from the studies of nonpremixed reacting systems are that the mixing models must satisfy 
realizability conditions and that they must represent the effect of turbulent mixing correctly 
at least for lower order moments. Furthermore, the current model for the turbulent time scale 
associated with the mixing process is rather crude as it neglects the possible modification due 
to heat release. Several approaches for overcoming these shortcomings are indicated. 

The treatment of premixed systems was restricted to the flame sheet regime. In this 
regime, the effect of turbulence is essentially through the distortion of the thin flame sheet 
leading to increased flame surface area. Predictive models for premixed turbulent flames are 
discussed. 

The interaction of turbulence and chemical reactions in nonpremixed and premixed com- 
bustion can also be described in terms of the effects of turbulence on surfaces, in particular, 
the flame surfaces. Hence, the basic dynamical and topological properties of surfaces were 
introduced. It was shown that for flame surfaces, both the diffusive flux and the chemical 
sources can alter the relative progression velocity of such surfaces. Finally, the topological 
classification of embedded surfaces was discussed briefly. 

We used nonpremixed turbulent methane-air jet flames as an example to illustrate the 
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mixing properties predicted by a pdf method. It was shown that the mean reaction rates 
can be larger or smaller than the quasi-laminar rates by orders of magnitude due to the 
effect of turbulence. From the computation results, the predicted mixedness parameters were 
examined showing that they are not positive-definite as in binary systems. 
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Fig. 1. Mixedness parameter otij in a turbulent methane- air nonpremixed flame at x/D 
20 using the four step mechanism of Peters and Kee (1987) for: 

1 = CH A and H. 3 = CJ7 4 and CO. 


2 = CH* and H 2 0. 


4 = C and H 2 . 




Fig. 2. Mixedness parameter a tJ in a turbulent methane-air nonpremixed flame at x/D 
20 using the four step mechanism of Peters and Kee (1987) for: 

5 = CO and H 2 0. 6 = CO and H 2 . 


7 = CO and C0 2 . 




Fig. 3. Mixedness parameter a, ; in a turbulent methane-air nonpremixed flame at x/D = 
20 using the four step mechanism of Peters and Kee (1987) for: 

8 = 0 2 and H 2 . 9 = 0 2 and H 2 0. 



Fig. 4. Mean kinetic source (full line) and quasi-laminar source (broken line) in a turbulent 
methane-air nonpremixed flame at x/D = 20 for reaction I of the four step mechanism of 
Peters and Kee (1987). 
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Fig. 5. Mean kinetic source (full line) and quasi-laminar source (broken line) in a turbulent 
methane-air nonpremixed flame at x/D = 20 for reaction II of the four step mechanism 
of Peters and Kee (1987). 
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Fig. 6. Mean kinetic source (full line) and quasi-laminar source (broken line) in a turbulent 
methane-air nonpremixed flame at x/D = 20 for reaction III of the four step mechanism 
of Peters and Kee (1987). 
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Fig. 7. Mean kinetic source (full line) and quasi-laminar source (broken line) in a turbulent 
methane- air nonpremixed flame at x/D = 20 for reaction IV of the four step mechanism 
of Peters and Kee (1987). 
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1.0 Introduction. 


Kraichnan’s idea to apply mappings as tool in constructing closures for pdf equations ( Chen 
et al. 1989, Kraiehnan 1990, Feng 1991, Pope 1991. Valino et al. 1991) proved very successful 
for the case of a single scalar variable in homogeneous turbulence. It is not yet clear how 
powerful this approach is for the case of more than one variable (Pope, 1991). The rigorous 
theory on the functional level requires the determination of the probability measure governing 
the dynamics of turbulent flows. This measure is defined in a function space, that contains 
the set of all realizable flow fields. It can be constructed as a measure relative to a Gaussian 
measure or any other reference measure. It follows that the probability measure can be 
regarded as the image of a Gaussian measure. The underlying mapping contains, therefore, 
the essential information on turbulence. It is apparent that turbulence can be viewed as the 
mapping of appropriately defined function spaces. The structure of those function spaces and 
the mapping relating them deserve closer scrutiny. 

The notion of a mapping can be exploited in two distinct ways: 

I. The mapping is known and both the original and the image variables are the unknowns. 
The prime example for this case is the linear map provided by Fourier transform which can 
be extended to functionals. The resulting equation for the case of the probability functional 
is the Hopf/ Kolmogorov equation for the characteristic functional. 

II. The mapping is unknown and either the original or the image variable is known. The 
prime example for this case is Kraichnan’s method which is based on the well known relation 
between mapped and original pdf depending on a single variable 


dX 

dr} 



where p> = X(rj,t) and fdv) * s the Gaussian pdf. The extension of this method to the 
multi-dimensional case and to functionals is not obvious. The task of constructing a mapping 
between function spaces can be daunting and it is instructive to consider several examples to 
illustrate the properties of mappings. Fourier transformation, regarded as solution of the so 
far unspecified mapping equation, provides an example for the mapping of the space L~ of 
square integrable functions defined in R 3 onto itself 

^(£) = y <*£exp(i>.r)$(x) 

ft 3 

where E R 3 and $(r), , P(_r) E L 2 (R 3 ). The image field ^(c.) is a complex valued square 
integrable function defined on R 3 . The value of the image field $ at a location z_ E R 3 
depends on the values of the argument field $ at all locations r E R 3 ■ Hence has the 
mapping induced by Fourier transformation functional (or nonlocal) character. Mappings of 
this type are denoted by 


o 



where the semi-colon separates functional arguments from parameters. A different example 
is provided by 

T(c) = exp($(r)) 

which has local character. Locations in argument and image fields are the same and a change 
of the argument field <&(c/ ) at a location c/ 7^ r. has no influence on the value of the image 
field at z_. Local mappings are denoted by 

).;,<) 

where 

Y : R* — R 3 

is a mapping of R 3 onto itself. It is clear that the computational effort for functional and local 
mapping can be expected to be widely different. The comparison of a functional mapping 
with the probability functional indicates that the image domain for the mapping is a function 
space whereas the image domain for the probability functional is the unit interval. This 
indicates that the computational effort for the calculation of the mapping may be equal or 
larger than the effort for the probability functional. However, it should be noted that the 
Gaussian characteristic functional can be set up explicitely and the notion of the determinant 
can be extended to the case of countable infinite many variables (see Muldowney, 1987 and 
Skorohod, 1974), thus offering an avenue for theoretical investigations. 

Three aspects of turbulence involving mappings will be discussed in detail before the 
mapping equation is investigated. First, the basic laws (mass and momentum balances) are 
set up as equations determining a mapping of the flow domain at a reference time onto the 
domain at a later time. This mapping is a diffeomorphism as long as the smoothness of the 
solution of the Navier-Stokes equations is insured. In addition, incompressible flows generate a 
measure preserving diffeomorphism due to the particular form of mass balance. This mapping 
is, however, fundamentally different from the mapping employed by Chen et al. (1989) who 
map the phase space spanned by the solutions of the Navier-Stokes equations onto a reference 
space equipped with a Gaussian measure. Second, the properties of Gaussian measures are 
reviewed and finally transformations of the function space containing the solutions of the 
Navier-Stokes and Euler systems are discussed. 

The mapping method suggested by Chen et al. (1989) is the reviewed for the one- 
dimensional case. It is shown that it corresponds to a convolution of the characteristic 
function and the mapping equation for it is derived. The multi- dimensional case is then 
considered in detail. 
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1.1 The basic laws as mapping equations. 

The basic laws for a single incompressible fluid are introduced in the material (or Lagrangean) 
frame. The independent variables in the material frame are defined as time 0 < t < T and 
label a £ A, where A is the label space to be defined. The label identifies uniquely a material 
point in the flow field. There are many ways of defining a label and for each definition a 
different set of variables emerges. The present definition for the label space A is the position 
of all fluid material points at a reference time zero. The label space is now the fluid volume 
at the reference time zero 


A = {X. : Position of a material point at t = 0} 


( 1 . 1 ) 


and the mapping X_ : .4 — ► D{t) is assumed to have the following properties (Kreiss and 
Lorenz, 1989): 2L{o,t) is for a £ A and a nonzero time interval a smooth (continuously 
differentiable) function satisfying 

(i) 0) = a 

(ii) If a ^ b then X_ {a,t) ^ X(b* t) for t > 0. 

(iii) The mapping a — > Xfa, t) has a smooth inverse X ~ l . 

The set of independent variables in the material (sometimes also called referential frame for 
the present chioce of the label space) frame consists therefore of time t and the position 
a) at time zero. Note that the size of the time interval for which the mapping remains 
smooth is related to the existence of smooth solutions of the Navier-Stokes equations. There 
is no general existence proof for three-dimensional flows except for finite time intervals, whose 
length depends on the inital data, and we must keep in mind that the smoothness may break 
down in finite time and the conditions (ii) and (iii) for the mapping of the label space onto 
the fluid volume at later time may be violated. The position Xft.a) of a material point is 
considered a dependent variable in the material frame. 

The time rate of change keeping the label variable a constant is in the material frame 
the time rate of change measured by an observer moving with the material point. It follows 
that velocity is defined by 


dX 

s iit 


( 1 . 2 ) 


and acceleration by 


dV 


(1.3) 


Mass balance in classical mechanics is the requirement that mass cannot be created or de- 
stroyed and appears as .7 = 1 or 


dX Q dXj dA' 7 

Catf-f 


(1-4) 


where J denotes the Jacobian of the mapping A_( a, t). Momentum balance is the consequence 
of Newton s second law and can be shown to govern the temporal evolution of the mapping 
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X(a,t) by 


d 2 X a 1_ dX^dX^dP^ v dX < dX 0 0 , dA\- dX 0 d*X a ^ 

dt 2 2R ° J ' 1 6r> ' J da n da ^ das 2 fl9J d aj? do,^ da* d a-> da^dt + a 

(1.5) 

It follows that the material frame version of the Xavier- Stokes system set up in terms of po- 
sition and pressure fields can be regarded as the equations determining a measure-preserving 
mapping A_(a, £) : £)( 0) — ► D(t ) (where D{ 0) denotes a measurable subset of the flow do- 
main at time zero) of the flow domain at the reference time onto the domain at a later time. 
This mapping is smooth as long as the solution remains smooth. The role of the pressure in 
this system is to preserve the measure of any measurable subset of the flow domain. This 
property is lost in compressible flows where the volume of a materially invariant subset of 
the flow field may change in time. The notion of mappings can be applied to compressible 
flows with suitable extension of the variables to include mass and energy densities. 

The mapping equations (1.4) and (1.5) show’ that the dynamics of a nonlinear phe- 
nomenon can be viewed as the evolution of a mapping governed by a nonlinear system of 
equations. The solution of (1.4) and (1.5) allows to determine the image of any measurable 
subset of the original flow domain as a function of time. 

2*0 Basic Considerations. 

The notion of mappings can be extended to relations between function spaces and a general 
equation for such a mapping can be derived. It is worth noting that the mapping equation 
for a single variable governed by the linear diffusion equation and taken at a single point is 
linear (Chen et al., 1989 and Pope 1991) but not exact. It will be shown in chapter 4.0 that 
this is due to the requirement that the single point statistics of the image of the Gaussian 
random fields is equal to the single point statistics of the unknown turbulent fields. The 
extension of this closure procedure to multi-variables and multi- point pdfs and the passage 
to the functional level requires some preparations. This will be done in the present chapter. 
It is well knowm that the Lebesgue measure has no extension to infinitely many variables, 
hence there is no obvious extension of the cdf equation (derived in the next chapter) to the 
functional level since it contains multi-dimensional integrals whose limit for infinitely many 
variables is not defined. The pdf equation on the other hand can be (at least formally) 
extended to the functional level. It follows that the pdf equation and its Fourier transform 
(characteristic function) are the appropriate starting point for the development of mapping 
methods for the multi-dimensional case. Furthermore, it is possible to set up explicitelv the 
Gaussian measure for infinitely many variables in terms of its characteristic functional and 
there exists a well defined equation for the characteristic functional of the turbulence measure. 
It is possible to define properly measures relative to a Gaussian measure (Skorohod, 1974). If 
the turbulence measure is defined in such a way, all that is left to determine is the distortion 
of the Gaussian measure necessary to produce the turbulence measure. This is nothing but 
a mapping of the function spaces containing the Gaussian fields and the turbulence fields. It 
is clear that these two function spaces justify attention. 
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2.1 Function spaces. 

Consider a compact flow domain denoted by D(t) C /?* with boundary dD(t). The boundary 
is assumed to be orientable and sufficiently smooth such that a normal vector exists nearly 
everywhere. The surface area for the boundary OD must satisfy 


0 < 


/ 


dA < do 


and the volume of the flow domain is obviously bounded. D(t) is the domain of definition for 
the turbulence fields. The domain of definition for the Gaussian reference fields is D = i? 3 . 
Cartesian coordinate systems are introduced in both domains and denoted by r € D(t) and 
i € D respectively. Various functions of scalar, vector and tensor character will be defined on 
D and called turbulence fields. Functions defined on D will be called reference or argument 
fields. Various sets of fields will be considered and they will be embedded in various function 
spaces. The most important Banach and Hilbert spaces will be discussed briefly. A class of 
Banach spaces is given by 

L P (D) = {$(r)|$ : D — ► measure able, J drj<£(r)| p < oc} 

D 


with norm 

11*11 1.(D) 3 { f d t \Mx)\')i 

D 

which is a Hilbert space for p = 2 with scalar product 

($,'!') = J dr$(r)^(r) 

D 

A different class of function spaces can be constructed by requiring differentiability up to a 
certain order. These spaces are called Sobolev spaces and are defined by 


W m p {D) = {$(r)|D Q $ € L p ,\a\ < m} 
where D' = d/dx t denotes the differential operator and 

D° = 

“ dx°' • • • dx a n n 

with |a| = a » a i ^ 0. The norm is defined by 


||*||m,,D = E 

| a | < m 

For p = 2 a Hilbert space is obtained with scalar product 

($,*)„, d = f <lxD a *D Q * 

\ a \<m J D 

The general properties of these spaces can be found in the literature on functional analysis. 
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2.2 Local and global mappings. 

The aim of this chapter is to discuss the properties of mappings of a separable Hilbert space 
H = {$(£)!$ : D — + continuous} with scalar product onto another separable 

Hilbert space H (separable means that the space has a countably infinite basis). The elements 
of H axe the reference or argument fields and the elements of H are the image fields. The 
basic requirement for mapping methods is that the statistics of the image fields agree at 
least at one or several points with the turbulence fields. If they agree on all points of the 
flow domain an exact solution of the Hopf-Kolmogorov equation (see Vishik and Fursikov. 
1988) is obtained. This aspect will be discussed in chapter 4.0. The space H containing the 
image fields is embedded in the Sobolev space H m (D). The mapping A : H — ► H(D) will 
be time dependent because the image fields must be time dependent in order to simulate the 
properties of turbulent fields. The mapping is denoted by 


$(£,/) = A'[$(.);x,t] 


where the semicolon separates argument fields from parameters. The argument field 4*(r) G 
H, the values of the parameters x G D(t) and t > 0 determine uniquely the image field 
$(x, t) G H(D ) and the argument field, the location in the flow domain and time can be 
varied independently. Consider now a modified argument field ^(r) + — x_ 0 ) where 

hs £ H , e > 0 and 


M£-£ 0 ) 


> 0 for |x — xj <6 > 0 
0 otherwise 


and h £ C°°(R 3 ) such that 

J d£M£ - £J = 1 

R 3 


holds. The image fields at the same location x and the same time t are denoted by $(x . t ) — 
X[4(.);r, t] and t) = *Y[4(.) + cht(.)\ r, #]. Two cases are now possible: 

(A) *(x,<) / *(£ , t) for nearly all r £ Cl and e > 0. The value of the image field at a 
given location x_ G D(t) depends on the values of the argument field at nearly all locations 
r G fL This implies that X depends on <$(.) in functional fashion. 

(B) There exists a subset Cl a C Cl such that '{'(r, t) = $(r, t) holds for x_ £ Cl 0 an< l 

< C8 3 with C < oo as 8 — * 0 for all r G D(t). This implies the existence of a 

relation 

£ = Y(L,t) 


such that the mapping X is local, i.e. 

4>(r,t) = A'($(r(r, t));x,t), x £ D(t) 


The special case of a local mapping 


$(£,<) = X(MY(x,t));t) 


is called r-autonomous and 

$(x,#) = .Y($(r(x,t))) 

is called fully autonomous. 
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2.2.1 Properties of global maps. 

Global maps are characterized by the property that the change of the argument field in a 
small neighbourhood of any point € D leads to a change of the value of the image field 
at a fixed location x_ € D. This functional dependence indicates that it must be possible 
to express the mapping A' in terms of a functional. This can be achieved using a special 
construction which leads to a subset of the set of global mappings. Consider an arbitrary 
functional 

A[$] : H(D) -► R 1 

The value A[4>] is then independent of i £ D and derivatives with respect to location vanish. 
Suppose now that A is Frechet-differentiable in H 


r 4 ei 

( — , h) — lim — A[$ + eh] 

6*(i) t-ode [ j 


for <5, h € H . Then is the first Frechet-derivative of A a generalized function of the location 
i. It follows that 

$ A * A 

TJ— =.Y-[#(.) :i ], «€/f 

provides a mapping of the argument space H into some function space H*(D) because keep- 
ing the argument field fixed and varying the location x G D produces a scalar field whose 
smoothness properties depend on the functional A. A second step is required for the construc- 
tion of global mappings based on functionals. We need to define a mapping Y_ : D{t) — ► D 
which is bijective, local and sufficiently smooth. Then we can express the location in D(t) in 
terms of the location in D by 

i = Y_(lJ) 

If the functional A is defined in such a way that the fields generated by its Frechet-derivative 
are always contained in the space H m (D) for m > 0 and combining the derivative with the 
mapping F(r, t) of the domains we can define a mapping by 


_ 6A[4( . 

= my(xj)) 

It is clear that global mapping do not require a mapping of the domains and for this reason 
is the present construction rather special. 


-v[4( 
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2.2.2 Properties of local maps. 

The value $(x, t ) of the image field depends on the value of the argument field at a unique 
location x = F(x, t). If the argument field is modified at any other location no change of 
$(x, £) is observed. Hence is 


$(£' t) = X( $(F(x, t ). t ) 

a strictly local map. Differentiation of the image field can be carried out using the rules of 
standard calculus. For instance, the time derivative emerges as 

<9d> , dX <9$ d¥ a dX 

dt ” d$ dYa dt + dt 

and likewise for the spatial derivative. Higher derivatives follow from repeated application of 
these operations. 


2.3 A reduction property for multi-dimensional mappings. 


The closure problem for mappings can be viewed as the construction of a global map corre- 
sponding to the local map which is to be determined as the solution of the mapping equation. 
The relation of the global map to the local map needs clarification since they act on the 
same class of reference and turbulence fields. If the global map is known we can calculate 
the statistics at any number of points as the image of of Gaussian statistics. Hence, we can 
calculate the statistics at a single point which implies that there exists a relation between 
the mapping for a single variable and the mapping for many variables containing the single 
point. This relation can be obtained as follows. Pdfs posses a well known reduction property 
given by 


fi(VN) 



OO 

J dpN-ifMVl, ' ' ’ iVn) 

— OO 


which must be retained if the pdfs are transformed Gaussians. Let the local mapping be 
X : R 1 — > Q, where denotes the range of the scalar defined at a single point in the 

flow field D(t), and the global (N-dimensional) mapping X V : R N —* f!‘ v , then are the 
one-dimensional pdf and the N-dimensional pdf given by 


fiy) = 


fcin) 

a.v 

dt) 


and 


/n( j/i i ■ ■ ■ i y n ) - 


/g ( * 7i ' ■ ■ ■ • *?iv ) 
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respectively where y = X(y) and y = X[y hold. Application to the reduction 
property for pdfs leads to 


/ g ( y ) 
ax 

dr) 



30 



where the subscripts on the Jacobian matrices indicate the rank. Since the Gaussian reference 
measure is the product of N one-dimensional measures it follows that 


dX~ l 

dy 



N - 1 


n fG(Vi) 


1=1 


holds. Denoting by dG, = dy , fo(fji) the differential of the standard Gaussian measure (zero 
mean and unit variance) we get 



as reduction property for mappings. It is worth noting that this relation can be extended to 
infinitely many variables since the Gaussian has a well defined limit. It is easy to show that 
the reduction property appears in the form (note that y = y ,v) 


(^r) = I dG{Ih) "I 

— oo — oo 

if the Jacobian matrices are triangular. This particular form of the reduction property 
expresses the one-dimensional map for the N th variable y t v in terms of the global or N- 
dimensional map for all variables rji, • • • . yy . It depends obviously on the Gaussian reference 
measure. 
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3.0 Transport equations for pdfs and cdfs. 

Transport equations for pdfs and cdfs can be obtained using the notions of step functions 
and pseudo- functions. The starting point are the basic differential equations set up in the 
spatial (Eulerian) frame with respect to a Cartesian system of coordinates. Mass balance is 
given by 

dv a 


dx. 


and momentum balance by 


dVg dl'c, 

dt +VJ dxj 


= 0 


1 dp d 2 v 

+ v 


:3.i; 


P dx a 


A passive scalar obeys 


dxjdxj 

d 2 $ 


d<S> d<$> ^ „ 

tjt — — Q($) + r _ „ 

ot OXj OX iiOX3 


!3.2) 


: 3.31 


The flow domain is denoted by f2 and its boundary by dQ. Consider now N distinct points 
g ft. We derive first the cumulative distribution function (cdf) of the values of velocity 
ans scalar at those N points in the flow domain. We define the step function 

Fv = n S^Hivi - *(z li \t))H(vi-v{x li \t)) (3.4) 

where N denotes the number of points in the flow field, H denotes the Heaviside function 

H(x)= { 


1 for x > 0 
0 for x < 0 


and the step function of a vector argument is the product of the step functions of the com- 
ponents. It is easy to show that the expectation of Fv is the cdf Fv 


(F;v) = F/v(tL!, ■ • • , , v^v ;^. 1 


(i) 




: 3.5) 


of the variables v(x}'\ 1), $(x ( d, t) at the jV points in the flow field. The necessary tools for 
the derivation of the transport equations are the spatial and temporal derivatives of the step 
function Fv- Implicit differentiation leads to 


OF 


N 


dt 




and 


OF 




dFy 


dXa ) Ola’ dVj dx a 


(J) 


(x} } \t) 


d<Pi 

dP N 


! 3.6) 


(3.7) 


where the summation convention applies to Greek subscripts. 
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Pressure. 

The pressure p(x_,t) can be eliminated from the momentum balance according to 


— A p = 


dv a dr,) 
dx y dx Q 


with boundary conditions 


dp_ 

dn a 


= h a {v) 


on dSl. The solution of this Poisson equation for p(x_, t) can be given in the form 

dv a dry 

a yt-r (£_,!£, 

n 


P( 


£,,)= -h j iy - GU ' y - ) wM +B{ -' ,) 


:3.S) 


where G(x, y) = |x — y\~ x is the Poisson kernel and B(x_,t ) is a harmonic function ensuring 
that the boundary conditions are satisfied. It follows from this solution that pressure depends 
in functional form on velocity, it depends on time in autonomous form via velocity and it 
depends on location x parametrically. Hence p(ji,t) = p[u(., f); .x] and 


dp dv 

and the Frechet- derivative of pressure with respect to velocity is given by 


(3.9) 


S_ 

bv a (y,t) 


pbd-Jhi :} = - 


1 d 2 G(x,y) 

2tt dy a dyy 




if no boundaries are present. This derivative represents the change of pressure at x due to a 
change of the velocity component v Q at y. 


Transport equation for Fv. 

The transport equation for Fs follows from the differentiation rules and the Navier-Stokes 
system. We get 


dF N 

dt 


N 

i=i 


dF\ 

dXa ] 


N 

+ ]T{rA (,) $ 

i=\ 


dF„ 

dy> t 


+ i/A (,) t> 0 


dF s 

dv' a 


} 




;=i 


dFs 6 

d<p t dx 


p[£(..<);£]|~} = o 


(3.10) 


This equation can be viewed as the condition that the value of the cdf Fs remains constant 
for points in phase space (the product of flow domain, scalar space and velocity spaces) that 
move with the velocity (v 0 (x*‘\ <), TA^$ + £?($), uA Xt ^v a — dp/dx a ). 
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Preliminary form of the transport equation for the cdf Fv. 

Averaging of the transport equation for Fy leads at once to the equation for the cdf Fy. Its 
preliminary form is given by 


dFy 

dt 


N 


i = l OX a i = l X>1 


N 


dFy 


dFy 
dv ‘ 


)} 


,v 




1 = 1 


( 3 . 11 ) 


All expectations appearing in this equations need to be evaluated to provide expressions 
containing the cdf F*. If k < N the term is called closed otherwise it is nonclosed and 
represents an additional unknown. 


Flux in physical space. 

Mass balance holds at every point and this implies 

= -4 T (o<.(£ <i ’,()CN) 

U X a UX (y 

The flux is now according to the definition of the angular brackets given by 

(v a (x}' ] ,t)Fy) = 

J ■ • J dv(F N) )v a {x (,) )Fyfy(v{x {l] ), ■ ■ ■ , $(r (iV) ); r (1) , • • • ,x ( *Kt) 

The definition of F/v implies 

iV 

<P 1 v 

(v a ( xit \t)Fy) = j </$!••• J dv»v' a fy(i 


The relation 


leads to 


, d* N Fy 

* N ~ ^,,1 


dv\ ■ ■■d$ N 

t i dFy 


(v a (x U \t)Fy) = J dv' Q v' a -I^ 

— oo 

and partial integration results in 


(v a {x^\t)Fy) = v l a F t 


u a 

N ~ j dv' a 

— OO 


Fy 


( 3 . 12 ) 
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which is closed. It can be shown easily that it is consistent with the corresponding expression 
in the pdf equation. 


Diffusive flux in scalar space. 

The flux in scalar space is given by (A which is by definition 


(A (x) <t>F N ) = dtp J d{A^‘^) 


-£ ; )/:V+i($\--,£ j ,A (,) <M) 

Introducing the conditional pdf f c defined by 


leads to 


(A (,) $Fiv) = J dQ 1 

— oo 


■ J dv N { A (,) $|$(£ (1) ) = $V- 

— CO 


£(V N) ) = £' V ) 


cH v F v 
. . . 0)£ V 

(3.13) 


as final result. A new unknown (A^'^l • • •) appears and the flux in scalar space is, therefore, 
nonclosed. Consistency with the pdf equation follows at once from the equation above by 
differentiation. It is important to notice the fundamental difference between the cases A = 1 
and N > 1. The term actually appearing in the transport equation for the cdf is 


A(A {i >*F N ) 


which is for .V = 1 obviously a differential term 


d 

dp. 


\ Pi 

(A = J d*'(A {l) $Wx i ' ) ) = <l> t )^ 


given by (see Pope, 1991) 


JLiA^QF:) = <A (l) $|*(x (,) ) = ^,)|4 

0<P , Oip 

whereas for N > 1 the integral character is preserved since the single derivative removes only 
one of the integrals. 

Viscous flux in velocity space. 
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The viscous flux in velocity space is completely analogue to the diffusive flux in scalar space. 
It follows from (13) that 


(A { ' ) v a F s ) = J 

— OO 



holds. The flux is nonclosed and consistent with the pdf equation. 


(3.14) 


Source flux in scalar space. 

The source term (?($,£) is assumed to be a local function of the scalar and the 

velocty v(x}'Kt) It follows from (3.13) that 

= j <»'■■■ j 

— oo — oo 


holds. This expression can be simplified if the source term Q depends only on and r 1 
because partial integration removes the differentiation with respect to all the other variables 
and we get 

w n' 

(QF N ) = J dV J d £ ’Q(* l ,v 4 )-^ (315) 

— oo — oo 

which is a closed term. For the special case that the source Q depends on only we obtain 
for the derivative appearing in the cdf equation 


d 

&V>i 


(QF N )=Qivi) 


dFy 

d<pi 


Consistency with the pdf equation becomes evident by differentiation of (3.15). 


Pressure flux in velocity space. 

It follows from the functional dependence of p on the velocity v_ that the pressure gradient is 
a nonlocal function of v_. Hence we can apply (3.13) and get 



(3.16) 


This result is nonclosed. 


Cdf equation in terms of conditional fluxes. 
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The expressions (3.12) to (3.16) can be used to obtain a new version of the transport equation 
for the cdf Fy 


dF* 

dt 


A r , <9 F v f , dF w , 

+ S j ' v axL" f 


1=1 

;/ V 


N V?1 — 

f d* 1 -- f dv y { a ,, ' ) $|$(£ ,1) ) = *',••• 

1=1 Vi J d 


<9 4 - v F.v 

as 1 •••av‘ v 


— oo 
,v 


+l, M f <'*'•-/ <fe‘ V (A l "l'a|*(£ (11 ) = 


, v v?i n 

*»&/ d<F' J dv l <?(*•.£ 


. .... a 4 ** 


il 

7 / - = *',•• > a|-^N ) = 0 < 317 > 


d 4N F N 


dv ' 0 

This version is not suitable for the limit iV — ♦ oo since it contains the A^-dimensional Lebesgue 
measure for which no limit exists. However, the transport equation (3,17) has several note- 
worthy properties. First we note the special case of a single scalar in homogeneous turbulence 
with £?($). It follows that the cdf for the single scalar satisfies 


dF dF 

— + { T ( A ^=^ + Q (^)}—=0 


(3. IS) 


which is non-integral with respect to the scalar value <p . The hyperbolic equation (3.18) 
expresses the fact that the value of the cdf remains constant for points moving with the 
velocity r(A$|$ = tp) + Q{p). We conclude that the structure of the cdf equation for A = 1 
and N > 1 is fundamentally different and methods developed for N = 1 cannot be expected 
to carry over to N > 1 without major modifications. 

Pdf equation in terms of conditional fluxes. 

The transport equation for the cdf can be used to obtain the equation for the pdf /.v by 
differentiation according to 

d iN F N 

/ v ~ dv\ • - • <9<$ iV 


It follows that the pdf is governed by 

df N 

dt 


,v 


* d/,y 
+ ^'' Y ar (,) 

1=1 UX Q 

16 

a- A 



(3.19) 


N d 1 

V{r — [(A ( ‘ ) $|$U< 1) ) = • ■ -)/.v] + u — [(A (l, r Q |$(x (1) ) = • • ->/.v 

l=1 0/1 a 


]} 




which is formally of local (non-integral) character in scalar- velocity space. It reduces for a 
single scalar in homogeneous turbulence to the well known equation 


^ + ^{[r<A$|* = *)+<?(*)]/} =0 (3.20) 

which can be regarded as the condition that the divergence of a flux in the phase space 
spanned by time and the scalar space is zero. Equation (3.20) and its generalisation to higher 
dimensions (3.19) can be regarded as balance equations for the ’’mass” per unit "volume 
fw. The coefficients of / can be interpreted as velocity components and fs will remain 
nonnegative and its integral unity. 


4.0 Mapping method for the one-dimensional case. 

The mapping method suggested by Chen et al. (1989) has two distinct advantages over the 
previously constructed models. It produces the Gaussian pdf as asymptotic limit for decaying 
turbulence and agrees very well with DNS results for pure mixing in homogeneous turbulence. 
The derivation of the mapping equation will be reviewed and its variant for the characteristic 
function will be discussed. 
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4.1 The Chen-Kraichnan method. 

The starting point for the Chen-Kraichnan method is the pdf equation for a single scalar 
$(x, t) given by 

I + <’•>& + - l/ {rAm = M) (41) 


which reduces for homogeneous turbulence to 


¥ + |(^ |=o 


The equation for the corresponding cdf is according 


dF „ 
dt +F *dip 


= 0 


The velocity in scalar space is defined by 


F^ = Q + (rA$|$(x,t) = <p) 


(4.2) 


(4.3) 


(4.4) 


It governs the dynamics of pdf and cdf. The probabilistic argument space of pdf and cdf 
is the range of values the scalar $ can assume at any point (x, t). Let this space be the 
unit interval [0, 1]. There are two slightly different ways for the development of a mapping 
method. 

Method I.: A mapping X : R l — ► [0,1] is defined by two conditions. First, the value of 
the cdf at the image variable X(rj) is equal to the value of the standard Gaussian cdf at the 
argument variable q 

F(X(r],x,t)) = F G (r)) (4.5) 

and second, the mapping is monotonically increasing 


X(rfi) < X(t} 2 ) for 2 ( 4 - 6 ) 

Mapping the domain of definition of the single point cdf implies that the value of the scalar 
$(x, t) is the image of a variable ranging on (— 00 , 00 ). Hence will the mapping in general 
depend on the location in the flow field and on time. The first condition (4.5) implies the 
existence of a Gaussian random variable 'I' such that 


F g (tj) = Probi^ < *?} 

holds. The argument variable is now extended to a Gaussian random field such that the 
Gaussian random variable is the Gaussian field at a location £ = !_(£, t). Note that this 
extension is completely arbitrary and the mapping X does not determine the relation between 
the locations of the image variable $(x, t) and the argument variable 'J'(C). The second 
condition (4.6) is a direct consequence of the monotonicity of the cdfs. 
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The probabilistic interpretation of the first condition (4.5) using 

F(X(q,x,t)) = Prob{$(x,t) < A'(//, x,t)} 


and 

is given by 


F G (n) = ProbiViC) < t /} 


Prob{$(x_, t) < A’(r/,x, t)} = Pro6{^(C) < r ]} 
which can be regarded as 


Prob{$(x,t) < 9 } = Prob{X{*(Q,x,t) < <p) (4.7) 

It should be noted that the mapping thus constructed is local in the sense that any change 
of the Gaussian argument field 'J'(C) at any location £ ^ !_(x., t) has no influence on the 
mapping. We note that there are three fields involved in this version of the mapping method: 

1) The turbulent field 4>(x, t). 

2) The Gaussian argument field 'J'(C)- 

3) The (local) image field <£(.r, t) = X(\&(C where £ = K(x, t) holds. The image 
field is called surrogate field (Pope, 1991). 

It follows from the fact that the mapping relates only single point statistics to Gaussian 
statistics that no scale information will be determined by this version of the mapping method. 
The mapping Y_ relating the locations of turbulent and argument fields is so far undetermined. 

The mapping equation is now derived from the first condition (4.5) by differentiating 
with respect to time. It follows that 


d¥_ dX dF 
dt dt dip 


(4.S) 


must hold since the reference distribution is time independent. 
(4.4) leads to 


^ = Q + *) = ¥>) 


Comparison with (4.3) and 


(4.9) 


The closure model is completed by requiring that the conditional expectation of the turbulent 
field is equal to the conditional expectation of the image of the Gaussian reference field 
(surrogate field) 


|A=g(9) + (rA$|$(x,#) = 9) 


(4.10) 


The most important property of this mapping equation is the fact that the Gaussian 
reference measure allows explicit calculation of the conditional expectations in terms of the 
mapping and correlations of the Gaussian reference field. The derivatives of the image field 
can be expressed in terms of the reference field as follows 


dx a 


r\ 

—xmru,t),t) 

OX 
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and thus 


d$_ = d^d^_m 

dx a dip dQj dx Q 

holds. If the so far undetermined relation between the locations in physical and reference 
spaces Y_ is restricted to a stretching transformation uniform in physical space we get 

| — =8 aij m{t) (4.12) 

ux a 


where m is time dependent. Denoting the derivatives of the mapping with 



a A' a a' 

(4.13) 


A"=— , -r— = A 

dip dt 

we get 

d$ , ^ dY a 

(4.14) 


dx a " do dx a 

and for the Laplacian 




2 ( Y I,0* W , a 2 * \ 

dx a dx a m \ d( Q d( a " dCadw 


(4.15) 


The conditional expectation of the Laplacain can be established explicitly in terms of the 
mapping X and correlations of the Gaussian reference field if and only if the derivatives 
of the mapping are completely specified by the condition $(.r) = p>. This implies that the 
mapping must be local, i.e. the variation of the reference field at locations <p ^ C has no 
influence on X('$(£),t). It follows then from $(£) = H being the unique inverse image of 

$(x, t) = that 


( 


a 2 $ 

dx Q dx ( 


■|$(*) = <p) = m 2 


3<i> a$ a 2 * 


(4.15) 


holds. The correlations of the Gaussian reference maesure can be 


evaluated as follows 


and 


(^ 7 - — 1^(0 = '/) = ( — —) 


aco aca - 


' aCa dCc 


T / cW dty \ 

= r>) = 


1 ac^aco 1 


<* 2 > 


(4.16) 


(4.17) 


The variance {'I' 2 ) of the reference measure can be set to unity without restricting the gen- 
erality of the result. This completes the closure of the mapping equation. The resulting 
equation determining the mapping is thus given by 


dX 


= Q{i l) + m 2 r( 


a$ a* 

aCa d( a 


){X" - T}X') 
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This result was first obtained by Chen et al. (1989) and generalized by Pope (1991) to 
several variables using the fact that the .V- variate pdf can be represented as the product of 
N conditional pdfs. 

Method II.: A general transformation of the domain of definition of the pdf equation 
(4.2) is considered. The domain of definition is the space 5 = 3ft x [0. T } which is the image 
space of the mapping X_ : R l X [0, T] — » S (T denotes the time interval considered and )ft is 
the range of the values of 4>). The mapping is then written as 

t = X t (r),T), p> = X*{r],T) (4.18) 

The mapping must possess a unique inverse satisfying the same smoothness conditions as A_ 

r = A7 l (sM), *? = AV(vM) (4.19) 

and it must have a positive Jacobian. The mapping equations are established in two steps: 
First the pdf equation in the reference domain is set up and then the mapping condition that 
this pdf is Gaussian is introduced. The function f T defined by 

f T {r], t) = f(X<t>(r], r), X t {rj, r)) (4.20) 

is the pdf in the reference variables, but it is not pdf with respect to r/ since it it is not 
necessarily normalized. The transport equation for f T requires the transformation of the 
derivatives in (4.2) given by 


d_ _ dAg 1 d dX~ l d 

dt dt drj dt dr 

and 

d dAV d dX~ l d 

d’^> dip dr] dip dr 

The inverse transformation can be expressed in terms of the mapping itself (see C'ourant 
(1968), vol.II, ch.III, sect.4) and we obtain 


d 1 / dX$ d dX# d \ 

dt J \ dr dr] drj dr ) 

and 

d_ _ 1 / dX t d dX t d \ 

dp> J \ dr dt] dp> dr ) 

where J denotes the Jacobian defined by 

= dX^dXt_ _ dX* dXt 

dr dr] dr] dr 


(4.21) 


(4.22) 


(4.23) 
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The pdf equation in the reference variables follows now in the form 


dx* 0f T dX, 


df T ,dX t d 
+ {■ 


dX t d 

dr dr/ dr) dr 1 dr] dr dr dr] 


}(Fj/ T ) = o 


(4.24) 


where Fj = F^(X${r), r),Xt(r], t ) ). The mapping X$,X t is now in part fixed by the require- 
ment that the pdf with respect to the reference variables is Gaussian 


dX 


* r T 


dr] 


f iV X) = fo(r]) 


(4.25) 


or with (4.20) 


/(v,0 


fain) 

dr) 


where <p = X<t,(T),r) holds, 
obtain 


The derivatives of f T can be expressed in terms of fa and we 


df T fg [ dX* d 2 X* 

dr] (^^) 2 ? dr) dr] 2 


(4.26) 


and 


df_ _ [g &X* 

dr r d_X± 2 dndr 

' dt) 


(4.27) 


It follows that the Gaussian cancels out of the pdf equation which can be recast as 


dX* &X* dX* d \ ( dX$ T dX t_\ dX* d f T dX t \ T d^&X* 

^ dr] dr ] 2 dr] dr] ) \ dr * dr ) dr] dr \ * dr] ) * dr] dr\dr 

(4.28) 


If we set (compare to (4.10)) 

dX$ _ 'pdXt 

~dT ~ 

it then follows from the pdf equation that 


(4.29) 


d_ 

dr 



(4.30) 


must also be satisfied. There are now two distinct possibilities for the mapping: We can 
restrict the temporal map to Xt(rj,r) — Xt(r) which then implies that dXt/dr] = 0 and the 
second mapping equation (4.30) is trivially satisfied. Hence, the condition that the reference 
measure is Gaussian does not determine the scale factor dX t /dr. In the second case where 
the dependence of X t on r] is retained we find that the second mapping equation (4.30) 
determines the temporal evolution of Xt but not the variation with respect to r] which is set 
arbitrarily by the initial condition. 
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4.2 Mapping method for the characteristic function. 

The characteristic function corresponding to the pdf /(^>, t) is defined as the Fourier transform 


= — i= f difif(?,x,t)exp(ifO (4.31) 

\/27r J 

If we regard the pdf as the image of a Gaussian via = X(r 7 , t) we can transform the integral 
and obtain 

= ~i= f drj^-f(X(r],t),x,t)exp{iCX(rj)) 

v27r J otj 

It follows from the properties of the mapping X(r]) that 


m(C,x,t) = J drjf G {r])exp(i(iX(r])) (4.32) 

holds. If we introduce the characteristic function of the Gaussian m G ( 0 we find ( hat the 
characteristic functions are related by 






dwm G (uj)exp{i(£X(r}) - u ;/?)} 


(4.33) 


and we can regard 

Y(C,uJ) = -J== J dT]exp{i(CX{T],t) - ur))} (4.34) 

as mapping function. It is straightforward to derive the transport equation for 1 from (4.34) 
and the equation for X given by (4.10) 


— ( ^' 0 = m 2 (V* • V 4 ){-ui 2 r(C,u;) + ll(uiy)+ 



d{«fr(c,K)y(c,{)r(-c.« 


* — 0} + *c 



-w)y(c,«) 


(4.35) 


where Q is the Fourier transform of the source term Q. This equation is apparently nonlinear 
and its usefulness depends essentially of the way the limit of zero fluctuations is approached. It 
follows from the definition of Y (4.34) that the Gaussian characteristic function is approached 
if the mapping Y approaches the Dirac-function 


F(C,u;,<)-+6(C-u>) as *(7?) -+ >7 

This property eliminates Y from consideration since it leads back to the awkward properties 
of the time inverse diffusion equation whose solution must reduce the width of the solution 
profile as time evolves. However, a simple modification of the definition of the mapping 

Z(C,/?,t) = -^=exp(iC(X(j 7 ,<) - rj) (4.36) 
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leads to the limit 


2(C ,nJ) 


\/2k 


, as / ->o o 


and the corresponding transport equation appears in the form 


dZ r\i ssr, 2^,5* w d 2 Z 1 

^ = tQ(r/)CZ + m 2 r{— . . 

at d(, a dc, a drj 2 i(,Z\di]J 


( d J \ 


2 • -<7 

-*r&Z - i?-r-} 


C?77 


[4.31 


We note the absence of integral terms but the presence of nonlinear terms. The characteristic 
function can be recovered from 


m(C; x, t) 


J diomci^) 


J drjZ(CTj,x,t)exp{i(C ~ uj)tj) 


(4.38) 


In conclusion we note that there is considerable freedom in setting up a mapping procedure. 
The usefulness depends on two aspects: The approach of zero fluctuations and the number 
of independent variables, which becomes critical if more than a single probabilistic variable 
in the pdf is considered. 


5.0 Mapping method for the multi-dimensional pdfs. 

The extension of mapping methods to multi-dimensional pdfs requires some preparations. 
We begin with a fundamental property of the pdf equation and its consequence for mapping 
methods. The pdf equation is regarded as a first order pde and a slight modification of the 
characteristic theory of this class of equations (see Courant and Hilbert, vol.II (1962), ch.II) 
leads to the basic result. It will be shown that it is not necessary to resort to the cdf equation 
(as done by Pope, 1991) for the development of mapping methods. 
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5.1 Fundamental mapping equations. 

Suppose the pdf = Prob{p t <$,<>?, + dp> t , i — 1(1) A/ } depends on M 

probabilistic variables, i.e. f m integrates to unity with respect to , <pxi* and time. The 

transport equation for f is then given by 

=^‘. •■■)/, h]=° (5.1) 

1=1 

where the fluxes R t are subject to the conditions that random variables assume the values 
ipi for i = 1(1 )M . The particular structure of the conditional expectations {R, \ • • •) is not 
important for the following assertion. Consider now a local mapping X_ : R V? —* where 
denotes the range of the variables for i — 1(1 )M, 

9. = -V,( 77! , ■ ■ -,r)M ,t), i = 1(1 )M 

such that the Jacobian J defined by 


j = det (W >0 

remains positive. Let /g(^i, • • • ,rj m) be M-v ariate Gaussian and let 


(5.2) 


f M ( V 1 ( ^7 1 1 ' ' ’ i H Mi t ) i ' ' ' i V \{ ( 1 ' ' ' ' i HM i t) 


fG(Vir‘- RUf) 

J(r) i,- • • 


(5.3) 


be the pdf defined by fa and the Jacobian .7. We will prove that 
equation 



satisfies the transport 


(5.4) 


where <fi = X,(r} 1 ,---,rj M ,t). 

Proof: The time rate of change of /*j for 77 kept constant follows at once from implicit 
differentiation 



where </?, = X;(? 7 i, • • • , was used. Differentiation of the right hand side of (5.3) leads to 


Using (5.3) we get 
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The time rate of change of the Jacobian is given as a sum of determinants by 


d.J 

dt 



dXi 

. . d.\ k . 

d\ vf 

M 

dr) i 

dm 

dm 

E 


* 


k= 1 

d\\ , . 

. a A*. . . 

ax M 


di)M 

df)M 

dj)\t 


where the notation 


A', = 


dX, 

dt 


was used for the time derivative. Consider now a differentiable function F(X i, 
where the arguments are functions Xi(rj i, • • • , 17 ,v/ , f ) and set up 


■ A.u 


M 


dF _ dF dx k 

drjt ^ 


k = 1 


dX k drji 


Cramer’s rule leads to an expression for the derivatives with respect to the A, 

j*L 

ax, 

where the derivatives of F appe 
i = 1(1 )M produces 


dXi 

dF 

BX M 

dVi 

dm 

dm 

d\\ 

BF 

BXm 

df)M 


df)M 

in the ith column. 

Setting . 




dX t 

d\ k 

dXxt 

,ydX, 
^ dX, 

M 

dr) i 

dr) i 

dr) i 

= E 


• 

• 

1=1 * 

k=\ 

dXx 

dt)M 

ax k 

dr)M 

dX M 

dr)M 


which is identical with the time rate of change of the Jacobian. It follows that 


/ 1 dJ \ dX i 

\7~3t)rh' §T > 


holds. Combining the results for the left and the right hand sides and using 

<P, = X t (Tj u - ■ ■ 

leads to the conclusion that /^, as defined by ( 5 . 3 ), satisfies 

df\f 


M 


M 




dt ) ^ dip, 

/ * 1=1 r 


1=1 


dip. 
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as claimed. 

Note that the particular properties of the Gaussian reference measure did not enter the 
proof, only its time independence was used. It follows that any other time independent 
reference measures such as the measure with beta-function density, suitable for bounded 
scalars, could be used. Comparison of (5.4) with (5.1) shows that this result allows the set 
up of the mapping equations for any number of variables. It follows that /v/ = f\i holds if 


= 1*1 = <**!,•••), *' = l(l)Af (5.5) 

and the initial and boundary conditions for (5.1) and (5.4) axe the same. The relations 
(5.5) are the central result for mapping methods. It is instructive to compare (5.5) with the 
dynamic equations for the scalars t). The scalars axe governed by 

d$ 

-~(x,t) = J?i(£,^i(x,<),---,$A/(x,t),t) 

for i = 1(1) A/ and the right hand sides R, do not depend on the parameters v?i, • • ■ . y\Y/- 
Note that the dynamic equations may be taken at different points in the flow field and the 
location vectors x_ are then labelled accordingly. The mapping equations (5.5) contain the 
conditional expectation of the same right hand sides Ri but the expectations depend on the 
conditioning parameters v?i, • • • The dependence on the location is now parametric if 

all scalars are taken at the same point. 

The generalisation of (5.3) to time dependent reference measures 


(5.3') 

1) ' ' ' > Wd) 

leads to a modified equation for /^. The dependence of the reference measure on time is 
established for the case of a non-degenerate A/-variate Gaussian given by 


? • • • , t)i • • • t Xm(v i — 


fair 

J(n 


M M 

= {(2n) M det(M,j)}-% exp{-- ~ - Pj(t))} 

“ .=i j=i 


in terms of the time dependence of its mean value vector //(<) and covariance matrix 
The equation for the pdf defined by (5.3') can be shown to be 



— fli 


3 log if a) 


(5.4') 


dt 



I 


The derivative for the logarithm of the reference density /g can be established if the particular 
form of f G is known. It follows for the Gaussian that 

d log? /g ) _ dlogjfo) dp, ,sr"ST dlogt/cldA/.j 1 

at a,., at aw,-' a 

holds. It follows from the pdf equation (5.1) that a time dependent reference measure is 
inappropriate if the right hand side of the pdf equation is zero. 

The properties of the flux equations (5.5) depend essentially on the formulation of the 
basic laws. Their particular structure will be analyzed in the following sections for the spatial 
and material frames. The case of two-point pdfs will then receive special attention to illustrate 
the properties of mapping methods that are able to produce scale information. 


5.2 Mapping method in the spatial frame. 

The evaluation of the conditional expectations requires now the knowledge of the the partic- 
ular properties of the fluxes R,. The random variables <$, are now regarded as the the values 
of stochastic fields at one or more than one points x} 3 ^ in the flow field D. At each point 
— 1(1 )N in D a set of I\ variables $ j , j = 1(1) A’ consisting of velocity, scalars and 
other variables is taken as the probabilistic variables in the pdf. The notation is modified 
to indicate the location in the flow field. Accordingly are the fluxes and the 

components of the mapping denoted by R^'\ The fluxes R^ can be split into a local 

and integral contribution 





a**," 




d$ k 


(5.6) 


where the lack of a superscript in the integral contribution indicates that it depends on any 
location in the flow field. It should be noted that the presence of spatial derivatives in the R } 
required the extension to stochastic fields. Examples for the local and integral contributions 
can be found by inspection of the basic laws (3.1) to (3.3) and (3.8). It is easy to see that 



= v 


dx^dx^ 


is local and 


N j' ] = " 7 “ / d y~u)G(x {,) ,y) 

1 4 tt J - g x [ ■> 

D 3 


dv a dvQ 

'-d^dy Q 


is integral and, therefore, nonlocal. The Gaussian random variables representing the argu- 
ments of the mapping are also regarded as the values of stochastic fileds at N locations £ • 

i = 1(1 )N in the domain of definition R 3 of the fields j = 1(1)A. The argument fields 
(&)•(£), j = 1(1) A' are homogeneous Gaussian fields with time independent statistical prop- 
erties. The extension of the random variables to stochastic fields taken at A points in the 
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respective domains of definition implies that there must be a relation of these N points in 
R 3 for the Gaussian argument fields to the corresponding .V points in the flow domain D for 
the image fields. This mapping V : [£)] ;V —* [i? 3 ]‘ v is denoted by 


£ 10 (5.7) 

It is time dependent and was defined as pure stretching in the case of single point pdfs (see 
(4.12)). It is important to notice that Y. i s not determined at this point and the subsequent 

development will show that Y_ may depend on the mapping A’j*\ 

The fundamental requirement of mapping methods is now that the conditional expec- 
tations of the turbulent fields $0(x^*\f) is equal to the expectations of the images of the 
Gaussian argument fields i = 1(1)A, j = 1(1) *V . If we denote the fluxes taken at 

the image fields with R 


R 


(*) 




(1) (< ) , 


•), 




■u 1 


(1) 




(5.8) 


and 


$ ( .° = xj ,) (^ < 1 1) (r (1) (x (1) , • • • , x} N \ t)), ■ ■ • , 4' ( / ^ ) (F <iV) (x (1) , • • • ,x (N \ t )), x (1 \ • • • , £. {i \ t) 

3 3 (5.9) 

for i = 1(1 )A\ j = 1(1) A', we can write the mapping closure as 


<i?y’i«- < 1 1> ■ ■ •) o.io) 

where the fact that the mapping has a positive Jacobian was used to express the conditions 
on the image variables in terms of the argument variables. The equations 


dX) 

~df 


(*) 


-tol 


( 1 ) 


^A V) - 


.( 1 ) 






= 


- J 1 * 

— r}\ , 


•), j = 1(1) A*, 


i = 1(1 bV 
(5.11] 


(with appropriate change of notation compared to (5.5)) are called flux equations. Introducing 
the the representation (5.6) for the fluxes leads to 


dX) 


(«) 


f<‘) 




= 




(5.12) 


The presence of nonlocal contributions iVj** needs some attention, because they contain the 
values of the turbulent fields at locations y ^ for all i = 1(1 )iV, where they are not 
image of a Gaussian argument field since the mapping (being local i.e. only defined for the 
= l(l)JV) is not defined. Extending the mapping to nonlocal (or functional) character 
does not make sense since no tractable mapping equation would emerge. The only avenue 
open on the level of A -point pdfs is the construction of an additional closure for the nonlocal 
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terms (such a closure was outlined by Chen et al. ( 1989)). This aspect of mapping methods 
will not be discussed in the present chapter. 

It is important to notice that the pdf defined by (5.3) is not solution of the pdf 
transport equation (3.19) for N > 1 or non-homogeneous turbulence if the mapping Xj 
is only applied to the probabilistic variables $(* . This is a consequence of the fact that 
the locations x} 1 ^ in the spatial frame are parameters and not probabilistic variables since 
there is no transport equation governing them. However, it is clear from the structure of 
the convective term in (3.19) that the notion of the mapping can be extended to include 
Z_ : [i? 3 ]‘ v — ► [£)]^ which is determined by 


dZ 


(«) 


(«)/„(!) 


dt 




-r. 


id 


.(AO 


,0 


15.13) 


for j = 1,2,3 and i = l(l)iV, where the supscripts of A’j** are arranged such that j = 1,2.3 
correpond to the velocity vector. The closure is completed by requiring that Z is the inverse 
map of Y_ introduced in (5.7) 



(5.14) 


It is now apparent that we must require that Y_ : [,D]' V — ♦ [f? 3 ] lV has a unique inverse. It 
follows then that the mapping Y_ for the domains of definition is in general a function of the 
same set of independent varaiables as the mapping A °f the range of values. 


5.3 Mapping method in the material frame. 

The basic laws in the material frame can be given in mixed (spatial-material) notation as 
follows 

■X a {a,t) = V n (a,t) 


dt 


as kinematic condition and 


dVg 

dt 

d$ 


dVg 

dX a 

1 dP 


= 0 


+ V 


d 2 V a 


p dx Q dXfidXff 

d 2 $ 


dt Qi * ) + F dXpdXp 


(5.15) 

(5.16) 

(5.17) 

(5.18) 


where A(a, t), F(a, t), P(a,t) and $(a, t) denote now position, velocity, pressure and scalar 
in the material frame as function of the label variable a = A(a, 0) and time t. The time 
derivative is now the substantial derivative with the label a kept fixed. It was shown in 
section 1.1 that the implicit derivatives with respect to the actual location can be expressed 
in terms of derivatives with respect to the independent variable a, but the resulting expressions 
are highly nonlinear and will be invoked only if necessary. The set of dependent variables 
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consists now of position X, velocity V and scalar $ and the mapping equations developed in 
section 5.1 can be applied without modifications. It is worth noting that the minimal set of 
dependent variables appearing in the spatial frame (section 5.2) requires an extension of the 
mapping and results in the same set of variables for the mapping as in the material frame. The 
components of the mapping AT* 1 ' (not to be confused with location X in the material frame) 

axe ordered such that A’j 1 corresponds to location for j = 1,2,3, to velocity for j =4.5.6 
and to the scalar for j = 7 and the superscript indicates the material point a (,) € D( 0). The 
flow domain D{t ) is a function of time and represents the range of the dependent variable 
position X whereas D{ 0) is the domain of definition of ail dependent variables X, V and 
The random variables <$, appearing in the pdf equation (5.1) are now regarded as the 
the values of stochastic fields at one or more than one labels aS 1 ^ € -D(O). At each label 
aS'\i = 1(1 )N in £)(0) a set of K variables = 1(1)A' consisting of position i, velocity 

V_ and the scalar $ is taken as the set of probabilistic variables in the pdf. The notation is 
modified as in section 5.2 to to indicate the location in the initial flow field. 

The fluxes and the components of the mapping are denoted by R y *\ Xj 0 as before. The 
fluxes can be split into a local and integral contribution 


R 


(0 




dX'a" ' 


d$k 
' dX a ' 


(5.19) 


where the lack of a superscript in the integral contribution indicates that it depends on any 
location in the flow field ( mixed spatial- material notation is used for the fluxes). The Gaussian 
random variables representing the arguments of the mapping are also regarded as the values 
of stochastic fields at time t and at j V labels * = 1( l)iV in the domain of definition R 3 
of the fields j = 1(1) A'. The argument fields ^(C, 0, j — 1(1) A’ are homogeneous and 
stationary Gaussian fields. The mapping F introduced in the spatial frame by (5.7) appears 
naturally in the material frame. Recalling that the mapping X j 1 * represents for j = 1.2.3 
the position of the material point ( i ) 


• • • , £ (1) , • ■ • , a (N \t) 

we note that at time t = 0 

.Yj'W, •••, *'«»;*<», ..•,«'">,(»« af ( 5 . 20 ) 

must hold. The condition of a positive Jacobian for the mapping Xj ^ implies that there 
exists a unique inverse which is for j = 1,2,3 the Gaussian distributed argument position 
field at 

*,■( c (,) ) = • ■ • . 

(j = 1,2,3). For t = 0 we set 




0 ) 


(5.21) 
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and (5.20) implies that {Xj'*} *(• • • , 0) depends only on the label a (, \ Hence, we have found 
a mapping 



[5.221 


defined by 


jfVlspry’rv. 


(1)1-1/ (1) 




,£ ( ‘ V) ,0) 


[5.23) 


which is independent of time. In summary we note that ail argument fields are stationary and 
homogeneous Gaussian fields. The deterministic conditions (5.20) and (5.21) at time zero are 
enforced as initial conditions for the mapping Xj^ which approaches appropriate constants 
for t — ► 0, thus producing marginal Dirac pdfs for position. There is no need to introduce 
time dependent reference measures. 

The fundamental requirement of mapping methods is now that the conditional expec- 
tations of the turbulent fields is equal to the expectations of the images of the 

Gaussian argument fields * = 1(1 )/v, j = l(l)iV. If we denote the fluxes taken at 

the image (surrogate) fields with 


R 5° = • • • ,a iN) ,t) (5.24) 


and 


for i = l(l)iV, j = 1(1)A\ we can write the mapping closure as 

= tiV-> 


o.2o 


(5.26) 


where the fact that the mapping has a positive Jacobian was used to express the conditions 
on the image variables in terms of the argument variables. The equations 


dX) 


(0 


, V* \ « (1 \ • • • , « (N \ 0 = <^‘ ) ^ ( » 1) = ^ • • •>, J = KV)A', i = l(l)N 


^ (5.27) 

are called flux equations as in section 5.2. Introducing the the representation (5.6) for the 
fluxes leads to 


^ i*! 11 =»iv-> 


(5.28) 


The presence of nonlocal contributions ivj 1 ^ presents the same difficulty as in the spatial 
frame since they contain the values of the turbulent fields at labels a ^ <£*** for al 1 i — 1( 1)A . 
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5.4 Mapping method for two-point pdfs. 

The mapping methods developed in the previous sections will be applied to a special case of 
particular importance. We consider the pdf taken at two points in homogeneous turbulence 
and restrict attention to position and velocity desribed in the material frame. The pdf 
equation (3.19) appears now in the form 


9 -t +2: ■ ->1^> - 0 

(5.29) 

where the subscript indicates that two points are considered. Homogeneity in physical space 
implies that f*i depends on r e x — x} x ^ but not on x = an< ^ ^ follows that 

the pdf equation is reduced to 


i=l 


(5.30) 

The pdf fi depends on relative position r, velocity at two points v^'\ i — 1,2 and time t. 
The image space S of the mapping Jfj 1 * is, therefore, spanned by the range of these variables. 
Since no boundaries are present for homogeneous flows it follows that the image space is given 
by S — R 9 . The argument variables are all Gaussian random variables and this implies that 
the domain of definition of the mapping Xj 1 * is also given by R 9 . Hence, X } : R — * R ■ 
We denote the mapping of the relative position with Z_ : R 9 — *■ f2 3 and the mapping onto the 
velocity space with A ( ^ : R 9 -* R 3 , i = 1, 2. The image variables are now regarded as values 
of stochastic fields at two labels a (l) , i = 1,2. 

*(,) _ f Xj(s£ 3 \t) = AA :>(«<*>, 

* } { Vj(a^,t), i=1 ’ 2 

and likewise for the argument variables 


' \ i=1 ’ 2 

The mappings relate the argument and image variables 

AI(a (1| ,a (2) ,f) = Z(A£,£ (1 \* (2) ;a (1 \a (2 \*) 

and 

V(a {i \t) = X (0 (A£,£ (1) ,iL (2) ;a (1) ,a (2) ,t) 

where the label a in physical space and the domain of definition of the argument fields are 
related by the time independent mapping Y_ given by (5.22). The flux equations follow at 
once from section 5.3 in the form 


5.31) 

5.32) 


= A' (2) - A' (1) , j = 1,2,3 
dt 3 3 


(5.33) 
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and 


dX) 

dt 


(.) 


= (IflA*! = ,„•••) +( ; Vj"|A*, = ! = 1.2 


CO, 


The implicit form of the local and nonlocal parts of the conditional fluxes follow from 
as 

d 2 X‘„“ 


=!?„•••) = (» 


dZjidZa 


|A«, = Vu---) 


and 


1 d 


TO” ia*, = 

P a 


(5.34) 
(5.17) 

(5.35) 

(5.36) 


where the pressure is a functional of position and velocity according to (3.8) which can 
be translated into the material frame without difficulty. However, position and velocity in 
P\X_{-, t), V(., f)1 cannot be expressed in terms of any argument fields because the mapping 
is defined for two labels only. 

The evaluation of the local part of the flux equation is rather complicated and will be 
outlined without explicit calculation of all expressions. Furthermore, it will be assumed that 
the mappings do not depend parametrically on the labels i — 1,2 and the dependence on 
the labels is via the map Y_ only (the maps Z_ ans X}^ are autonomous with respect to label). 
The position field is dependent variable and the mapping is defined for two different labels. 
This implies that the implicit derivatives with respect to actual position must be expressed 
in terms of derivatives in terms of labels. The Laplacian is thus according to section 1.1 given 
by 


d 2 x { a ,] _i dz^dz* d ( dz c dz* dxl ' ] . 

dZjdZfi ~ 2 dar, da^ da/ dap da-, da 9 


[5.37 


The argument fields A.)P and = 1,2 depend on the label a via the map i_(a) and 

parametrically as indicated in (5.31) and (5.32). It is now more convenient to expand (5.37) 
into 


d 2 X 


(*) 


1 


dZ < dz+ dZ < dZ+ d 2 X 


Ai) 

Of 


dZ^dZ* dZ c dxl l) d 2 Z 


dZpdZp 2 da,, da M dap da da 9 das 


u ; 


da v da u 


dap da 9 


da y das 
(5.38) 


Implicit differentiation leads to 

d_Za_ = dZ^d^dYy dZ^d^dYy dZ^d^dYy ( - 3Q) 

dap d<ps < 9( 7 dap d<p^ ^ap 

where Z_(g_, £)* The mappings Z_ and X } * * were assumed to be local. This implies 

that the derivatives dZ a /dps* • • as well as dY^/dap (since it is the inverse of Z at time 
zero) are completely determined by the conditions A$ = ^7**' an d the conditional expec- 
tation applies to the derivatives of the Gaussian argument fields only. The second 

derivatives follow from repeated application of (5.39). The conditional expectation of the 
Laplacian emerges as complicated combination of products of first and second derivatives of 
the mappings Z and X}*\ The degree of nonlinearity is clearly given by the nonlinearity of 
the momentum balance (1.5). 
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6.0 Conclusions. 

Several conclusions can be drawn from the mapping methods discussed in the present paper. 

1. Mapping methods applied to velocity pdfs require additional closure assumptions 
since the pressure depends in functional (integral) form on velocity. 

2. Random variables which are regarded as values of stochastic fields at particular 
locations require an additional mapping for these locations in the domains of definition of 
image (turbulent) and argument (Gaussian) fields. 

3. The computational effort for the flux equations is approximately A/-times the effort 
for the pdf equation where M is the number of probabilistic variables of the pdf. Mapping 
methods are, therefore, not competitive for M > 1. 
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1.0 Introduction. 

The prediction of compressible turbulent reacting flows requires careful consideration of all 
possible formulations of the basic laws to produce a set of equations suitable for pdf methods. 
The first step is, therefore, devoted to the study of possible formulations of the pdf method 
for compressible turbulent flows with combustion reactions. This aspect of the research work 
was carried out in a previous grant and is documented in the report by Farshchi et ah ( 1991 ). 
The second step is the development of closure models for this type of flow with particular 
emphasis on the effect of compressibility. The pdf method can be based on the transport 
equation for the pdf of thermo-chemical scalars plus variables measuring the rate of relative 
volume expansion or the material derivative of the pressure (see Farshchi et ah, 1991). The 
pdf approach offers the possibility of treating chemical non-equilibrium in rigorous fashion, 
which is particularly important for high speed flows characterized by high shearing rates and 
short residence times. The progress achieved in the development of a closure model for the 
pdf equation valid in this situation and the successful application of this model to supersonic 
hydrogen flames will be discussed in detail. 


2.0 Pdf transport equation for compressible flows. 

The single point pdf equation for scalar variables determining the local thermodynamic state 
is considered. Turbulent flow at supersonic speed can be modified significantly by compress- 
ibility and the interaction with shocks created outside the turbulent flow field and random 
shocks (shocklets, Johnson et al., 1973) generated in supersonic turbulent shear layers. Pdf 
methods can be adapted to cope with the effects of compressibility including random dis- 
continuities and combustion. We consider the case of infinitely fast reactions, in which three 
variables determine the local state: Mixture fraction, pressure and enthalpy or any other 
equivalent set of thermodynamic variables. Pressure can vary significantly in supersonic 
flows and enthalpy is not conserved due to frictional heating in high shear regions. Hence, no 
further simplification, as in the case of low Mach number subsonic flames, is possible. The 
single point pdf f\ is then set up for the velocity £, density p, or a local function of density 
such as log(p) which will be used below, internal energy u, relative rate of volume expansion 
D and mixture fraction £ (a choice that was found to be advatageous by Farshchi et ah, 1991) 


/i(u, d, = (<5(v — v)Hp — d)S{e — u)6{D — C)<5(£ — //)) 


T 


The transport equation for this pdf can be obtained using standard methods and emerges in 
the form 

9f ' 3 {-(£./) + i<§^/> + Bd(fJ)) 


UJ\ , UJ In 

«( ~ + v a ) — — o l \ * 

at OX a OV a OX a 
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Mean thermodynamic properties follow from the pdf /i by integration using the local relation 
determined by equilibrium considerations. The mean pressure for instance is given by 


(P) 



<%>((/, u . i) )/i ( d. u,Tj ) 


where p(d , a, p) denotes the local relation of pressure to density, internal energy and mixture 
fraction. The calculation of this type of local relation is straightforward. The relations for 
other thermodynamic variables such as composition and temperature to the pdf variables 
density, internal energy and mixture fraction were established using the equilibrium code 
STANJAN ( Reynolds, *1986). 


2.2 Closure model for the pdf equation. 

The pdf equation to be considered is the result of integration of (2) over velocity space. It is 
given in terms of the density-weighted pdf / 1 defined by (Kollmann, 1990) 


? _ P(V i. 

/i = — 


PH 


{ P ) 


/l ( P 1 


rP-Lit) 


where p corresponds to the scalar variables (p, e. D, (,') and l — 4. The integrated pdf transport 
equation for the set of / thermo- chemical variables follows from (2) for high Reynolds numbers 
in the form 


(P){ 


dfi 

dt 


+ ~ ^jUKh)} + 

dx 3 v. a o u c 


- 1 


«C/i + 7(7 - 1 


Re 


<*/>} 


OTj OX 3 OX a OX a OX a p OX a 
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■«p)K = 


d 7 ,dg a a. _ d_ J_ 
OuPe { dx Q J) 


_(1 dT ?£)f) - — { — (pT^-)f) 

d( Re ' dx a p dx 3 ' dp dx Q dx a 


14) 


The terms on the right hand side can be shown to contain dominant terms describing turbulent 
mixing in scalar space which has the well known structure of a time-inverse diffusion process 


° 4 ) -EE 


di 


1 1 


d 2 


j=2 S 3 d ^ k 


((e J k\* J =Y})fi) 


[ 5 ) 


and the scalar dissipation rates f tJ in the conditional expectations are defined b\ 


e 0 


dx a dx a 


( 6 ) 


with equal diffusivities T, = Tj = T for simplicity. Note that no such term acts in the first 
scalar direction which corresponds to the variation of density. 
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2.3 Mixing Model. 

Any closure model for the mixing process described by (5) should share as many properties as 
possible with the exact term. It should preserve normalisation and mean values and decrease 
variances and covariances. The pdf should remain nonnegative and should not spread outside 
the domain of realizable states. The pair interaction model for the / — 1 (note that one of the 
scalars does not mix as it corresponds to the density, which does not diffuse) scalar variables 
is defined by 


= -{ / dp' f dp"f l (p')f l (p")T(p'.p".p)-f 1(£)} (7) 

It is assumed that all scalars are appropriately normalized such that the scalar space (set of 
all realizable states) 1? is a subset of an / — 1-dimensional unit cube. It should be noted that R 
may have intricate boundaries as a consequence of realizability conditions that mass fractions 
cannot become negative or exceed unity and that they must add to unity. The transition pelf 
T(p' ,p" ,p) must satisfy the requirements 

2V.£".£) = T ( £'.£".£'+£" -£) (8) 



and 

r(£',£".£) = 0 for p 2 AV,£") (0) 

The central part of the condition (9) is the construction of the neighbourhood A(r ) 
which is the interval [p\p"\ in the single scalar case. N can be at most the cube O-i = 
{ip : i E [y>', y 1 '/],/ = 2,/} defined by p' and p" for pairwise interaction according to our 

assumption of normalisation of the scalars. Realizability requires that the mixed states are 
in 3ft, hence 

V c C,., n ft (10) 


must hold. Symmetry 

£ 6 N <**£' + £"-£ E N (IK 

must be imposed to insure the properties of a mixing model. Furthermore is T pdf with 
respect to p 

f d£r(£ , .£".£) = l < 12 > 

Conditions ( 8 )-( 12 ) do not define the mixing model uniquely but represent a class of models. 
It is important to realize that the structure of the scalar domain -ft modifies the neighbourhood 
N unless N is reduced to the line connecting £' and £ /f . If pf and £" are close to the bouudan 
of the neighbourhood is essentially the connecting line due to (10), but if the pob aie 

inside 3? and far away from its boundary then may .V be the cube C/_i- The transit. pdf 

T determines the particular form of the mixing model and is set up in the present case as 


T(£',£".£) = G(Otf(£'<£".£) 


13) 
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where 




1 


for £6 .V( £\ £ ") 
ot herwise 


( 14 ) 


and m-i(X) is the /— 1-dimensional volume of „Y(y\ y”) and C denotes the centered variable 


= 




+ /')]• t=2./ 


It follows that T satisfies 


and 


/ v <KG{ O 


dCG(0 = 2 / - 1 


G(0 = G(-0 


(13) 


16) 
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The present choice for the function Cr(0 is a constant determined by the condition ( 16) which 
can be regarded as the condition to assign equal probability to all possible outcomes of the 
mixing interaction of two fluid elements. The mixing model for (5) is thus set up. 

2.4 Compressibility Effects. 

It is advantageous to set up the closure model representing the effects of compressibility 
in the Lagrangean frame as stochastic differential equations. The basic laws are written in 
abbreviated form 

d log/? 


for mass balance 


dt 

dD 

dt 


= -D 


= Qd 


IS) 


19) 


for the balance equation for the relative rate of volume expansion 

dE 


dt 


= Qi 


20) 


for the energy equation in terms of the internal energy per unit mass and finally 

ti-Q. 

dt ~ ^ 


21 ) 


for the mixture fraction, where d/dt denotes the material derivative and D = V e the relative 
rate of volume expansion. The right hand side terms are conveniently set up in the Eulerian 
frame (which can be considered implicit Lagrangean expressions). The basic laws combined 
with the constitutive relations for Newtonian fluids lead to the explicit form of the Q t given 

by 

Q = J d 1 0t q3 d 1 dp dv Q dv 3 dg a (2 2) 

D ~ Re dx a p dx j dx^ p dx a dxj dx a dx a 


o 



where q Q denotes the energy flux vector. 


Q e = -*,(7 - l)Mt(l+p)D + -^7(7 - DM'ft ~ 


7 d f ln 
Re Pr d.r n 


23) 


and where $ is the dissipation function. 



_2 o oc_ 

Re Sc d r a P 0x a 


(24) 


The general form of the closure model (24) for the pdf equation given above (4) (which con- 
tains the dynamics of the variables density, internal energy, relative rate of volume expansion 
and mixture fraction) is set up using the form of stochastic differential equation 


clY x = A t dt + btjdWj + dJ t (*-°) 

The stochastic nature of (25) is given by dW t , which is the increment of a normalized random 
process (such as the Wiener process), and dj ln which is the increment corresponding to a 
jump process. The closure for the equations (19) - (21) will discussed for each of them in 
some detail. 

I. The time rate of change of the relative rate of volume expansion consists of three 
contributions: The increment due molecular transport which is regarded as mixing, the in- 
crement due to the passage of isentropic compression and expansion waves past the material 
point considered, and the passage of random shock waves past the material point considered. 

The first contribution i\D mtT is represented by the mixing model (7) together with 
the requirement that T is constant in its domain of definition. It can be shown that the 
viscous term in (22) implies indeed that D is subject to diffusion. Hence will D participate 
in the mixing model described in the previous chapter. The second contribution is modelled 
according to an Ornstein-Uhlenbeck process 

A D IS = {c p if{M a )—^-}*ri - c p2 f(M a )^(D - (D)) (20 

where c p \ = 1.0 and c p 2 = 0.5 are constants (the values given here are arbitrarily chosen and 
a systematic variation is discussed in a later section) and 

f{M a ) = M* 

is an empirical function of the local Mach-number. It ensures that the increment of D vanishes 
as the Mach-number goes to zero. The first part of A D is a W iener process ( 7 / is a Gaussian 
random variable with zero mean and unit variance) representing the random stirring effect of 
isentropic compression and expansion waves moving past the material point considered. The 
second part is a drift term ensuring the existence of a steady state. Finally, we note that r is 
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the turbulent time scale provided by the second order closure (see Dibble et al.. 1986). The 
model for Q o has now the form 


Q £> Af= AD, n IX + ‘AD,, 4 - AD,* 


where the first and the second contributions have been established. The last term represents 
the random occurrence of shocks. This contribution is nearly singular and corresponds to the 
derivative of a Dirac-pseudofunction in the inviscid limit. There is no model for it at present 
and a way of treating random shocks will be discussed in the section II below. 

II. Mass balance (IS) does not require closure and contains only a drift term 


d log p = —D clt 


(28) 


as long as the relative rate of volume expansion remains sufficiently smooth. The case of 
random shocks leads to a singularity for D and will be treated as separate contribution to 
d\og(p)jdt in the form of a jump process. If the local Mach-number is greater than unity, 
shocks may appear with the maximal strength given by the normal shock relation 


G(M a ) 


V V- 

T+^aTJ 


for M a > 1 
otherwise. 


(29) 


and the increment d.J\ for the jump process representing the random shocks is modelled by 


d.J x =G(M a )jN s {-)ii (30) 

k r 

where N 3 (6) denotes a nonegative integer random variable representing the number of shocks 
arriving at the material point in 6 dimensionless time units and 0 < i] < 1 is the random 
variable giving the dimensionless shock strength. The current model for A .,(£>) a Poisson 
process and r) is a random variable with uniform distribution. The complete increment for 
the logarithm of density is thus given by 

A log p == — DAt -f dJ\ + A Ddts ( 31 ) 

where the last contribution is due to frictional heating at constant pressure. This contribution 
is given by 

A D^s = piC « + Au,p) - p(C P) 

where Au is the increment of internal energy due to frictional heating. Finally, we note that 
logp does not participate in the mixing process. 

III. The increment for the internal energy consists of several contributions 
Q E At=Au miI + Au„ + {- 7(7 - l)M*(l+p)D + -^ 7(9 (32) 
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The first part is due to heat conduction and is therefore part of the mixing model applied 
to internal energy. The second part of the increment is due to the isentropic expansion and 
compression waves passing the point considered and can be written as 


— m(C. p 4- Ap) - »/((,'•■*< p) 


where Ap = —Di\t denotes the change of density as a result of the change in the relative 
rate of volume expansion. The third part contains the pressure work term and the frictional 
heating contribution. The dissipation function consists of 


= u 





for flows of boundary layer type. The last contribution is due to the random arrival of shocks 
at the material point considered. 


S 



3.0 Prediction of supersonic hydrogen flames. 

The closure model developed in the previous chapter was applied to the prediction of super- 
sonic hydrogen flames burning in coflowing stream of air. The flow configuration was a round 
H 2 jet with a coflowing stream of air at higher temperature than the fuel. The prediction 
requires accurate initial data which will be discussed next. 

3.1 Initial conditions. 

The flow conditions of the first test case of Evans et al. ( 1978) are shown in Fig.l. 



Fig.l Flow geometry for the supersonic H ^- air flame. 

Cold hydrogen at T h, = 251 A’ is injected at the axis of a circular supersonic air flow 
generated by a convergent -divergent nozzle. The air temperature is given by T air = 1495A . 
This temperature is achieved by burning hydrogen upstream of the nozzle and then adding 
oxygen to the hot products to produce A*o 2 = 0.21 mole fraction corresponding to air. The 
air flow contains, therefore, a high percentage of water as product of the heating process 
(Xh 2 o — 0.281). The boundary conditions and the nozzle geometry are summarized in Table 
1 . 
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Exit condition 


H 2 -Jet Outer Jet 


Mach Number 

Ma 

\ 

2.0 

1.9 

Temperatur, K 

T 

251.0 

1495.0 

Mean Velocity, m/s 

U m 

2432 

1510 

Pressure, MPa 

P 

0.1 

0.1 

Mass Fraction 

Yh 7 

1.000 

0.000 


Yo 7 

0.000 

0.241 


Kv, 

0.000 

0.478 


Yh 7 o 

0.000 

0.281 


Table 1. Initial data. 



The calculation of the turbulent nonpremixed flame is carried out with the hybrid method 
developed by Chen and Kollmann (1988). The first step in the solution procedure is the 
calculation of the thermo-chemical properties, which are stored in a table for the later use in 
the solution of the pdf equation. The results of this calculation were reported by Farshchi et 
al. (1991). The next step is the set up of the initial (or entrance) conditions appropriate for 
the first test case of Evans et a. (1978). The initial velocity profile is shown in Fig. 2 where 
the symbols indicate the initial location of the grid points. The H 2 stream contains 30 grid 
points and the coflowing air stream 18 points. 



Fig. 2 Initial velocity profile. 
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4,0 Results and comparison with experiments. 

The thermo-chemical properties of the reacting mixture of H 2 and air are determined from 
the condition of chemical equilibrium constrained with pressure and internal energy. This 
assumption is unrealistic for many situations in supersonic flows, but it is the logical first 
step in the development of prediction models for such flows. The extension to chemical non- 
equilibrium has been carried out for zero Mach-number flames (see Chen and Kollmann. 1 OSS, 
1990) and, once the questions concerning compressibility effects on the turbulence structure 
have been sorted out, the results obtained for zero Mach-number flames can be applied to 
supersonic flows. 

The fuel considered in the present prediction was a mixture of hydrogen and nitrogen 
(Yh2 = 0.22335 and 1> 2 = 0.77665, in order to raise the stoichiometric value of mixture 
fraction from — 0.0283 for pure hydrogen fuel to £ 3t = 0.113. The pure hydrogen case was 
also considered but only results for the former case will be presented. 

The pdf equation is solved using a stochastic simulation technique (see Pope, 1985 for 
details) together with a second order closure model for the first and second order moments of 
the velocity field (Dibble et ah, 1986) including modifications accounting for compressibility 
effects (Zeman, 1989) on the dissiaption rate based on direct simulation results by Lele ( 1989). 


4.1 Parametric study of the compressibility model (26). 

The compressibility model (26) represents the random passage of compression and expansion 
waves passing the material point considered. The function f(M a ) = is an ad hoc model for 
the unknown dependence of this model on the Mach-number. It is clear from the consideration 
of the low Mach-number limit that this function must be nonnegative and vanish as the 
Mach-number approaches zero. These two properties are obviously satisfied by this function. 
The model (26) contains furthermore two constants c pl and c p2 for which no information is 
available at this time. Hence two reference values c p i = i.o and c p2 = 0.5 were chosen and a 
systematic variation of the constants was carried out to learn how the solution depends on 
them. The results for one of two sets of runs are presented in Fig. 3 to Fig. 16 varying the 
constant c pl . There is only a limited amount of experimental data available in Evans et al. 
(1978), which consist of Pitot pressure measurements and some composition information at 
two cross sections. The Pitot pressure results at the first station x/D = 13.8 in Fig. 3 shows 
cleaxly that increasing the value of c p i leads to significant improvement. It should be noted, 
however, that perfect agreement is not to be expected near the axis where the temperature of 
the fuel is below the minimal temperature of T = 296 A for which thermodynamic data were 
available and the equilibrium relations could be established. The initial temeperature of the 
fuel stream had to be set to this temperature and not the temperature of the experiments. 
The profiles for the mean velocity (which is calculated using a second order closure model 
solved parallel to the stochastic simulation procedure for the pdf equation) in F0g.4, the mean 
density in Fig. 5, the mean temperature in Fig. 6, the mean internal energy in Fig. 7, the mean 
mixture fraction in Fig. 8 and the mean value for the relative rate of volume expansion D in 
Fig. 9 show the corresponding variation with c p i. The mean value for D in Fig. 9 is apparently 
zero, with some numerical noise which is typical for stochastic simulation techniques, \isible. 
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™ ai sta * ion at ,/d = ^ «>■*“ <>* ^ 
in Fig 10 is Quite en l §reement between the calculated and the measured Pitot pressure 

temperature ?n Fie n ^ ! be aX1S for reasons gained above. Only the mean 

mean dilation ' V' i 6 mean ^nsrty in Fl g* 12 show a strong dependence on c pl The 

the mean stress occmrine the ‘° ‘“““'‘i 1 * * hs * here is “° turning „f 

Fiff 18 shows rP ki g ‘ e mean com P 0!,1 Gon for fuel and oxidiser in Fig. 17 and 

tig- 18 shows reasonable agreement between calculation and measurements. 

4.2 Mean values and correlations. 

^FklTto'trJ 0 ^/ SCa J ar , fiel ' ls and ,h<f velocity field are presented at v/D = 10.5 
variables of th/n If ' fm ' dh ;’ ds allow ,lre ‘'alculation of any moment of the probabilistic 

7e2 Zu s t P ,t ,n d 7 der ° f ,he is <»'>■ btnited by numerical accuracy The 

the thermodynamte variables show that the mean internal energy I Fig 20) 
s m mrnal at the ax.s s.nce the enthalpy of formation for the fuel is negative. Mean delitv 

for 1<w Mach. 

loTaTon h o e ft S L Pe fl r niC flame aS Can , be Seeu m Fi S-2ri P smTp;essu'e Ip^sl^tTe 

pressure The mlT 18 ap .f lrent J lnd 4 the Pressure in the cold fuel is higher than the ambient 

difference has he "' T, ^ ^ 24 Sh ° WS that at x/D = 155 most of the initial velocity 
t en smoothed out and only small mean strain rates are present . The Reynolds 

the J° mPOnen S m Fl ?'“ 5 indlcate the presence of two shear layers between the fuel jet and 

n YfZl T Jet , an k the air Jet a “ d Stagnant — ' ding - The mean dis,p aim “ ‘e 

or^fa L 0 7 h ° " f imier kyer f ° rmed bet — fuel and coflowing air jets. The 

ZtlZ therm ° dyn T C VaHableS haW — 1 “^cresting features. Mixture fraction 

fTet s ne7a ive e b S i T ^^elated (Fig.28) since the internal energy for the 

h^el is negative but he mixture fraction for the fuel is maximal (unity). The correlations of 

mixture fraction and density (Fig. 29) and internal energy and density (Fig.32) change sign 

small andTatV^ 11011 Fhe ^ rrelatlons with the relative rate of volume expansion are rather 
small and rather noisy Fig.30, 33, 35 and 36). 

4.3 One-dimensional pdfs. 

-7 7- rr pd r f ° r tbe thermodynamic scalars is contained in fig.37 
Fiff 41 {u i *° °l ^ Ve ra< ^ a ^ stat ions. The pdf for mixture fraction in Fig.37 to 

inteial ° f Slgn ° f the skeWneSS aS the flow is traversed. The pdf of the 

opposite^kiiff lg i-~ t0 Flg has a sha P e similar to the pdf of mixture fract ion with the 

of entrained he^ A S " eW ' nes ®' < be pdf for densit y in Fig.47 to Fig.51 shows the appearance 
ted air at radia l stations greater than .r/D = S4 (Fig.49) which emerges as 

ch!se to the Ga ~ ° 23 rY ? <lf ° f the relative rate of volume expansion in Fig. 52 to Fig.56 is 
process " ^ t0 the m ° del eqUati ° n (26) Emulating an Ornstein-Uhlenbeck 


4.4 Two-dimensional pdfs. 



The effect of compressibility becomes apparent if two-dimensional pdfs are considered. The 
comparison of the one-dimensional pdfs for mixture fraction and internal energy (or any 
other thermodynamic variable except enthalpy, which is a linear function of mixture fraction 
at M a ~ 0.0, and mixture fraction) does not lead to unambigous conclusions, because the 
local relation between those variables, that holds at zero Mach-number, is nonlinear. The 
Ornstein-Uhlenbeck process described in chapter 2.4 as model for the random fluctuations of 
the relative rate of volume expansion leads to a broadening of the pdf for thermodynamic 
variables and mixture fraction, which would be related locally in incompressible flows. The 
pdf of mixture fraction and internal energy shows some broadening due to compressibility as 
Fig. 57 to Fig. 61 prove. However, the pdf of density and mixture fraction in Fig. 62 to Fig. 66 
exhibits a much more pronounced broadening in particular in shear layer between fuel and 
heated air at x j D = 0.27 in Fig. 62. The pdf of density and internal energy in Fig. 67 to Fig. 71 
confirms this fact. The pdf containing the relative rate of volume expansion D as one of the 
variables allow some insight into the properties of the model suggested in equation (26). The 
pdfs for internal energy and D in Fig. 72 to Fig. 76 and in particular the pdfs for density and 
d in Fig. 77 to Fig. 81 show that The statistics of those variables are not Gaussian but the 
marginal pdf for D is close to Gaussian. 

5.0 Conclusions. 

It was shown that pdfs for three scalar variables describing the local thermodynamic state 
in a compressible reacting flow can be determined as solutions of model equation that sim- 
ulates the effects of convection, turbulent diffusion, chemical reactions and reversible and 
irreversible compression and expansion processes occuring randomly in a turbulent flow at 
high speed. The limited amount of experimental information does not allow to draw a final 
conclusion concerning the accuracy of the calculations, but it is clear that pdf predictions 
of compressible reacting flows are feasible. There are, however, several problems awaiting 
solution. In particular the role of the fluctuating pressure containing several different modes 
(acoustic mode, entropy mode) and the significance of chemical non-equilibrium need to be 
investigated. Pdf methods are especially w^ell suited for the latter because they allow rigorous 
treatment of nonlinear and local processes. 
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Fig. 3. Coaxial turbulent supersonic jet flame burning H 2 with air. Parametric study of the 
compressibility t<ini foi the stochastic simulation of the relative rate of volume expansion 

(equation (26)). Mean Pitot pressure at ., / D - 13.8. Open circles are the measurements of 
Evans et al. (1978). 
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Fig. 4. Coaxial t urlmlenf supe*r.semie* jet Hame burning fF with air. Parame'tric study of the 
compressibility term for the stochastic simulation of the* relative* rate* of volume expansion 
(equation (20)). Me*an velocity at x/D = 13. S. 
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Fig. o. Coaxial turbulent supersonic jet Hume burning H_> with air. Parametric study of the 
compressibility term for the stochastic simulation of the relative rate of volume expansion 
(equation (20)). Mean density at .r/D = 13.8. 
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Fig. 6. Coaxial turbulent supersonic jet flame burning with air. Parametric study of the 
compressibility term for the stochastic simulation of the relative rate ot volume expansion 
(equation (2G)). Mean temperature at x/D — 13. S. 
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Fig. 8. Coaxial turbulent supersonic jet tiame l)uniing H > with ail. Parametric stu<ly ol the 
compn ibility term for the stochastic simulation of the relative rate* of volume expansion 
(equation (2G)). Mean mixture fraction at xjD — 13.8. 
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Fig. 9. Coaxial turbulent supersonic jet flame burning H> with air. Parametric study of the 
conipiossibilitv term for the stochastic simulation of the relative rate of volume expansion 
(equation (2G)). Mean relative rate of volume expansion at r/D ~ 13.S. 



Fig. 10. Coaxial turbulent supersonic jet Maine hurtling H-» with air. Parametric study of the 
compressibility term for the stochastic simulation of the relative rate of volume expansion 
(equation (20)). Mean Pitot pressure at x/D — 20.2. Open circles are the measurements ol 
Evans et al. ( 197S). 
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Fig. 11. Coaxial turbulent supersonic jet flame burning H i with ail . Parametric stud\ of the 
compressibility term for the stochastic simulation of tin* relative rate of volume expansion 
(equation (20)). Mean velocity at xjD = 20.2. 
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(equation (2G)). Mean density at ./• / D — 2G.'j 
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Fig. 13. Coaxial turbulent supersonic jet flame burning H-> with air. Paramet ric study of t 
compressibility term for the stochastic simulation of the relative rate of volume expansion 
(equation (2G)). Mean temperature at x/D — 2G.2. 
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Fig. 15. Coaxial turbulent supersonic jet Haine burning with air. Parametric .study of tin* 
compressibility term for tin 1 stochastic simulation of tin* relative rate of volume expansion 
(equation (26)). Mean mixture fraction at x / D — 2G.2. 
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Fig. 1G. Coaxial turbulent supersonic jet Hame burning tF with air. Parametric study of the 
compressibility term for the stochastic simulation of the relative rate of volume expansion 
(equation (26)). Mean relative' rate of volume expansion at x/D = 26.2. 
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'ig. 17. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean mass fractions at 
/D — 15.5 for ( jj | = 1.0 and < f ,> = 0.5 compand to the experiments of Evans ct al. ( 197S). 
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Fig. 21. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean density at 
x/D = 15.5 for c p 1 = 1.0 and c p 2 = 0.5. 



Fig. 22. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean teniperatuie at 
x/D = 15.5 for c p i = 1.0 and c p 2 = 0.5. 




Fig. 23. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean pressure at 
x/D = 15.5 for Cp\ = 1.0 and c p2 = 0.5. 
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Fig. 24. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean velocity at 
x/D = 15.5 for c pl = 1.0 and c p2 = 0.5. 








Fig. 25. Coaxial turbulent supersonic jet flame burning H 2 with air. Reynolds stress com- 
ponents at x/D = 15.5 for c pl = 1.0 and c p2 = 0.5. 



Fig. 26. Coaxial turbulent supersonic jet flame burning H 2 with air. Mean Dissipation late 
at x/D — 15.5 for c pl = 1.0 and c p2 = 0.5. 
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Fig. 27. Coaxial turbulent supersonic jet flame burning H 2 with air. Variance of mixture 
fraction at x/D = 13. 3 for c p \ = 1.0 and e p i = 0.5. 



Fig. 28. Coaxial turbulent supersonic jet flame burning H 2 with air. Covariance of mixture 
fraction and internal energy at x/D = 13.3 for c pl — 1.0 and c pi = 0.3. 






Fig. 29. Coaxial turbulent supersonic jet flame burning H 2 with air. Covariance of mixture 
fraction and density at x j D = 15.5 for c p \ — 1.0 and c p 2 = 0.5. 



Fig. 30. Coaxial turbulent supersonic jet flame burning H? with air. Covariance of mixture 
fraction and relative rate of volume expansion at x/D = 15.5 for c p \ = 1.0 and r f) > = 0.5. 
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Fig. 31. Coaxial turbulent supersonic jet flame burning H 2 with air. Variance of internal 
energy at x/D = 13.3 for c pX = 1.0 and c p2 = 0.5. 



Fig. 32. Coaxial turbulent supersonic jet flame burning H 2 with air. Covariance of internal 
energy and density at x/D = 13.3 for c pl = 1.0 and c pl = 0.3. 





Fig. 33. Coaxial turbulent supersonic jet flame burning H 2 with air. Covariance of internal 
energy and relative rate of volume expansion at xjD — 15.5 for c p t = 1.0 and c oi = 0-3. 



Fig. 34. Coaxial turbulent supersonic jet flame burning H 2 with air. Variance of density at 
x ! D = 15.5 for c p x = 1.0 ancl c /yi = 0.5. 
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Fig. 35. Coaxial tubulent supersonic jet flame burning H 2 with air. Covariance of density 
and relative rate of volume expansion at x/D = 15.5 for c p \ = 1.0 and c p i — 0.5. 
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Fig. 36. Coaxial turbulent supersonic .jet flame burning with air. Variance of relative rate 
of volume expansion at xjD — 13.3 for c fjl = 1.0 and c ^ — 0.3. 








Fig. 37. Coaxial turbulent supersonic jet flame burning H -2 with air. Pdf of mixture fraction 
at x/D = 13.3 and r/D — 0.27 for c p \ = 1.0 and c p i = 0.5. 



Fig. 38. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of mixture fraction 
at x/D - 13.3 and r/D = 0.33 for c {)[ = 1.0 and c /t , - 0.3. 





Fig. 39. Coaxial turbulent supersonic jet flame burning Hj with air. Pdf of mixture fraction 
at x/D — 15.5 and v/D = 0.S4 for c p i = 1.0 and c f) i = 0.5. 



Fig. 40. Coaxial turbulent supersonic jet flame burning Ho with air. Pdf of mixture fraction 
at x/D — 15.5 and r/D = 1.14 for r /A = 1.0 and v p -> = 0.5. 
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Fig. 41. Coaxial turbulent supersonic jet flame burning H 2 with air. Pelf of mixture fraction 
at x/D = 13.5 and r/D = 1.42 for c lA = 1.0 and c p > - 0.5. 
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Fig. 42. Coaxial turbulent supersonic jet fianie burning H> with air. Pdf of internal energy 
at x/D = 15.5 and r/D = 0.27 for r n , = 1.0 and c.., = 0.5. 
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Fig. 43. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energy 
at x/D = 13.5 and r/D = 0.35 for c p \ = 1.0 and r,, 2 — 0.5. 



Fig. 44. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energy 
at x/D = 15.5 and r/D = 0.84 for c f){ = 1.0 and c p2 = 0.5. 










Fig. 45. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energy 
at x/D = 15.5 and r/D = 1.14 for c pl = 1.0 and c p2 = 0.5. 
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Fig. 46. Coaxial turbulent supersonic jet tiame burning H 2 with air. Pdf of internal energy 
at x/D = 15.5 and rf D — 1.42 for c pl = 1.0 and c p2 = 0.5. 
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Fig. 47. Coaxial turbulent supersonic jet flame burning H* with air. Pdf of density at 
x/D = 15.5 and r/D = 0.27 for c p \ = 1.0 and c p i = 0.5. 
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Fig. 51. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of density at 
x/D = 15.5 and r/D = 1.42 for c pl = 1.0 and c p2 = 0.5. 



Fig. 52. Coaxial turbulent supersonic jet flame burning H> with air. Pdf of relative rate of 
volume expansion at x/D = 15.5 and r/D = 0.27 for r,„ = 1.0 and r„, = 0.5. 







Fig. 53. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of relative rate of 
volume expansion at x / D = 15.5 and r/D = 0.55 for c p \ = 1.0 and c p 2 = 0.5. 



Fig. 54. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of relative rate <>t 
volume expansion at v/D = 15.5 and r/D = 0.S4 for e pl = 1.0 and c f ,i = 0.5. 






Fig. 55. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of relative rate of 
volume expansion at r/D = 15.5 and r/D = 1.14 for c p 1 = 1.0 and c p 2 = 0.5. 





Fig. 56. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of relative rate of 
volume expansion at r/D = 15.5 and r/D = 1.42 for c pi = 1.0 and c p 2 = 0.5. 










Fig. 60. Coaxial turbulent supersonic jet flame burning H -2 with air. Pdf of mixtvue fraction 
and internal energy at r/D = 15.5 and r/D = 1.14 for c {) \ = 1.0 and c p > = 0.5. 






Fig. 62. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of mixture fraction 
and density at .r/D = 13 . 0 and r/D = 0.27 for c p \ =1.0 and c a i = 0.5. 
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Fig. 69. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energ} 
and density at x/D = 15.5 and r/D = 0.S4 for r pl = 1.0 and c p2 = 0.5. 
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Fig. 71. Coaxial turbulent supersonic jet flame burning H 2 with air. Pelf of internal energy 
and density at x/D = 13.3 and r/D = 1.42 for = 1.0 and c p 2 = 0.3. 





Fig. 72. Coaxial turbulent supersonic jet flame burning H -2 with air. Pdf of internal energy 
and relative rate of volume expansion at v/D — 15.3 and r/D = 0.27 for c p i = 1.0 and 







Fig. 73. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energy 
and relative rate of volume expansion at x/D = 15.5 ancl r/D = 0.55 for c p 1 = 1.0 and 
c p2 = 0.5. 



Fig. 74. Coaxial turbulent supersonic jet flame burning H? with air. Pdf of internal energy 
and relative rate of volume expansion at r/D = 13.3 and r/D = 0.S4 for c p) =1.0 and 






Fig. 75. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of internal energy, 
and relative rate of volume expansion at x/D = 15.5 and r/D — 1.14 for c p \ = 1.0 and 
c P 2 = 0.5. 
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Fig. 77. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of density and 
relative rate of volume expansion at x / D = 13.3 and r / D = 0.27 for c p i = 1.0 and c p > — 0. 
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Fig. 78. Coaxial turbulent supersonic jet flame burning H 2 with air. 'Pdf of density and 
relative rate of volume expansion at x/D = 15.5 and r/D = 0.55 for c p \ = 1.0 and c p 2 = 0.5. 



Fig. 79. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of density and 
relative rate of volume expansion at x/D = 15.5 and r/D = 0.84 for c p \ = 1.0 and c p i = 0.5. 












Fig. 81. Coaxial turbulent supersonic jet flame burning H 2 with air. Pdf of density and 
relative rate of volume expansion at x/D = 15. 5 and rjD — 1.42 for c pi = 1.0 and c pi = 0.5. 
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